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The literature contains two forms of the two-nucleon potential computed to fourth order in the 
coupling constant from the gradient coupling of the pi-meson to the nucleon field. These are com- 
monly referred to as the Taketani-Machida-Onuma (T.M.O.) and the Brueckner-Watson (B. W.) 
potentials. The merits of the controversy surrounding this schism are reexamined from first principles 
starting from the covariant equation for two-nucleons, and it is concluded that the conditions for the 
applicability of the method leading to the T.M. O. potential are never satisfied in practice. On the 
other hand, the B. W. potential, suitably altered following recent suggestions by Miyazawa and the 
author may well yield a suitable approximation in the low energy region. 


SI. Intreduction 


Recently the author’? and Miyazawa” have suggested a new approach to the problem 
of deducing from meson theory the interaction between two nucleons. In particular they 
have demonstrated how within a framework, rigorous except for the neglect of a possible 
meson-meson coupling, the asymptotic form of that interaction for large internucleon 
separations—usually termed the potential—may be tied completely and directly to para- 
meters deduced from pion-nucleon scattering. 

In view of the immodest strength of the above assertion, it may well appear con- 
tradictory to have to admit that the functions proposed as “potential” by the two sets 
of authors differ in an important way. Moreover the nature of that disagreement, 
usually referred to as the difference between the T. M.O.” and the B. W.” potential, 
represents the essense of the dispute of longest standing among the cognescenti of this 
subject. The main objective of this note is to attempt to shed definitive light on this 
controversy. It was felt that the most satisfactory way to proceed. would be to reformulate 
fundamental notions in the soundest manner available, and thus to be able to examine 


critically alternative and prior attempts. Many of our arguments have been given previous- 


ly, though we can claim novelty for several observations and for the attempted coherence 


of the account. 
We take as starting point the assertion that the most satisfactory method that we 


know of in principle of studying the properties of a two-nucleon system is to deduce 


k : ’ ere 
them from the solutions of the covariant equation for two-fermions 
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Here & (x,x,) is the sixteen component amplitude dependent on a pair of space-time 
-coordinates. Since there is a spinor index for every space-time coordinate, the former 
will be suppressed throughout. G,(x,, */), i=1, 2, are the Feynman or causal propagators 
for a physical nucleon and thus comprehend fully the self-field of each particle. The 
‘exact field theoretical definition of the quantities in eq. (1) as well as a convenient 
derivation of the exact form of the interaction I will be described in a succeeding paper. 
If we restrict ourselves to interactions via the pion field above, I is nevertheless represented 
by an infinite series which is a dual expansion both in the number of mesons exchanged 
and in the non-linear behavior of the mesons. In the applications so far contemplated, 
‘the non-linear aspect is neglected. On the other hand, the expansion of I contains in 
maximally summed form those self interactions which can be ascribed to an individual 
‘nucleon. 

The outstanding advantage of eq. (1) are that it can be renormalized and that the 
necessaty choice of an approximate form of I can be made so as not to destroy either 
the covariance or the renormalizability. It has been argued against this equation that 
‘the amplitude 7 (x,x,) does not admit of any simple statistical interpretation and that 
it does not lead to a standard boundary value problem, but neither objection can be 
maintained on principle. Against the former one, it is only necessary to exhibit formulae 
by which all observables can be related to the amplitude 7. This has been done by 
several authors." Against the latter objection, it need only be shown that the boundary 
‘value problem is well-defined, for example that a suitable normalization condition exists 
for the discernment of members of the discrete spectrum and that scattering solutions 
‘exist. Recent work®’' indicates that this is probably the case. 

Indeed the most serious objection against eq. (1) is our impotence to produce solu- 
‘tions for physically applicable situations. In view of this, various efforts have been made 
‘to replace eq. (1) by a more Schrodinger-like equation, containing in principle still com- 
plete information on the system in question, but one in which consequent approximations 
would result in a problem for which physically useful solutions could be obtained. To 
‘characterize these efforts, let us introduce total and relative coordinates 


X=4(x,+%),° x=%,—%, (2) 


‘and further suppose that we are dealing with an eigenfunction of total energy-momentum 
vector P,=(P, W), in particular that P=0. If we then define ¢/(x, x,) and (x), 
irespectively, by the equations 


P(x, %) Se AG (x, %), (3) 
9 (x) =9(x, 0), (4) 


we first ask whether eq. (1) can be replaced by an equivalent equation of the form 


[W —H,(p) —H,(—p)]x(x) = \ D(x, x’) 2(x!)d°x! (5) 
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where %(%) either equals (x) or is related to it in some known manner, 
i, (p) =[—ia V+ pM] (6) 


is the free-particle Dirac operator for the i” nucleon, and D(x, x’) is the full three- 
dimensional kernel. Equivalently if we think of 7(x), or rather of its Fourier transform, 
as decomposed into its constituent positive and negative energy free particle projections 
and all negative energy components somehow eliminated, we may consider in place of eq. 
(5), an expression of the form 


[W —E,(p) —E,(—p) |Z(x) = W—2E(p)) £(x) 


= \ Ux, x) 1a )id* x’, (7) 
where 
AG) ak 6 dia 200 ied Goel AL © 


and it is henceforth understood that by %(x) we mean a wave-function with positive 
energy components only. 

Many authors” have emphasized that a transition from eq. (1) to eq. (7) is always 
possible, where, moreover, U(x, x’) may be chosen to be independent of the total energy 


of the system. It is essential here to resume the common thread of the argument. 


§ If. Method based on scattering matrix 


Let the S matrix to which eq. (1) gives rise for the scattering of two nucleons 
from four-momenta p,, p, (and suitable other coordinates which are suppressed) to momenta 


pi’, po’ be represented in the form 


( pi’ po’ |S |p; po) =0(p, py 0 (po! — Po) 
2210 (P/—P;) 6(P'—P)T (pp, p; P,P), (2) 


where we have introduced the total and relative variables 


pi=sP+p, pi =sP' +p’, 


Po=eP—p, pl =sP’—p’. (10) 
In the center of mass frame, P=P’=0, we have 
P,=2E(p) =2E(p’). (11) 
We may therefore introduce a reduced T matrix t(p’, p) defined by 
T(p’, p; P=0, P,)=t(p’, Pp), (12) 


regarding the P, dependence as indistinct, according to (11). 
With 


u(p’, p) = (27) | exp(—ip!-»/+ip-») U(x, x')d?x d*x’, (13) 
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the corresponding T matrix associated with the solution of eq. (7) for a scattering situa~ 


tion is the solution of the well-known Lippman-Schwinger integral equation 


pout Ep che, be oe ee ee ee y 
PRD) MPA? | 2E(p) +ie —2E(p”) (14) 


If we now take t(p’, p) as given from the covariant eq. (1) by eqs. (9) and (12) 
and look upon eq. (14) as an integral equation for u(p’, p) then for given t(p’, p) 
that equation presumably provides a unique determination of the interaction kernel. An 
ambiguity which may be serious for bound state problems arises from the fact that eq. 
(9) determines T only up to expressions which vanish on the energy shell, which has a 
similar consequence for u(p’, p). We may at least assert, however, that there exist a 
class of interaction functions, non-local but energy independent, which are equivalent with 
respect to scattering consequences. In our future discussion, we shall suppose a definite 
form of u derived from a definite choice of T. In practice such a choice usually presents 
itself in a natural way as the one which emerges directly from the calculation that leads 
to (9) without any reexpression by the gratuitous use of the energy conservation condition.” 

In practice again ¢(p’, p) is given as some form of power series simply related to 
the kernel I of eq. (1). Equation (14) is itself then solved by iteration methods 
giving rise to a new power series whose convergence properties are a matter of fundamental 
interest. Though we do not wish to enter here into the further questions, it should also 
be pointed out that other observables of the theory, for example, electromagnetic moments, 


are represented by similar power series. 


§ If. Method based on approximate relative time dependence 


Some light is shed on questions of convergence by the discussion of an alternative 
mode of obtaining a three-dimensional equation. Its end result replaces eq. (7) by the 
form 


(W—2E(p)) o(x) = |V (x, x"; W) o(x') dx (15) 


in which the kernel V(x, x’; W) is in general a transcendental function of the total 
energy W and again a power series of some kind. The basic assumption of this method 


is a separability hypothesis, that the relative time amplitude ¢/(x, x,) may be expressed 
in the form 


Sin ee | (x, ati x; W) ete) deat (16) 


By employing eq. (16) on the right hand side of eq. (1), the latter itself may be used 
in an iterative procedure to determine J once the zeroth approximation has been stated. 
The starting point is taken to be the relative time dependence characteristic of a situation 


in which the interaction I is itself instantaneous, in which case the separability hypothesis. 
(16) is known to be fulfilled exactly.” 
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We shall omit further details since they have been recounted exhaustively in the 


literature.” The resulting equation (15) has however a pair of properties which are 


essential to the present discussion. For the first of these we consider the Fourier transform 
of V(x, x’; W) which we denote by v(p’, p; W). The associated transition amplitude 
t(p’, p; W), as determined from the analogue of (14), 


cr se t he ; W d° 


agrees for W=2E(p)=2E(p’) with t(p’, p) of eq. (12). From this point of view 
eq. (15) is fully as satisfactory as eq. (7). 

The second point concerns the relationship between eq. (7) and eq. (15). We 
here describe a procedure by which the transition between the two can be effected. We 


ask that eq. (15) be equivalent (in momentum space) to an equation of the form 


(W—2E(p)) Lp) =| ve(p, PY LCP) ap, (18) 


where %(p) is taken to be a normalized amplitude. We are also free to require for the 


amplitude of (15) 


|e) ¢(p)d*p=1, (19) 


though this necessarily implies a non-unit norm for the state vector from which @ was 


‘derived. Under these circumstances we can seek a unitary connection 
e(p) =| Alp, ps W)Z(p') d°p’ 
At (W) AWW) =1. (20) 
A(W) and v, are then determined by the matrix equation 
(p\ {di W) o(W) AW) +A W)2E AW) } |p’) 
—2E(p) 0(p—p') =w.(p, Pp’), (21) 
and by (20). To eliminate the explicit energy dependence of the left hand side of 


(21) we carry out the power series developments 


AAT 
Ai (W) =A (2E(p)) + W—-2E(p)) - aa paar ants? (22) 
AW) =A(2E(p')) + W—9E(p’)) - sa Prono Sr ae ' (23) 
| 
v(W) =v(E(p) +E(p’)) + W—-E(p) —E(p')) -- eeege ns miGaa) 


a ue E(p) + Ep") 


Inserting (22) —(24) into (21), and noting that the expectation value of (21) in 
the state %(p) is implied, we may replace W7—2E(p) on the left and W—2E(p’) on 
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the right by v,. We obtain an equation for the determination v,(p, p’) and the matrix 
elements of and its energy derivatives to be solved in conjunction with (20). In 
practice the power series method in the coupling constant or some small modification 
thereof is the only one available. For example, an equation which is correct to fourth 
order in the coupling constant (ignoring the terms involving E(p)) but which also 
implies selected higher order effects is given by the expression 

Ov(W) | 


ow E@)+E@n!) P’)- (25) 


ve(Ps P') =0(P, p's E(p) +E(p')) +(p —— jes 
An equation such as (25) will be compared below with eq. (24) which suggests itself 
as the expansion to be used directly in (15). 

To pursue the discussion further at this stage, we must ask for the connection between 
v.(p, p’) and the kernel u(p, p’) associated with eq. (7). Our general discussion 
should make it clear that there is no a priori reason that they be identical term by term. 
What is evident is that they again lead to the same scattering matrix on the energy 
shell. For the same reason that this is true they must also yield the same form of the 


adiabatic potential, to the discussion of which we must now turn. 


§IV. Comparison of methods and conclusions 


To our knowledge, no fruitful solutions of eqs. (7) or (15) have been obtained 
for a non-local kernel. Though this reality must eventually be faced, efforts to date have 
been directed toward replacing U(x, x’), for example, by a static potential, 


U(x, x') >U(x) 3(x—2’) (26) 


a replacement that is justifiable asymptotically and combining this with a phenomenological 
treatment of the interaction at small distances. In momentum space this corresponds to 
the reduction of u(p’, p) to a function v(p’—p) of the momentum transfer alone. 
This can be brought about and is justified whenever both p and p’ can be considered 
small compared to M and thus differences such as E(p)—E(p’) can be neglected." 
Since this same neglect is possible on the energy shell, where u(p’, p) and x,(p’, p) 
coincide, it follows that these kernels also give rise to the same adiabatic potential. If 
we further restrict ourselves to the coupling of p-wave mesons (gradient coupling) with 
nucleons, work only up to the exchange of two mesons, ignoring all radiative corrections 
other than those which produce the renormalized coupling constant from the bare one, 
the resulting local interaction constitutes what has been termed the T. M. O. potential. 

When we turn next to the corresponding specialization of the kernel (x, x’; W) 
of eq. (15), in addition to the requirement stated in the previous paragraph, necessary 
for the passage to the adiabatic-limit, two additional conditions must be met. These 
are that we may ignore the differences W7—2M and E( p) + E(p’) —2M, again compared 
to a typical meson energy. With these neglects, we find as before 


V(x, x's; W) >V (x) d(x—x’), (27) 


ina 
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Under the same coupling conditions as stated in the previous paragraph, the resulting 
local interaction is what has been termed the B. W. potential. 

We now compare the two potentials. We shall be forced to conclude that under 
conditions for which the adiabatic approximation is valid, if the T.M.O. and B.W. 
potentials differ in their predictions of the interaction in any state of the two-nucleon 
system, then if either is correct, it must be the B. W. potential. This conclusion follows 
from the considerations associated with eqs. (18) —(25). Thus the expansion (24) is 
justified only if the “adiabatic” hypothesis is warranted for V(x, x'; W). This may 
well not be the case for energetic collisions of two nucleons, as an example for 100 Mev 
available energy in the center of mass system. As a prototype of what occurs at low 
energies, however, Brueckner and Watson’ have estimated the relative size of the second 
term of eq. (24) compared to the first for the Deuteron problem and found it to be 
quite negligible. However, when treated as in the transition to eq. (25), where it gives 
rise in a perturbation treatment to a contribution to the fourth order potential (exchange 
of two mesons), it is well known to yield a large, indeed dominant contribution to the 
central potential in the triplet even state. In this transition, however, the second term 
of (24) is replaced by a power series in the coupling constant of which the fourth order 
term is only the first one. The series in question must be either poorly or non-convergent 
in the region of interest. This is illustrated for eq. (25) in the Appendix. We further 
conclude that the method of treating the two-nucleon problem based on an energy- 
independent kernel, V(x, x’) where the latter is understood as a power series is not 
applicable to pseudoscalar mesons in any practical sense. 

To complete the discussion we must endeavor to dispose of the main argument 
remaining which is adduced in the attempt to refute the above conclusions.” It is based 
on the historically oft-used adiabatic method for the treatment of the two nucleon problem.” 
In this method one computes the interaction energy of two nucleons at rest and uses 
this as the potential energy in a non-relativistic Schrodinger equation. Now if one applies 
this to the neutral scalar meson theory, it is well-known that the second order Yukawa 
potential indeed provides the exact static interaction energy. If one now turns to the 
methods of this paper it is found that U(x, x’) of eq. (7) does indeed have this value 
in the adiabatic limit, whereas V(x, x’; W) of eq. (15) gives rise to finite corrections 
of fourth and higher order. From this it is concluded that the adiabatic limit has not 
been properly taken in the latter case. 

We now insist that the argument given above begs the question. For it employs as 
a club the adiabatic method, whereas it is precisely the conventional application of the 
latter which is in contention. It should be evident by now that the question just raised 
can be settled only by appeal to non-adiabatic methods such as we have exposed in this 
paper. Alternatively phrased, there is no a priori reason that the exact answer of the 
adiabatic method yield the most reliable simple approximation to the fully relativistic 


situation. Indeed, we believe that we have established that for pseudoscalar meson theory, 


this is definitely not the case. 
In the work of Miyazawa’ and the author” alluded to in the introduction, all 
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results were obtained through the intermediary of an S-matrix formalism similar to that 
which leads to eq. (7). However, the work of the author was modified so as to take 
account of what was asserted to be the known form of eq. (15) as applied to the 
formalism of that paper. In the paper which follows, we verify those assertions, at the 
same time extending the previous formalism. 

Our general conclusion is therefore that eq. (7), with U (x, x’) represented by a 
power series in the number of mesons exchanged is never applicable. We have, of course, 
not thereby established the validity of eq. (15). Arguments given elsewhere," however, 
lead us to believe that with V(x, x’; W) determined by the method of the following 


ce 


paper, it should have a useful domain of application and that its adiabatic form’ may 


be valid at low energies. 

Of course, the use of (15) with its energy dependent interaction implies that the 
formula for other observables must be suitably modified, and the utility of the complete 
program based on (15) becomes tied to a rapid convergence of the series representing 
these quantities as well as that which describes the interaction. This latter question will 
be dealt with on a future occasion. 

The author takes this opportunity to recall with pleasure stimulating conversations 
on the topic of this note with Drs. K. Brueckner, D. Feldman, N. Fukuda, S. Machida, 
and M. Ruderman. 


Appendix 


We shall illustrate here how the procedure leading from (24) to (25) of the text 
gives rise to results which gainsay the fundamental assumptions on which it is based. 
Working only to the second order in the coupling constant for v(p, p’; W) and passing 
to the adiabatic limit, we easily obtain from (25) the following equation for the determi- 
nation of v,(r)™, 


ve (r) =%%(r) —ry' (r) u(r) (A-1) 
‘where 
Vo (r) =4(f?/47) Ut, T.(e7*/x) X [0+ ox? (224 3x4 3) Sie], (A-2) 
vy! (r) =h(f°/47) t,- t(2/7) {F,- F9[ (K, (x) /x) —K,(x)] 
—Sjol (2K, (x) /x) + K(x) ]}. (A-3) 


‘We have chosen units in which b=c=1, S\, is the tensor operator, K,(x) and K,(x) 
are the Hankel functions of imaginary argument and x=yr is the separation measured 
in units of the pion Compton wavelength. 


We shall solve (A-1) in the given approximation for the triplet even central force. 
‘We write in this state 


u(r) =G;\(r) +G, (7) Si, (A-4) 
and further define 
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A= — (f?/42) p(e-*/x), 
A, = (f?/4m) p(e-*/3*) (x24 3x43) 
a= (f?/47) (2/7) [ (Ki (2) /x) —Ky(x)] 
A= — (f?/47) (2/7) | (2K, (x) /*) + Ky(x) ]- (A:-5) 
By means of the relation 
S\5=6+20,-0,—28,, 
> 8—25, (A-6) 


for the triplet spin state, we easily obtain two linear equations for G, and G,, with 


solution 


A, ¢! =a, 2a,) — 8a, A, 


G, (r) ac ee 

(T=@.)\ (ia, 22a) S802 (A-7) 

(LIA ve 
1—a,+2a, 


For a typical value of (f*/47)—=.1, two observations can be made. [Even at the force 
range |G,(r)|>|A,(r)|, the latter measuring the strength of the second order central 
potential. Moreover G,(r) itself blows up at x~.7. The comparison of G(r) with 
A,(r) though not quite as catastrophic, is qualitatively similar with respect to the second 
observation. In any event this is sorry behavior for the quantity |G,(r) —A,(r)| which 
should be small compared to |A,(r) | itself and which apparently” has this property when 
treated in the form (24), the approximate passage to (25) being avoided completely 


because of its demonstrated invalidity. 
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Exact plane wave solutions of the field equations in general relativity and those of the field 
equations in non-symmetric unified field theory are obtained under similar conditions for the case in 
which electromagnetic fields are present. Main properties of both solutions are compared with each 
other. The most important result is that, as far as the solutions dealt with in this paper are con- 
cerned, the space-time is closely connected with the electromagnetic field in general relativity, while 
in non-symmetric unified field theory the structure of the Riemannian space-time is determined quite 
independently of the electromagnetic field. 


$1. Introductory remarks 


The object of this paper is to compare some wave solutions of field equations in 
the general theory of relativity with those of field equations in the non-symmetric unified 
field theory of Einstein, and to offer some materials for the solution of the problem 
which of the two theories is more appropriate for describing natural phenomena. These 
two theories will hereafter be referred to as GR and NS respectively. 

The field equations in NS coincide formally with the purely gravitational equation for 
empty region in GR when the non-symmetric fields (fundamental tensor and coefkcients 
of connection) are reduced to the symmetric ones, and the difference between both 
theories appears only when these fields are not symmetric. To express it in terms of | 
physics, this means that the both equations differ only in some region where there 
exists some other field besides gravitation. Therefore we shall confine ourselves to some 
region in which both gravitational and electromagnetic fields are present. 

It is not easy, however, to compare the solutions of both theories ‘completely ’. 
This difficulty comes mainly from the fact that the fundamental ideas with which both 
theories are constructed are not the same and the characters of both field equations are 
different from each other. Main differences are as follows: (1) In GR, the source 
terms for gravitational field are contained in the right-hand side of the field equations, 
but on the contrary, all such terms are absorbed in the left-hand side of the field equa- 
tions in NS. (2) In GR, it is generally accepted that some kinds of singularities re- 
present the existence of matter and that solutions containing such singularities are admitted. 
On the other hand, only those solutions that are regular everywhere are required in NS. 

In this paper, therefore, as the first step of the comparison, we shall compare both 
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solutions ‘in the small’. Thus, we restrict the problem in some limited portion of the 
space-time, and do not deal with the problem concerning the singularities of the solutions, 
nor the one whether we can extend the solutions in question into the whole space-time. 
In making such a comparison it is very desirable to deal with exact solutions. The 
both field equations, however, are considerably complicated, therefore it is not so easy to 
obtain their solutions, especially exact ones. Up to the present, many investigations have 
been done in GR concerning its field equations and various solutions obtained. However, 
if we exclude those for cosmological problems, there have been obtained few exact non- 
static solutions. Especially, little is known concerning exact wave solutions for the case 
in which electromagnetic fields are present. As to NS, the circumstances are more un- 
satisfactory, and it will not be too much to say that we have no knowledge concerning 
such a wave solution. This is the present situation of the exact wave solutions. 
Recently the present writer has succeeded in obtaining some exact solutions of the 
both field equations corresponding to the system where matter does not exist at all and 
the electromagnetic field is composed of plane waves propagating in one direction. The 
‘details of the calculation are given in the references.”*) In what follows, we shall ex- 
plain these solutions and then compare some main features of the both solutions in 


conformity with the considerations given above. 


§ 2. Field equations 


(a) Field equations in GR 
As was stated in § 1, we are going to deal with a system composed of purely radia- 
tional fields alone. Accordingly, not only the energy tensor due to mass but also the 


-charge-current vector must vanish everywhere. Thus, the field equations become 


G,,= —87E F¥, ,=0, 


aj? 
(GRE) 
Fiat etl g=0, CG, jy S1) %-) 4); 

where G,, is the Ricci tensor of the space-time, E,, is the electromagnetic energy tensor, 
F,, is the antisymmetric tensor representing the electromagnetic strength, and a semicolon 
followed by an index denotes covariant derivative. It is needless to say that the last 
two equations of (GRE) are the generalized Maxwell equations and are reduced to the 
usual Maxwell equations in the Minkowski space-time. 


(b) Field equations in NS 
The strong field equations in NS are as follows : 


Fig. n=O, I./=0, 
+= Vv j 
a (NSE) 

R,;=9, OF 1 ee a  eaaS 4) 

where g,,, 1’, r i i 
Isp Py eae R,; are the non-symmetric fundamental tensor, the coefficients of non- 

Symmetric connection and the generalized Ricci tensor respectively. The notations + and 
— under indices indicate the kinds of covariant derivatives. 


In (NSE), in view of later convenience, we have used the same notation ¥/,, to 
‘ a 
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denote the non-symmetric fundamental tensor as the symmetric fundamental tensor in 
GR. But there will be no apprehension of confusion. 
It should be noted that, in NS, the field equations are given by the above (NSE) 


irrespective of the kind of system under consideration. 


§ 3. Identifications of mathematical expressions with 
physical quantities 


-In GR, the physical quantities are well defined and their mathematical expressions 
well known. The gravitational potential is expressed by the fundamental tensor, and 
further, as was stated in last section, the electromagnetic strength, the electromagnetic energy 
and the charge-current vector are expressed respectively by F,,, Eg=9aF gk" /4—F pF 9” 
and Fi, 53 

In NS, on the contrary, there is no rule of the correspondence between the 
mathematical expressions and the physical quantities, which is universally accepted. With- 
out this rule, however, we have no way of continuing our reasoning. Therefore, follow- 
ing the opinions of some authors, e. g. Einstein,” Hlavaty,” Ikeda,” we shall make the 
assumptions : 

(i) h;==9:; plays the same role as g,;; in GR. That is, h,,; expresses the gravita- 
tional potential. 


(ii) The electromagnetic strength F,, is given by 


Bysavs =9 Sapa Gly 25 (3-1) 


where g=det.g;,, 7” is the tensor defined from g,; by 9,,9=0,, and €,,, is the relative 
tensor whose components are antisymmetric with respect to every pair of indices and 
€ 1934 1. 

(iii) The charge-current density vector is expressed by 


PHO On 1/33, (3-2) 


where a comma followed by an index denotes usual partial derivative and «”” is the 
relative tensor similar to the €; i. 

Here we shall give two remarks concerning these assumptions. 

(1) We start with the identification (3-1). However, even if we adopt the 


relation 
Fiys=9is (3.3) 
Vv 


in place of (3-1), we arrive at the same results so far as the solutions dealt with in 
this paper are concerned. Of course, we can not expect such a coincidence in more 
general cases. . 

(2) Strictly speaking, F,, defined by (3-1) is not an ordinary tensor but a idee 
tensor. In order to make Fj, an ordinary tensor, we have to multiply the right-hand 


member of (3-1) by a suitable factor. Detailed investigation concerning this fact has 
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§5. Exact solutions 


In this section we shall explain exact solutions which are obtained from Assump- 


tions 1 and 2. 
(a) Solution of (GRE) 


Under our two assumptions, the general solution of (GRE) is as follows : 
(GRS) gi; is given by (4-1), and F;, is given by 
(E: Fy=o, Fy=—p, Fu=0, 


Fg (5-1) 
H: Fy=p, Fy=¢, F,=0, 
where p, o as well as A, B, C, D are functions of Z. Six functions are required to 
satisfy only one condition, 

P= — 87 (Av? + 2Dv0 + Bo*) m, (5-2) 
where m= AB—D*, (5-3) 

P= (Av—2Dw-+ Bu) /m 

=m/2m—m 4m? —mC 2mC— (A B—D*) /2m, (5-4) 


( 2u=A— (BA? + AD*—2DAD) /2m—AC/C, 


20=B— (AB? + BD’—2DBD) /2m— BC/C, (5-5) 
2w=D—(BAD+ ABD—DAB—DD*) /2m—CD/C, 


and a bar over a kernel letter means derivative with respect to Z, e.g. A=dA/dZ, 
A=@?A/dZ?. 

We can consider the relation (5-2) as an algebraic equation to determine F,, (i. e. 
(, ©) when g,; (i.e. A, B, C, D) is given, or alternatively, as a differential equation 
to determine 9,; when F;; is given. Anyhow, it is only one condition for six functions 
and accordingly it admits solutions with great arbitrariness. Let (9:33 Ps ©) be a solu- 


tion of (5-2). Then, since (”, ©) is contained in (5-2) in the combined form 
a(p, 0) =Ay*+2Dp0-+Bo® (>0), (5-6) 


there exist other (”, o)’s which satisfy (5-2) together with the same g;;. Detailed 
mathematical research concerning such problems is given in the reference.” 
(b) Solution of (NSE) 

Using the assumptions stated in the last section and one additional technical assump- 
tion which is explained in the Appendix we can solve (NSE). The solution is as 


follows” : 


(NSS) bys =Ies and f;;==9,; are given respectively by (4-1) and 


fs=—fu=0'> fa=—fu=o', fio=fa=9, GG-7) : 


/ . : 
where p’ and o’ are arbitrary functions of Z, and A, B, C, D are subject to only one 
condition 
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P=0, (5-8) 
P being the quantity given by (5-4). (p', 0’) and (A, B, C, D) are quite independent 
of each other. 
For this solution, F,, is of the same form as that in the (GRS), 1. €. 


Fe=Fy=0, Fy=F,=0, —Fy=Fy= (349) 
where we have put 
p= (Bo! —Do"’) //m, o= (Ao!—Dp’) /Vm. (5-10) 
Further, if we calculate the charge-current density vector j' defined by (3-2), we have 
identically 


Here we shall add two remarks. (i) A Riemannian space-time satisfying (5-8) is 
not Minkowskian in general.” 


Gi) If we use the identification (3-3) in place of (3-1), we have 
Fye=Fy=0, Fy=Fy=—', —Fy=F,=0' (6712) 


i 


in place of (5-9). (5-12) is of the same form as (5-9), for if we replace arbitrary 
functions 9’ and o’ by —o and p respectively, (5-12) becomes (5-9). 


§ 6. Comparison of both solutions 


At first sight the exact solutions (GRS) and (NSS) explained in last section bear 
some resemblance to each other. If we compare them in detail, however, they differ 
remarkably in some points. We shall give an account of these circumstances in the 
following : 

(3) Sek, 
coincide with the corresponding one in the ordinary Maxwell theory. However, since 


and J’. F,;s in both solutions are of the same form (5-1) and 


both Riemannian space-times are not the same in general as will be seen in (ii) below, 
F,?’s, F’’s ahd F%’s are not necessarily of the same forms. 

Next we shall consider the charge-current vector J’. Js are zero in both solutions. 
But the circumstances concerning this vanishing of J’ are not the same. In solving 
(GRE) we have put J’=0 from the outset. This is seen in the second equation of 
(GRE). On the contrary, in solving (NSE) we obtained (NSS) without assuming 
J’=0, and obtained it as a result. (See (5-11).) 

Anyhow F,/s and J”s (=0) are of the same forms in both solutions. 

(ii) Relation between F,, and the Riemannian space-time. First we shall consider 
(GRS). In this case g,; and Fi; are connected by (5:2). When F,, is given, (5-2) 
is a differential equation for 7,,, and conversely if Y,; is given, it becomes an algebraic 
equation concerning Fy. Thus (5-2) indicates that the space-time and the electromagnetic 
field cannot be independent of each other. In other words, (5-2) shows that the 


electromagnetic field participates in the structure of the space-time through its energy 
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Ey=a(~, o)/m. When and only when p=o7=0, i.e. when the electromagnetic field is 
absent, //,; becomes a solution of the purely gravitational equation for empty region. 

If we consider the same problem in the case of (NSS) the situation ts quite different 
from the above. In this case, h,;, the quantity which determines the Riemannian space- 
time, is a solution of (5-8) which has no bearing on F;,, and » and o are quite 
arbitrary. There is no relation between h,; and f;;. This circumstance bears some resem- 
blance to the case of the ordinary Maxwell theory which admits various plane waves in 
the same Minkowskian space-time. After all, gravitational waves and electromagnetic 
waves which are propagating in the same direction and are independent of each other 


can coexist in NS. 

As was stated above, g,; in (GRS) is a solution of (5-2) and 4,, of (NSS) is a 
solution of (5-8). Since the quadratic form @(, 7) is positive definite, both Rieman- 
nian space-time can never be the same so far as F;;#0. Moreover, both Riemannian 
space-times can never be the same no matter what /’s and o’s we may take in both 


solutions. 
(iii) Superposition of the electromagnetic waves. A remarkable feature of the ordinary 


Maxwell theory is that the field equations are linear and that it is possible to superpose 
any electromagnetic fields. On the contrary, both field equations (GRE) and (NSE) 
which we are dealing with are very complicated non-linear equations for F,, and the 
fundamental tensor of the Riemannian space-times, and therefore we cannot expect any 
possibility of superposing electromagnetic fields in general. 

In spite of these circumstances, however, we can show that the superposition of the 
electromagnetic fields is possible in some sense concerning the solution (NSS). To see 
this we take two solutions (h,;; 0;, 7) and (h,;; ~:, 2) belonging to (NSS) and 
having the same 4,, in common. If we put (,;=¢,/;+0f and o,=c,o,+¢ , where 
c, and c are any constants, it seems most natural to interpret the words ‘the super- 
position of the electromagnetic fields is possible’ as that (h,;; 3, 03) satisfies (NSE) 
again. It is evident that this is the case with our (NSS). Accordingly, we can say 
in this sense that the superposition is possible. 

In the case of (GRS), on the contrary, the situation is quite different. Let two 
solutions having the same ,; in common be (His Ory 1) and (9453 2, 2). “Then 
these solutions are connected by @(p,, 7,)=a(f, o»). If vs and o, are any linear 
combinations of (,, (2 and o,, 7, as given above, we have a(i,, o,) a (¢;, %,) in 
general. Hence (/;;; (3, 73) cannot be a solution of (GRE) and the superposition in 
the sense above stated is impossible. Moreover, it is difficult in GR to make the super- 
position possible no matter how we may modify its definition. 

Summarizing the above, we can conclude that GR differs from NS considerably 
from the standpoint of the superposition of electromagnetic fields. 


$7. Concluding remarks 


We solved the field equations in GR and NS under the same conditions and by 


taking one of the simplest systems containing electromagnetic fields. Then we compared 
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the exact solutions thus obtained. Some different results as well as some common points 
were made clear. 

The most remarkable and important result seems to be the one stated in (ii) of 
the last section. That is, the structure of the Riemannian space-time is closely connected 
with the electromagnetic field in GR, while in NS the structure is determined quite 
independently of the electromagnetic field. The former is quite natural. On the contrary, 
the latter seems to be somewhat strange. The curvature of the space-time is independent 
of the electromagnetic field. What has caused the curvedness of the space-time ? 

Our present research is only a simple consideration concerning a system of a very 
restricted type. Therefore it is difficult to come to some definite conclusion concerning 
the problem that which of the two theories is more appropriate to describe the system 
containing electromagnetic fields. The writer hopes, however, that the present research 
will turn out to be of some use in solving this problem. 

The writer wishes to express his gratitude to Professor Y. Mimura and other 


members of the Institute for their helpful discussions. 


Appendix 


It is not easy to solve (NSE) under only two assumptions 1 and 2. Therefore 
we put an additional assumption which we shall explain in the following. 


It is known that F,, defined by (3-1) satisfies as a consequence of (NSE), 
FigntE pat Fu, j=9; (A-1) 


which is of the same form as the second half of the Maxwell equation given in (GRE). 
From Assumption 1, (A-1) and (4-4), we have 


F,=F,,=0, Fas —F,, 24; Fy= —Fy+b, (A-2) 


where a and 6 are arbitrary constants. If a or 6 is not zero, different from the case 
of the ordinary Maxwell equations or (GRE), the calculation to solve (NSE) becomes 
very complicated and it also becomes difficult to forecast the results. Because of this 
fact and the analogy with (GRS), we shall make the following assumption and_ shall 
deal with the simplest case : 


Assumption 3: 
a= b= 0k (A-3) 


As a consequence the electromagnetic field F,, becomes of the same form as (5-1). 


If we use this assumption we can prove (5:7), i.e. 


fo=fa=fatha=fotfu=0- (A-4) 


6),9) 


Proof. It is known that the following relation holds for non-symmetric metric 4;; 
Fy=afyt Besuf "/2; (A-5) 


where f* is the tensor obtained by raising the indices of f,; by means of h which is 
defined by h,,h*=0/, 
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a=—p/V —9, p=—b/y —-9, (A-6) 


(g=det. 9,,, h=det. hy; pHe fufu/8)> 


and p” is equal to det. f;;. 
We have from (A-2) and (A-3) 


Ofio+ Bf" = afytPf?=0, ey 
a (fs tha) +8 f 8—f) =4 (forth) —BF8—f) =0- 
Since the non-vanishing components of h'/(=h”) are 
b= —B/m, b?=D/m, P=—A/m, —h=h*=1/C, (A-8) 
we have 
baka Be PI Rie A Oe 
f= (Bhs—Dfis) /mC, f'=— (Bfha—Dfaw) /mC, (A-9) 
fP=—(Dfa—Afs)/mC, f= (Dfa—Afas) /mC. 
Then (A-7) becomes 
fa—B fin/C2=Bfn/m+afi.=0, 
(a—BD/mC) (fist fia) + (84/mC) (fos +fos) = 9, (A-10) 
(BB/mC) ( fistfia) — (@+8D/mC) (fos+fos) =0- 


Since cannot be zero, m>0O, and we are dealing with real quantities only, we can 


easily obtain (A-4) from these homogeneous equations for the components of f;,. 


Q. E. D. 
From (A-4) we can deduce 


p=0, and accordingly det. f,,=0, b=g, (A-11) 


which shows that our g,,; belongs to the third class in the sense of Hlavaty.” One of 
the great obstacles which appear when we try to solve (NSE) is the fact that gy and 
accordingly g and 4 become complicated. By the result above obtained, however, this 
obstacle is partly surmounted. 

The remaining calculations to solve (NSE) is given in the reference.” 

Lastly we add a remark. If we start with the identification (3-3) in place of 
(3-1), the relation (A-1) does not necessarily hold. But if we assume (A-1), (A-4) 
is self-evident from Assumption 3. 
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Along the same line on which we have previously discussed the cubic cerium metal, the magnetic 
properties of the succeeding rare earth metals, praseodymium and neodymium, are investigated. Because 
these metals have the structure of h.c. p., their crystalline field potentials cannot be uniquely determined 
so as to fit the experimental data. The analysis is, however, qualitatively satisfactory. 


$1. Introduction 


The magnetic properties of the rare earth metals have been discussed on the bases 
of Zener’s mechanism for indirect exchange interaction” and of the crystalline field 
potential” for 4f electrons. The characteristic features of these metals are thought to be 
due to the simultaneous effects of them and have actually been understood to some extent 
from this point of view. 

In the present paper we investigate the effect of the crystalline field on the magnetic 
properties of hexagonal rare earth metals. For this purpose we choose the praseodymium 
and the neodymium metals, because the effect of the crystalline field predominates the 
exchange field in these metals. The understanding of the crystalline field in these metals 
will be useful for further researches on complicated phenomena like ferromagnetic-anti- 
ferromagnetic transitions observed in other rare earth metals. 

Pr, unlike Ce and Nd, has no antiferromagnetic phase.” Pr shows a Shottky 
type anomaly in the atomic heat,”” whose maximum lies at a comparatively lower tempera- 
ture than Ce and Nd. The magnetic susceptibility of Pr at low temperatures seems to 
tend to a finite value at 0°K.” Nd has two anomalous peaks at 7.5°K and at 19°K” 
besides the same kind of anomaly in specific heat as Pr. Both peaks may be regarded 
-as the cooperative phase transitions and the former of them can be considered to correspond 
to the Néel point.” The latter peak is suggested to be due to the interactions among 
the multipoles of the 4f electrons.”» The anisotropic magnetic susceptibility of the Nd 
single crystal has recently been reported.” The electrical resistivity curves of Pr and Nd 
are concave downward over a wide temperature range and Nd shows another anomaly in 
the lower temperature range.” 

The phenomena mentioned above in the praseodymium and the neodymium metals 
can be explained semi-quantitatively by the view of the Stark splitting of the 4f electrons 
-due to the crystalline field. We shall show in § 2 the calculation of the crystalline field 
potential, the results of which will be applied to Pr and Nd in § 3 and § 4, respectively. 
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In the present paper, we confine ourselves to analyzing the specific heat and the magnetic 
susceptibilities above the cooperative transition points. In §5 discussion of our results 


will be made. 


§2. Calculation of the crystalline field 


While both metals, Pr and Nd, have the h.c. p. structure in the temperature range 
to which our interest will be confined, they have cubic closest packing structure at higher 
temperatures.” Calculation of the crystalline field is performed on the model similar to 
that for the cerium metal which has already been reported by us’: the valence electrons 
are treated as free electrons and smeared out uniformly over the whole of the metal, so 
that they have no contribution to the crystalline field discussed below. Therefore, the 
crystalline field can be obtained from the hexagonal closest packing of the ionized atoms. 

The crystalline field potential for 4f electrons in the h.c.p. structure is expressed 


as 
Vice ig ss re bs Vo = ya 6 ) 
where 


; lp» 
Vi?"= >} mx Pr” (cos4,) cosmo, , 
x 


—ve* 0 
ae Pe Aid P;’ (cos 4.) , (2) 


n 


lm 2 ese a is 
(im) tar RE 


P2" (cos4,,) cosm¢,,. 


In eq. (2), 7, 4, and g are the spherical coordinates of the /-th electron in the atom 
whose nucleus is taken as the origin, and R,,, 4, and ©, the spherical coordinates of the 
nucleus of the n-th atom, and v the valency of the ionized atom. On account of the 
symmetry of the crystal structure the sine term does not appear in eq. (2). 

In the lattice sum for />4, the direct sum method is applied over the range 
R, 4a, where a denotes the lattice constant along the a-axis. For /=2, the convergency 
of the series of the direct sum is rather slow. We therefore adopted the Fourier method,” 
whose convergency for /=2 is very rapid, and carried out for 4,<4/a in the sacieiaual 
space. The calculation is performed for c/a=1.579 and 1.613 and the results are shown 
in Table 1 and Fig. 1. To verify we calculated c, by the both methods for c/a=1.613. 
In Table 1 cell, n, indicates the contribution from (n—1)a<R,,< na for direct summa- 
tion and from n/a<h,< (n+1)/a for the Fourier method, except for n=1 where 
h,<2/a. As is seen in Table 1, the larger / is, the more rapid the convergency is. 

Cy needs some comment. As can be seen in Fig. 1, the value of cy, is approximately 
zero when c/a=y 8/3=1.633. This is due to the symmetry of the lattice. When 
we carry out the lattice sum for cy, the contributions cancel one another noticeably 
within the Ast cell and the total contribution from the Ist cell is reduced by a factor 
of 10°~10°, which is just the order of magnitude to be cancelled by the contributions 
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Table 1. The calculated values of Cim- (—ve?/a’*!) should be multiplied to each value in the Table. 
Ceg6= Cee X (2 xX 10395/12 !) ‘ 


cla Method cell C20 C40 C40 | cae 
1 +0.8355 — 4.688 +5.704 
Se 2 haa (202173 +0.196 —0.204 
3 | 0.0142 +0.011 —0.025 
sega 4 | +0,0195 +0,003 —0,001 
1.579 Total +0.5845 —4.478 +5.474 
1 _ +0.04799 oe 
Fourier 2 | +-0.07960 | 
3 |  +0.00155 
Total | +0.12914 
1 +1.0508 | —4461 | 45.753 
pits 2 —0.2260 EO. 74e | SY O497 
3 —0.0159 +0.010 |  —0,009 
a 4 —0.0234 FEO? —0,001 
yale Total +0.7855 —4.274 | 5.551 
7 > L Fae 
1 | —0.00598 0.5689 | | 
ae 2 | +40.05754 +0.2159 | 
3 +0,00152 —0.0062 | 
Total | +0.05308 +0.7781 | 


from other cells. (cf. Table 1). The 
convergency of the lattice sum is never- 
theless very rapid because the cancellations 
in the 2nd and 3rd 


in the lst cell 


of the lattice sum 
cells are less than those 
and they show the characteristic convergency 
of the Fourier method. 
of the energy levels of 4f-orbitals, the 


In the calculation 


contribution from the term V’,’ is equally 
important as the other terms, and ¢, is 
so small in magnitude and so sensitive to 
the value of c/a that it cannot have any 
rigorous meaning in view of our model 
used for the calculation of the crystalline 
field. In what follows, therefore, we 
modify D,, somewhat arbitrarily. (cf. eq. 
(3) and (4)). 

Eq. (1) can be conveniently expressed 
in terms of equivalent operaters as far as 
our interest is confined to the manifold 


of J2=const.,"” i. e., 


— 60 


€40 


0.9 


0.8 ++5.0 


0.7 


0.6 


1.560 1.580 1,600 1.620 1.640" 


c/a 
Fig. 1. The calculated values of c;,,. (—ve®/a’*") 
should be multiplied to each value of the ordinate. 
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V = Dy Joo+ DaJao+ DooJoo + Pov Soa. (3) 
where 
Ina 3le—J JED), 
Ja=35J2—30J (J +1) J2+25J2—9J J+) +3h ITD, 
Jog 231 Je? — 315J (J+) Je + 735Je' + 105)? (J +1)*Je 
—525J (J+1)J2+294J2—5]?(J+1)*+40f* J+ 1)* 
—60J(J+1), (4) 
Ja=3 Ui+iJ,)'+ e—td 
Dy, =? Cy Cm < (Coctheient of Fy) 


In eq. (4), @, for example, corresponds to the @ in Stevens’ notations,” and r is also 


due to his notation. 


§ 3. Praseodymium 


The most stable ionic state of Pr is Pr** whose ground state is “H,. The next 


excited state is *H, which is estimated theoretically to be higher by 2000 cm™’ 
) 


than the 
ground state.” Therefore, we discard the possibility of Pr°* being in the “H, state. 
The splitting of the energy level, “Hy, by the crystalline field (3) can be obtained in 
the same way as in the case of Ce”: In the present case, however, the difficulty arises 
from the fact that the relative magnitude of 1’, where / takes 2, 4 and 6, cannot be 
determined uniquely. There is no available calculation for the 4f wave functions in the 
rare earth atoms. We therefore employ for the radial part of the 4f wave function the 
Slater function containing the effective nuclear charge Z* as a parameter which will be 
determined below. In Fig. 2, the energy levels are shown as a function of Z* for the 
observed values of a4(=3.662A) and c/a(=1.613)." The | )’s in the figure indicate 
the following eigenfunctions. 


A1)=Na(|44)+Cai|#2)>) (doublet), 
A2)=N {| +4) +- C42| #2) (doublet) , 
B1)=By:(|+3>+Cm|—3)) (singlet), 


5 
B2) = Nyo(|+3)>+Cy2|—3)) (singlet) , ©) 
Ch=) ar) (doublet), 
D)=/0) (singlet) . 


In eq. (5), the state |+4), for example, represents the normalized wave function of 
Jz=+4. Ny is the normalization constant and is equal to (14+4-|C,,|*)~"?, where s 
stands for A or B and i for 1 or 2. As can be seen from eq. (4), only the states 
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that are different from each other by 6 in the 
value of J, can combine together, and the upper 


(or lower) signs of the former double signs in the 


right-hand sides of eq. (5) correspond to the upper 
(or lower) signs of the latter ones. The possible 
ground state is either |D) or |B1) depending on 
Z*. In either case, that the ground state of Pr’* 100 


is a singlet explains the fact that the praseodymium 
metal has no antiferromagnetic phase. Furthermore, 
the levels crowd relatively near the ground state 
and this fact is also consistent with the fact that 
the Shottky type anomaly of the atomic heat 
appears in a comparatively lower temperature range 
than in the cases of the cerium and the neodymium 


metals. 


E/ (10—5 ve2/a) 


With these theoretical backgrounds, our calcu- 
Jations are carried out for the atomic heat and the 
magnetic susceptibility. The fair agreements with 
the observations are obtained for Z*=10.5 with 
D,, multiplied by 2. The multiplication factor is 
necessary, for without such a modification for D,,, 
the results are not good regardless of the choice of 
Z*. In Table 2, we show the energy eigenvalues 


corresponding to the eigenfunctions in eq. (5) 


— 150 


together with coefficients C,, in the same equation. 
Fig. 2. The energy levels of the 4f 


electrons of Pr3+ (83H,) in the crystal- 
zero. line field as a function of Z*. 


The eigenvalue of the ground state is taken as 


Table 2. The energy eigenvalues and the coefficients of the eigenfunctions of the 4f 
electrons in the Pr metal. 


lai) | |42y | |Buy | 1B | icy | ID 
& : a 4 3 | | | cs 
(Energy eigenvalue) /(10~° ve*/a) 84.3 400 | 604 | 132 | 364 | 0 
Cai V eaen. 4 — 0.0865 | 1 | ge | 


(teh eee ein, Si eee a a A Es eel SS le 
By using the energy eigenvalues, E,, given in Table 2 and the multiplicity, ,, shown 
in eq. (5), we can calculate the atomic heat due to the localized 4f-electrons from the 
following equation : 
AC= Nk { (1/Z,) (3 nj E ge Pe 
— (1/Z,)? (S.mE se **)*}, (6) 


where Z,= 5,1; eo XE, 
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In above equations E,’s express the numerical values given in Table 2 (i.e. the value 
measured in the unit of (10~*ve*/a)) and x=10~'( ve'/a)/kT. The magnetic suscep- 


. ° - 2) 
tibilities are calculated in the same way as in the case of Ce”? and expressed as 


sie ee! : 
1, = Nts) (Fit —G,), 
0 


kT Zz (7) 
N(gpp)? 1 2 . 
ti CUE (Fi heat }: 
EY ya x / 
where 
F', =0.000444 (17240 e+. 70540 e-4™) 4 2e-¥e 
8432 _—-400x -604r 139, 
Gp 1.435 See eas ae ae (8) 
316 718 
0, 
—12r 400 —Siae 1s 
G,=0.0149 {331 SAC ETS EY, eM S 
268 47.9 
6042 44004 —60.4e —84.Sx -—M 4c —400r 
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Fig. 3. The part of the atomic heat due to TK as a 
the 4f electrons of Pr metal. The calculated Fig. 4. The magnetic susceptibilities calculated 
curve is shown together with the experimental for the Pr metal are shown together with the 


Ones. experimental curve. 
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In these equations, N is the number of atoms, J the Lande y-factor equal to 4/5, My 
the Bohr magneton, and Z, the partition’ function. 7% means the susceptibility when 
the external magnetic field is parallel to the c-axis and Y , the one when the magnetic 
field is perpendicular to the same axis. F , and 

F’, are the contibutions from the matrix elements Eo 

between degenerate states, and G, and G, are those 
from the matrix elements between non-degenerate 
states. The exchange effect is omitted here. 100 
The results are shown in Figs. 3 and 4 together 

with the experimental data. In the actual 
calculation v is taken as 3. In Fig. 3, the 
experimental curve indicated by (Pr-La) means 50 


the excess part of Pr atomic heat over La’s.” 


The other experimental curve is the one obtained 


with the Debye temperature 6,=138°K." In Fy ’ 
Fig, 4 ¥=2%,/3+2%,/3, which should be »% 
compared with the experimental data. 7%, and & 
7%, are very sensitive to the character of the zl 
ground state but there is no observation to be —50 


referred to. 


$4. Neodymium 
— 100 


We treat the neodymium metal in the 
similar way as for the praseodymium metal. 
The ground state is “I,. and the next excited 
Srey ‘Tnpes estates Di swtina Gi ah ene Fig. 5. The energy levels of the 4f electrons 
fore, only J=9/2 manifold is considered. The Re Nar Cine, i the cepstalling eld 
observed values of c/a are 1.613" and 1.609," 
but the former is taken here with the value of a=3.650A. For Nd, there remains 


— 150 


Kramers’ doublet on account of the odd number of the 4f electrons, and all levels are 
doublets as is shown in Fig. 5. The notations in the figure are similar to the case of 


Pr. Those eigenstates are represented by the following expressions : 
AIY=Nui(|+9/2)+Cu|+3/2)) (doublet), 
AL SSN 40 ( +9/2)+ Cy + 3/2)) (doublet) , 


B1)=Npi(|47/2)+Ce|F5/2)) (doublet), (9) 
B2)= Np (| £7/2) + Co) 3/2)) (doublet) , 
Gaal 12) (doublet) , 


where similar notations have been used as in eq. (5). 
The calculations have been performed for Z*=11.0 with D,, multiplied by a factor 
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2.5. The coefficients in eq. (9) and the energy eigenvalues are listed in Table 3. The 


zero of energy is set on the ground state. 


Table 3. The energy eigenvalues and the coefficients of the eigenfunctions of the 4f 


electrons in the Nd metal. 


(41) es aap Bi) B2) C) 
} = 
(Energy eigenvalue) /(10~* ve"/a) | 16.7 197.9 0 137.8 76.7 
Cx Lin SS 0.127 —3.70 0.270 


In this case, the specific heat of the Shottky type can also be calculated by eq. (6) 
with E,’s in Table 3 and n,’s in eq. (9). The magnetic susceptibilities are expressed 
by eq. (7) with g equal to 8/11 and with F’s and G’s given by the following : 


B= 8.77 19 et 394 ee eT en ot 
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Fig. 6. The part of the atomic heat due to the 
4f electrons of Nd metal. The calculated curve 
vis shown together with the experimental ones. 


Fig. 7. The magnetic susceptibilities calculated 
for the Nd metal are shown together with the 
experimental curves. 
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The results are shown in Fig. 6 and 7, where v=3 is also taken and the notations are 
similar to those in the case of Pr. Again the exchange interaction is not taken into 
account. 

As is seen from Fig. 5, there is one doublet near the ground state and this situa- 
tion does not change even if some modification of D,, is made. This explains the 
experimental fact that the Nd metal shows JS=Rln4 below the second cooperative 
transition point which lies at a comparatively low temperature.” In Table 2, we show 
the excess entropy corresponding to the anomalous atomic heat, JC,. In the last column 
' there are also shown the theoretical values to be compared with. As has been discussed 


in §1, the cooperative transition points are at 7.5°K and 19°K. 


Table 4. The entropies related to the excess atomic heats of the Nd metal. 


T’K | 4S (Nd La) |  48(@p=141°K) | 
10 1.18 1.2 | Rin2=1.38 


75) 2.02 2.2 Rin4=2.75 


The magnetic susceptibity has been measured for both single-” and poly-crystals.?™” 


The observed anisotropy of % is smaller than the result obtained here, but it is qualita- 
tively understood that 7; is larger than %, and that the sudden change of the slope of 
the experimental 1/%—T curve is probably related to the effect of the crystalline field. 


§ 5. Discussion 


So far we have used the three parameters, v, Z* and D,,. wv determines only the 
absolute scale of energy as can be seen from Figs. 2 and 5. wv=3 seems to be some- 
what large compared with the case of the cerium metal where v=2.5 was consistent.” 
But it would not be unreasonable. The values of Z*, 10.5 and 11.0, should be com- 
pared with 12.65 and 13.3 obtained from Slater’s rule for the Pr and the Nd atoms, 
respectively. Some discrepancy between the Slater function and the wave function 
calculated numerically by the Hartree method can also be found for the 4f wave func- 
tions of Au and Tl atom.’” It is difficult to make clear the meaning of the modification 
for the D,, on the basis of the present model of the crystalline field. We can probably 
say that a small departure from our model for crystalline field may be responsible for 
this. The results obtained in this paper should be understood to be rather qualitative or 
semi-quantitative. It cannot be uniquely determined whether the ground state of Pr is 
|D) or |B1). <A better agreement with the observations can be obtained if we choose 
|D) as the ground state and the realization with proper z* and v relates to the modifi- 
cation as D,,. To explain the susceptibility and the specific heat by taking |B1) as 
the ground state, a fairly large exchange interaction is necessary. A large exchange effect 


is not plausible, so that this situation is abandoned in this paper. 
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The discrepancies between the present analysis and the experimental data can be 
related to the following facts : 

(1) The modified hexagonal structure has been reported.”"? 

(2) The temperature dependency of c/a has not yet been observed, which sensi- 
tively affects Dy). 

(3) The conduction electrons may deviate from the free electrons and the polari- 
zation of the core electrons, especially those in 5s and 5p orbitals, may arise due to the 
4f electrons. 

(4) There is a lack of our knowledge about the 4f wave function in the metal. 
It is possible that there is a marked deviation from the atomic or Slater’s function. 

For example, the larger anisotropy of the susceptibility for the Nd metal which is 
obtained in the present calculation would be intuitively understood in view of (1) and 
(3) mentioned above, because these factors have a tendency to mix the levels further 
than those indicated in Figs. 2 and 5 and diminish the anisotropy. Detailed considera- 


tions on these points will be made in the future. 
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A method suitable for analysing the anomalous f-spectrum of RaE is developed. It is shown 
that, if the Gamow-Teller component of the f-decay interaction is of the tensor type, it is very difficult 
to explain the shape of the #-spectrum of RaE unless all of the following conditions are satisfied 
simultaneously: (a) the Fermi component is of the scalar type, (b) Cs/Cs’=Cr/Cr’, and (c) 
Cs/Cr=real. Alternatively, if the Gamow-Teller component is of the axial vector type, the condi- 
tions are: (a) the Fermi component is of the vector type, (b) Cy-/C’p=Cy4/Cy4’, and (©) Cy/C4= 
real. 


$1. Intreduction and summary 


Since the experiments by Wu et al’). showed the violence of the space reversal in- 
variance all our previous knowledge about the type of the /?-decay interaction has been 
undergoing a thorough re-investigation. Although many recent data on angular distribu- 
tion of /-rays from polarized nuclei, longitudinal polarization of f-rays, and angular 
correlations between -rays and circularly polarized 7-rays, etc. (except for the data on 
electron-neutrino angular correlations), gave us a fairly consistent picture of the inter- 
action type, they are still not sufficient to determine it definitively. 

Though the peculiar (-spectrum of RaE is seemingly difficult to explain, it is now 
well-known that it can be reasonably explained as an accidental case in which consider- 
able cancellations occur among several partial matrix elements”. The aim of this paper 
is to re-investigate the cancellation spectrum by taking into account the fact of parity 
nonconsetvation in [-decay and to obtain additional information about the type of the 
P-decay interaction. A similar analysis was made by Lewis”. 

Through out this paper the following assumptions are made: 1) The Fermi theory” 
in a broad sense is valid, in which the four spinor fields interact at a point without 


‘ : 5 ‘be ; 4 
derivative coupling”. Parity nonconservation is taken into account.” 


* A preliminary report of this paper was published in Phys. Rev. 108 (1957), 1104. 
* On leave of absence from Department of Physics, University of Tokyo. 
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2) The time reversal invariance of the Wigner type” is valid for strong interactions, so that 
, ne é : 
the relative phases of nuclear matrix elements can be fixed, while no such assumption 


. . . 9 
is made for the weak interactions’. 
If the parity nonconservation is taken into account, the interaction Hamiltonian for 


B-decay can be written as a sum of two parts, one with coupling constants C,’s and the 
other with Cj’s in the notation of Lee and Yang”, where i, j=S, V, T, A and P. For 
B-spectra there exists no interference term between these two parts. Therefore, the only 
quantities which can give us important knowledge about the validity of the two-component 


neutrino theory” and the time reversal invariance from study of f-spectra are 
a= Re (CC, + C/*C/)/{ (ICP+ ICLP) IGP + ICY. (1) 
Here Ste means the real part. a, is always equal to or smaller than unity. If we put 
i=S(T) and j=V(A) in (1), the @,; appears as the coefficient of the Fierz term”. 
The empirical fact that the Kurie plots of the allowed spectra are nearly straight can be 


explained by the Fierz conditions, 
Asp=0, Ap4=0. (2) 


On the other hand we can obtain other conditions from cancellation spectra of 
forbidden transitions such as the decays of RaE and C,'", because in these cases the 
interference terms involving as, and @;4 play the most important parts. By the analysis 
of RaE spectrum we are led to the conclusion that the peculiar shape can be explained 
only by the optimum effect of the interference terms and even a small deviation from 
the optimum situation makes the explanation remarkably difficult. The optimum situation 
is realized by the following conditions : 


Caf Cs = Ge Go, Sal Ge = GGL 
C;/C,=teal, C,/C,=real. (3) 


§ 2. Theoretical expression for the correction factor 


Deviations of -spectra from the allowed shape are expressed by correction factors 
introduced by Konopinski and Uhlenbeck™. Though their formulation is useful, several 
corrections must be taken into account in order to make it more accurate. Among them 
the effects due to the finiteness of nuclear size are considered to be the most important 
for RaE decay. They are composed of the finite nuclear size effect in the narrow sense’, 
the finite de Broglie wave length effect’, and the effect of the 8-decay from inside of 
the nucleus. These effects were originally proposed independently of each other. 
Takebe’ proposed a general formulation including the above corrections synthetically, but 
the number of unknown parameters seems too large for practical purposes. 

. ey two of the authors" proposed a formulation convenient for practical uses 
in saeenie B-spectra taking into account the above mentioned effects and the screening 
oe # arosnetically. It was shown that the effects due to finite nuclear size can be treated 
without detailed knowledge about the nuclear charge distribution if it is spherically sym- 
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metric, so far as we regard the nuclear matrix elements as unknown parameters. Namely 
the main part of this correction can be included legitimately by only replacing the 
ordinary nuclear radius in the Konopinski-Uhlenbeck’” formulas for aZ<1 (q@: fine 


structure constant, Z: atomic number) by the ‘“ 


effective nuclear radii”, which are ex- 
pected to be about the same as or somewhat larger than the ordinary one. Moreover, the 
unimportance of the screening effect on the RaE spectrum (except for the very low 
energy region of W7<1.1) is suggested even if the cancellation is remarkable. 

_ According to this theory the general correction factor for the RaE spectrum (spin. 


change 170°), Cispra(W) is given as follows : 


Cisrra (W) = Cierra (W) + Cerra W), (4) 
with 
aZ 
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and C(%»4(W) the expression (5) with all the coupling constants replaced by primed 
ones, C,/. In (5), ¢.c. means complex conjugate, p and q represent the electron and 
neutrino momenta respectively, and W the electron energy including its rest mass. We 


* These approximations correspond to the omission of aZHp in comparison with unity. F,(o) and 
In(r) are included only in the L, term!)') which has no importance in RaE decay. _Although the approxi- 
mation Ip(r)=I,(r) reduces the number of nuclear matrix elements by one, it is reasonably permitted because 
-of the absence of interference between the L, term and the main term. The approximation of, J#=y*=y* 
=1 and y~=0 concerns with the main term, but does not reduce the number of nuclear matrix elements. 
Thus the effect of this approximation is considerably smaller than the effect of the new nuclear matrix elements 
which will be discussed in 24 as neglected terms in the effective radii theory. 
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use the natural units (6=c=1) and also take the electron mass as unity (m=1). The 
first forbidden transition including only scalar, tensor and pseudoscalar types is explicitly 
treated in ref. 17. Expression (5) is obtained by the general prescription in ref. 17 with 
additional approximations of FOR @O=1, £0) =L(r) {=n =)" = and v-=0 
in the notation of ref. 17.% We can improve (4) by removing these additional approxti- 
mations, but the degree of improvement is so small and insignificant for the following 
argument that we use (5) in this paper. 

_ Each term of (5) clearly corresponds to that of the familiar Konopinski-Uhlenbeck 
formula for low Z. The energy independent function I,(r) is defined by “eq. (13) or 
tef. 17, and its behaviors in nuclei are shown in Figs. 3° to 6 of ref. 17. The effective 
nuclear radii, 0,(X)’s are defined by eq. (21) of ref. 17 as the ratios of nuclear matrix 
elements. Different ,’s have the same order of magnitude but may take different values. 
Their magnitudes are comparable with that of the ordinary nuclear radius and are likely 
to be somewhat larger. Since the exact values of the nuclear matrix elements are not 
known, we need to treat the effective nuclear radii as flexible quantities to some extent. 

It is convenient in our analysis to expand eq. (4) into a finite power series in the 
electron energy W by using the relations, q=W,—W (W,: the maximum electron energy 
and p=W?—1. 


Cisrra W) =a_4/W+a+aW+aW?, (6) 


where a,'s are energy independent coefficients. 


The explicit forms of a, and a, are 
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9 5 9 ( 9 

+ Aej(@ pt esljexrhO tt (8) 
where we defined €; as €&= (ete py 

It is interesting that 4, is always positive, 4, > 0. 


If we assume the Fierz condition (2), (6) reduces to a simpler form. 


§ 3. Analyses of the experimental data 
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ig. i Experimental correction factor for the RaE #-spectrum obtained 
4 19 ‘ 
ae ssman and Langer A The error, — 1%, of each count is 
taken into account in addition to the uncertainty, +: 1%, of the 
maximum electron energy. 


$3.1 The method of W expansion 
The -spectrum of RaE has been precisely measured by many persons’. The: 
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experimental correction factor obtained by Plassman and Langer’) is shown in Pigs thy 
in which we determined the width of the-curve by taking account of the uncertainties, 
+1% in the maximum energy W, and. +1% in each count. In order to compare it 
with our effective radii theory (eqs. (6) to (8)), we expand the experimental curve 
into the finite power series in W as shown in eqs. (6). By tedious but straightforward 
‘calculations we can prove that the experimental curve in Fig. 1 leads to 


yy) a7 ee ote, (9) 
which is implied in (8). 


On the other hand, using eqs. (7) and (8) and the Fierz condition (2), we get 
the theoretical expression for a,/ay as: 
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From (10) we can easily see that the interference terms between S and T and between 
V and A (asp and ay, terms) play important parts to make —d,/dy large enough to 
fit the experimental data. In fact, if we put @~=0 and ay,=0 in (10), we obtain 


the inequality relation, 
—4,/a,< W,=3.29, CL) 


unless we take negative values for p, (fr) and p,(0 Xr). 

Here, let us consider about these effective nuclear radii’. An effective radius is 
defined as a ratio of two nuclear matrix elements whose operators have the same angular 
part but somewhat different radial dependences. We can divide one of the operators into 
two parts, one (c) has the same radial dependence as the other operator and another 
(d) has an entirely different radial dependence. Let us call them part c and part : 
respectively. The magnitude of the part d as an operator, not as a matrix element . is 
about 1/25 of the part c in our case’. If we neglect the part d we get an effective 
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radius of about 1.1 times the ordinary radius. An effective radius much different from 
this value takes place only when cancellation occurs in the radial integral of the part c 
and not in that of the part d. In order to get a negative effective radius, cancellation 
dows forvabout 1/25) isiiecessary.| (Phat is ‘to sayvin yourmeases the radial integral of 
every concerning nuclear matrix element, jfrl,(r), |@Z/(2/: (8r)) (Pr @).§oxrl, (r), 
or [aZ/(2p,(6Xr)) Jo xXrl, (r), is required to take a valu2 of about 1/25 of the 
possible maximum value. On the other hand the nuclear matrix elements j/rJ,(r) and 
{oXrlI,(r) are essential to make the absolute value of a, large enough, as can be seen 
from (5) ot (7). In order to produce a strongly energy dzp2ndent correction factor as 
that of RaE decay with (uncorrected) log ft=8.0, either j3r I,(r) or foXrI, (7) must 
be fairly large. Magnitude of a nuclear matrix element |X can be expressed by a correct- 
ed ft value, f(X) t=27°In2/(&?||X|"), where €, is the relevant coupling constant. f(X)t 
is inversely proportional to ||X|*. The magnitudes of jr [,(r) and joxr L(@) 
necessaty to explain the absolute energy dependence of the RaE spectrum correspond to 
log f(Srl,(r)) 59.4 and log f(eoxXriE@)is9.4* 

If log f (Sr I,(r))é or log f(oXrlI,(r))t=9.6, we get a correction factor (log ft= 
8.0) whose lowest energy patt is at most only 150% larger than the highest energy 
part. Now all the known log ft values for the first forbidden transitions are not smaller 
than 5.222, The smallest values occur for JJ=0 transitions. If we take account of 
the spin statistical factor, (2J,+1)/(2J;+1) (=1/3 for RaE decay) or more con- 
servatively | (2J;+ 1) /(2J,+1) ]'"**, tae possible minimun log ft valu: for J=1~ 0° 
transitions is inferred to be about 5.5. It is very likely that this minimum value occurs 
only when there is no appreciable destructive interference in any respect. Thus we arrive 
at the conclusion that for a J=1~—>0* transition the possible minimum values of log 
f(Pr[,(r))t and log f(oXrlI,(r)) are equal to 5.5+2.3=7.8, where the 2.3=log 
[aZ/ (2) P came from the Coulomb factor. Comparing this 7.8 with 9.4 which we 
obtained as the possible maximum value for RaE decay, we see that the radial integrals 
of [fr I,(r) and joXrl,(7) take the value of at least 10°°""°"<1/6 times the 
possible maximum value. This result excludes the possibility of negative effective radius, 
for which the value of <1/25 is necessary instead of 1/6. Thus the variations of 
o,(8r) an ~,(0Xr) are confined to 


* : 
The rough evaluation of these upper bound can be made by using the approximate form of the 
correction factor, 


AG Ae x oa ceo 
Stas Wi xle+ 9 Mex, 


where & and ja|<1 are energy independent and x is €s \ er I(r) or &4 \oxr I). The corresponding 


f(x) values is 
( o if. -~ 4 » o 
f(X) =\FoW, Z| P+esWwe + 9 Ww |pira? dW. 


The uppper bound given in the text is taken somewhat larger than this rough estimate making allow 

for the incorrectness of the above formulas. y ow 
’ We cannot define the general statistical factor uniquely. The former example in the text is onl 
the ratio of the number of final substates of the decay in question to that of a hypothetical 4J=0 nS. 


starting from the same nucleus. The latter is su db iti 
) ; ggested by the fact that the reduced ility i 
(2J¢+1)/(2J;+1) times that of the reverse process. a ade da sean 
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(1—6/25) 1.19 >,(8r), p,(oXr) < 1.1o/(1—6/25), 


eee 
0.8¢ <1 (3r) > C4 (o me r) ded 0. (12a) 


For the later purpose (see §4) let us consider the case where we take 9.6 as the: 
possible maximum value of log f(rI,(r))t and log f(oXrI,(r))¢ for RaE decay. In 
this case the radial integrals of |r I,(r) and {oXrI,(r) take the values of at least 


1/8 times the possible maximum value. The variations of (i's are 
0:7 Pa (Pr) A, (OXt) <1.7/. (12b) 


Thus we were able to establish (11) on solid foundation. 
The relation (11) clearly contradicts with the experimental fact (9). Thus we 


can reasonably exclude the case of 
Wn 0 and ay,—0. (13) 
It is interesting to note that the above relation holds for the following special cases : 


a) Only one type of S, V, T or A. 

b) VT type. 

c) SA type. 

d) ST type with C;/C,;,=C,’/C,’=pure imaginary. 
e) VA type with Cy/C,=C,’/C,’=pure imaginary. 
f) ST type with C,/C,’=—C,'/Cp=real. 

g) VA type with C,/C//=—C,'/C,=real. 


In the cases d) and e) the two component neutrino theory” is valid, but the time reversal 
invariance” is not. In the cases f) and g) the two component theory is not valid. In 


the next two subsections we investigate the special cases ST and VA in detail. 


§3-2 ST type 


In this section we examine a possible special case, ST type, suggested by the data 


on electron-neutrino angular correlation of He’ *. In this case (5) can be written as. 


hy 20,5 a) yo LF AI ee Cy SG ieee aN 
See [Poxeh(n| x hee (8r) 3 (9 W  20,(8r) 
G+tP 24f"| yz, {2 (4 haa Vesey 
Gime ae 4 pie os are 
p’ EE aa a i Bl +2 |, (14) 
eola ig )\X+4 3 ree she ranean 
where 
(eam) as ees ya fel Prh (15) 


| jfoxrl, (r) 2 ( (80Xr),” E,{ Bo XrI, (7) 
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According to the usual method for analysing /-spectra, we choose the three electron 


energies, W,=1.72, W.=2.25 W,=2.75, arbitrarily. The requirement that the theore- 
tical correction factor (14) runs through the experimental region of C(W) on Fig. 1 
at W=W, and W,(or W=W, and W.,) restricts the allowed region of the parameters, 
’X and Y. In the actual analysis we normalize C(W) at W=1.72 according to Plass- 


‘man and Langer’. We obtain the ratios 
C(W,) /C(W.) =1.40 to 1.50, (16a) 
C(W,) /C(W,) =1.96 to 2.32, (16b) 
taking into account the experimental errors. The allowed region of the parameters, x 
and Y, lies between the four curves given by eqs. (16a) and (16b) with the limiting 
values in the right-hand sides. 


Ww 
20 Y 
jg 
U/ J 
ia 
4 £4 
of 7 
1.5 - Uj Y 
ly \ 
bye 
, 
i, 
Ly 
na le 4 XL 
a a 
We 
Y 
if 
03 + rp xy 
LY 
Vy 
Oe | 1S : | Ve ia J 
: 10 20 30 x 


Fig. 2. The shaded area is the allo i 

Ph wed region for the parameters, X and Y, i 

oe ore. Y and LW are the points used by Yamada’) and Lee-Whitin; pide bie ‘Fes 
a an the correction factors are shown on Fig. 3. aie ake 
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Fig. 3. Correction factors for ST with as7=1. 


In order to compare the effective nuclear radii theory with the previous analyses, 
we first treat the case, @sy—1. For simplicity, we adopt a fixed value, W,=3.29, for 
the maximum electron energy in the theoretical expression (14), because the uncertainty 
is already taken into account in Fig. 1. Furthermore we put aZ/| 20, (Sr) ]|=13.25, in 
which the effective nuclear radius /,(j3r) is assumed to be 25% larger than the ordinary 
nuclear radius, ¢=1.2A'” 107" cm. The region of the parameters thus obtained is 
shown as the narrow shaded area in Fig. 2. The previous parameters used by Yamada” 
and Lee-Whiting” are plotted as the crosses in the figure. Our allowed region may be 
shifted to include these crosses by choosing a somewhat different value for /,(fr). To 
visualize our results we plot in Fig. 3 the correction factor for special choices of the 
parameters, X=16.7 and Y=1, and X=24 and Y=1.5. It should be noted that the 
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_Y. Therefore the sign 


above results are unchanged for @s~= —1 if we replace Y by : 
know the relative signs 


of Asp cannot be determined from the present analysis unless we 


of the concerning nuclear matrix elements. 
In this case it is more difficult to explain 


The restrictions for |@s,| are sum- 


Now let us examine the case, |@sp|<1. 
the experimental shape than in the case, |@s7|=1. 
marized in Table I for a few assumed values of p,(8r)/p which cover the allowable 


variation discussed in the previous subsection. 


Table I. Restrictions for degrees of interferences between S and T, |@sr'; and between V and A, |@r Als 
which are obtained from analyses of RaE spectrum with the effective nuclear radii theory. Only 9;/p values 


satisfying (12a) can actually occur. 


type Siu VA 
Eq. used | (16a) | (16b) (16a) (16b) 
0;/0 \esr|> |@sr| > aya\> lar al > 
0.4 0.9991 0.9993 0.9995 0.9996 
0.6 0.9980 0.9985 0.9988 0.9990 
1.25 0.992 0.994 0.995 0.996 
2 0.982 0.987 0.989 0.991 
3 0.966 0.975 0.979 0.983 
Ce a ee eee 


For the excluded values of |a@s,|, at least one of the left-hand sides of (16a) and (16b) 
is always smaller than the right-hand side. The allowed deviations of |@s,| from unity 
are quite small irrespective of the values of the effective nuclear radius, suggesting the 
actual occurrence of the maximum interference, |@s7\=1. 


§3-3 VA type 


Another possible type besides ST is VA as suggested by the recent experiments on 
electron-neutrino angular correlation of A35° and on circular polarization of 7-rays follow- 
ing the electron capture of Eul52". The expression of the correction factor for VA 
can be obtained only by reversing the sign of the neutrino momentum q in the correction 
factor for ST in § 3—2, hence they are quite similar. 
and Y, as 


? 


Defining the parameters, X 


az 
20, (r) ” 


= jal, (r) 
ijrl, (r) 
— &4{O Xr, (7) 


i€) (rl, (r) (17) 


we get the following theoretical correction factor : 


Cip4(W) =E2|\r 1, (A) [2° | X24 ae te is (etl Se ame 
W) =e | [rnp | Xe+ | ( rg oe eae, 
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Fig. 4. The shaded area is a part of the allowed region for the parameters, X and Y, in the 
VA case with @y4=1. For c and d the correction factors are shown on Fig. 6. 


Analyses were made in the same way as the ST case. First, for a@y4,=1 allowed 
regions for the parameters, X and Y, are shown in Figs. (4) and (5). The only main 
difference in this case is that the allowed parameter regions exist for both signs of Y 
contrary to the ST case, two typical correction factors are shown in Fig. 6. The fit is 
as good as that for the ST type. For a@y,4=—1 the results are the same if we reverse 
the sign of Y. 


Then we calculate the case of |a@)4|<1. The allowed regions for |@y4| are sum- 
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Fig. 5. The shaded area is a part of the allowed region for the parameters, X 
and Y, in the VA case with ay4=1. 


marized in Table I for several assumed values of ¢,(0 Xr) / which cover the allowable 
variation of the effective radius. The smallness of the allowed deviations of |a@;, from 


unity suggests the actual occurrence of |@p4|=1. 


§ 34. SVTA type 


The case of SVTA mixing is qualitatively surveyed in this subsection. It can be 
inferred from the argument of § 3-1 that the RaE spectrum favors |@s7|= /@y4 =1. In 


such a case all the Fierz terms (terms having 1/W factor in (5)) vanish, and the cor- 
rection factor can be written as: 


Cispra (W) = Cisr (W) =P Cir (W) : (19) 


Therefore, in the case of SVTA mixing, independent cancellations must occur simul- 
taneously in both parts of ST and VA in order to explain the RaE spectrum. Although 


we cannot exclude this possibility definitely, the occurrence of such an event is expected 
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Fig. 6. Correction factors for VA with ay4=1. 


to be highly accidental. Hence, the RaE spectrum favors either ST or VA, but not 
SVTA mixing, from the view point of probabilism. 
It can be shown from analyses of Figs. 2,4 and 5 that if the relation, 


EsE4/EpEp=YY' > 0.16, 


is satisfied for the case of |@s7|=|@y4|=1 there certainly exist suitable parameters to 


explain the experimental correction factor. 


§4. Evaluation of the effects of neglected terms 


Although we obtained very stringent restrictions for |as,| and |a;4| in § 3, the 
arguments were based on the validity of the effective radii theory. This theory seems 
to be sufficiently accurate for usual analyses, but in order to deduce definitive conclusions 
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we must examine the terms neglected in that theory. 

Aside from the radiative correction, three approximations are made in the effective 
radii theory: 1) approximation to the inner wave functions of electron, 2) approxima- 
Fon to the normalization of the wave functions of electron, and 3) restriction to the 
number of partial waves to be taken into account. The effect of the approximation 2) 
does not seem to be significant, because no new nuclear matrix element is concerned and 
the effect is not magnified even when there is a large cancellation’”. On the other hand 
new nuclear matrix elements appear if we eliminate the approximations 1) and 3), and 
their effects are magnified when cancellations occur among the main terms. Therefore, 
we evaluate the effect of the approximations 1) and 3) in this section. 

If we denote the correction factor Cis7(W) or Ciy4(W) in the effective radii theory 
C,(W), the correction factor that includes the above mentioned corrections are approxi- 


mately given by the following expression : 
CW) ~C,W) £2|GW) |V CF.W) + |G) ', (20) 
where the G(W) is 


Ww? hog e 
Csr) =—" he, ([Bar | dr, ri (oh) 


ZW TL) Tae | K® 
ar 3W De 7 ( pa Sees), 


ad eee os =& 3ar°I, (r 
reer r(\2ar°I,()) 


eae aW eld 
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tee 

2 ates 1) _ His (jark, Y(r)) 

Si aL &((r(Sa-r)1 
+7 Ce any PIT CaO), (218) 
“for ST and 
Ww? niece 
GraW) =F he, (lar dr, nV) LH) 
W Jo 


2(2W?+1) 71 ¢ : 
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ER CLL 8) iy mags | i 
—_ IP UT em K& 
3W Due rin (7))) 
24 az 

3 20; (r (a i ry) 


&_((r (a-r)I,(r)), (21b) 


for VA. K,(r) and K,°(r) are defined in eqs, (3) of reff 17), .D,, in eq. (12) \or 
ref. 17) and 7,=)/1—a@’Z?. In (20) the second term in the right-hand side is the 
main correction term, and shows strong energy dependence. Its form, given by (20) 
and. (21a) or by (20) and (21b), is approximately the upper bound (in magnitude) 
of the exact and more complicated form. The last term in (20) is generally small, 
and has a negative effect to produce strong energy dependence as is required for RaE 
decay. Its form given by (20) and (21a), or by (20) and (21b), is the lower bound 
of the exact form. Thus in the approximate form (20) with (21a) or (21b) the 
energy dependence due to the correction terms is fully taken into account. 

In order to see how much effect the correction terms could have, we first consider 
the upper bounds of the new nuclear matrix elements. Generally it seems reasonable to 
set the upper bound of the nuclear matrix element \WU, (r), where 7 is independent of 
r, as Re (r) dr/ Wales (r)dr times the upper bound of \nU, (r). Using this relation 
and assuming the uniform charge distribution for the nucleus, we can relate the upper 
bounds of the nuclear matrix elements, except \r (8a-r) LG) and|r (a-r)I,(7), with 
the familiar first forbidden matrix elements. We do not know much about the new 
matrix elements, \r (Sa-r)I,(r) and \r(a-r) I(r). We assume that their upper bounds 
are equal to those of j BarI,(r) and \ orl, (r) respectively. According to the argument 
of § 3.1 we can set the upper bounds of Jal, G); \ eal, (ar \rl, (r), \Prl, (r), \oxrl, (r) 
and \3o Xr, (r) as those corresponding to the corrected log ft values (log f(a@I,(r))¢ 
etc., (see § 3.1) of 5-5,5-5,7-8, 7-8, 7-8 and 7-8 respectively. We take these upper 
bounds of the nuclear matrix elements and give them appropriate signs and put |\a@s7|= 
|\@y4|=1 so as to make |Gop(W)| and |Gy,(W)| maximum. In this procedure the 
relative signs of the first three terms in the right-hand sides of (21a) and (21b) were 
not adjusted because they are similar and not independent of each other. Using those 
IG(W)|’s and the experimental C(W) (including normalization) , we can calculate C,(W) 
from (20). This C,(W) has the possible least energy dependence. C,(W)’s thus 
obtained for ST and VA are shown in Fig. 7. To obtain them only the C(W) with 
smallest energy dependence allowed by Fig. 1 was used. We can see in Fig. 7 that the 
correction factor of the lowest energy part is still more than 150 % larger than that of 
the highest energy part. This fact was already taken into account in estimating the 
flexibility of the effective radii in § 3-1. From Fig. 7 we get 


MIG=t or)’. 
C,(W,) /C, V2) = (22a) 
1.22 fot VA, 


cand 
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Fig. 7. Corrected correction factors C,(W7) which are obtained by subtracting 
possible maximum effects of (theoretical) correction terms from the experimental 
one. For detail see 3 4. 


1.54 for ST, 
C,(W,) /C,W,) = (22b) 
1.45 for VA. 
The right-hand sides of (22a) and (22b) are valid only as the lower limits. 

Using (22a) and (22b) we can calculate in the same way as in § 3.2 and § 3.3. 
The limitations for |a@s,| and |ay,4| are summarized in table II for several values of 
(,/P which cover possible variation of ,, (12b). The restrictions are somewhat weakened 
from those in Table I. However, the interference terms must be still about maximum. 

Although we made somewhat doubtful assumptions in estimating the upper bounds 
of the nuclear matrix elements in the correction terms, it will be well compensated by 


the assumption that all the matrix elements take the possible maximum values simul- 


=) 
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Table II. Restricions for degrees of interferences between S and T, |asr|, and between V and A, |ay-4|, 
which are obtained from analyses of RaE spectrum with the effective nuclear radii theory including possible 
maximum effect of correction terms. Only p,/o values satisfying (12b) can actually occur. 


type | Syl VA 
i Se | a ee oer eee ae ee 
Eq. used (22a) (22b) (22a) (22b) 
0;/0 | lasr|> lasr| > laral> lev a|> 
0.4 0.997 | 0.998 | 0.998 0.998 
0.6 | 0.995 0.995 0.996 0.996 
| 

1.25 0.982 | 0.985 0.987 0.987 
2 0.961 | 0.968 0.973 | 0.974 
3 | 0.911 | 0.942 0.950 0.950 


taneously and are added constructively. Suppose that if this occurs for the ordinary 
first forbidden matrix elements the log ft is about 4.2 for a J—>J+0 transition, which 
corresponds to a fairly fast allowed transition. Thus we can regard the values in Table 


II as giving safe limits of interference terms. 


§ 5. Discussions and conclusions 


In this paper we have assumed the following points : 

a) The experimental correction factor for the /-spectrum of RaE which was measured 
by Plassman and Langer’. 

b) Validity of the generalized Fermi theory for /-decay, in which neither derivative nor 
non-local interaction is included. 

c) Applicability of the “effective nuclear radii theory’ with some correction terms. 
d) Validity of the time reversal invariance of the Wigner type for the strong nuclear 
interactions. 

We shall discuss the validity of these assumptions in the following. 

Our results will not be changed significantly even if we use other experimental data. 
Almost the same spectrum were obtained by the independent experiment of Wu and 
Lidofsky™”. 

We have no definite evidence that any theory of /-decay other than the Fermi 
theory is successful. The derivative interactions as proposed by Konopinski and Uhlen- 
beck”) might be successful for explaining the RaE spectrum, but contradict with 3-spectra 
of most of other nuclei. 

In connection with the assumption c), we should mention that the radiative correc- 
tion” is omitted in this paper. Although this correction seems to be negligible from 
both the theoretical arguments and the absence of positive experimental proof of its exis- 
tence in any (-decay, we cannot make definite statement about it. Also we cannot be 
completely sure that our estimation of the upper bounds of the nuclear matrix elements 


are all correct, though the failure of our argument about it seems to be highly im- 
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Fig. 7. Corrected correction factors C,(W7) which are obtained by subtracting 


possible maximum effects of (theoretical) correction terms from the experimental 
one. For detail see 24. 


1.54 for ST, 
C, (W,)/C, W,) = (22b) 
1.45 for VA. 
The right-hand sides of (22a) and (22b) are valid only as the lower limits. 

Using (22a) and (22b) we can calculate in the same way as in § 3.2 and § 3.3. 
The limitations for |a@gs,| and |ay4| are summarized in table II for several values of 
(,/? which cover possible variation of ,, (12b). The restrictions are somewhat weakened 
from those in Table I. However, the interference terms must be still about maximum. 

Although we made somewhat doubtful assumptions in estimating the upper bounds 
of the nuclear matrix elements in the correction terms, it will be well compensated by 


the assumption that all the matrix elements take the possible maximum values simul- 


1 


* 
- 
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Table II. Restricions for degrees of interferences between S and T, |as7|, and between V and A, lay 4|, 
which are obtained from analyses of RaE spectrum with the effective nuclear radii theory including possible 
maximum effect of correction terms. Only ,/p values satisfying (12b) can actually occur. 


taneously and are added constructively. Suppose that if this occurs for the ordinary 
first forbidden matrix elements the log ft is about 4.2 for a J]->J0 transition, which 
corresponds to a fairly fast allowed transition. Thus we can regard the values in Table 


II as giving safe limits of interference terms. 


$5. Discussions and conclusions 


In this paper we have assumed the following points : 
a) The experimental correction factor for the §-spectrum of RaE which was measured 
by Plassman and Langer’”’ 
b) Validity of the generalized Fermi theory for /-decay, in which neither derivative nor 
non-local interaction is included. 
c) Applicability of the 
d) Validity of the time reversal invariance of the Wigner type for the strong nuclear 


° a0 17 < . 
“effective nuclear radii theory” with some correction terms. 


interactions. 
We shall discuss the validity of these assumptions in the following. 

Our results will not be changed significantly even if we use other experimental data. 
Almost the same spectrum were obtained by the independent experiment of Wu and 
Lidofsky””. 

We have no definite evidence that any theory of /-decay other than the Fermi 
theory is successful. The derivative interactions as proposed by Konopinski and Uhlen- 
beck®) might be successful for explaining the RaE spectrum, but contradict with ,3-spectra 
of most of other nuclei. 

In connection with the assumption c), we should mention that the radiative correc- 
tion® is omitted in this paper. Although this correction seems to be negligible from 
both the theoretical arguments and the absence of positive experimental proof of its exis- 
tence in any [-decay, we cannot make definite statement about it. Also we cannot be 
completely sure that our estimation of the upper bounds of the nuclear matrix elements 


are all correct, though the failure of our argument about it seems to be highly im- 
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probable. It should be noted that the mesic corrections”” can be regarded as being 


included in our theory, because the nuclear matrix elements are treated as adjustable para- 


meters and their upper bounds are determined on the empirical basis. 

From the assumptions a), b), ¢) and d) we are led to the following conclusions : 
a) In the case of SVTA mixing, @sp=@y,—=0 can be definitely excluded (see § 3-1). 
b) For ST mixing 0.97;< |@s7|<1 (see Table Il). 

c) For VA mixing 0.97)< |ay,4|<<1 (see Table II). 

d) SVTA mixing seems to be unfavored from the view point of probabilism, although 
the RaH spectrum can be explained in principle if |agp|1 and |ay4\~1 (see § 3-4). 
These conclusions are consistent with most of other recent data” on the interference 
terms. Among all, our conditions are the most stringent and strongly suggest the actual 
occurrence of |@sp|=1 or |ay4|=1. The equations \@sp|=1 and |@y4\=1 are equivalent 
to eqs. (3), from which follow the validities of the two component neutrino theory” 
and the time reversal invariance” in /-decays. 

Similar results were obtained by Yoshizawa et al2® with similar analyses on Cs 
decay. 

The authors express their sincere thanks to Dr. M. Morita and Dr. T. Kotani for 
repeated correspondences. They are also indebted to Dr. H. Takebe for valuable discus- 
sions, who has investigated the RaE spectrum using his theory” mut 
(MY) thanks Dr. K. Way for her kind hospitality. 
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In this paper the mesic corrections to the beta-decay interactions in a nucleus are investigated. 
It is found that the Gamow-Teller coupling constant should apparently become about 8% smaller 
by exchanging pions among nucleons. Also the mesic effect in the Lforbidden case is discussed. 


§ 1. Intreduction 


It is known” that the strong interaction between a nucleon and pions somewhat 
changes the type of the beta-decay interaction of the nucleon. Although for the beta- 
decay of a bare nucleon Gamow-Teller and Fermi coupling constants have just the same 
magnitude, that is, |Jer/ g'p|=1, it is expected that it deviates from unity owing to the 
mesic corrections. In fact, this theoretical expectation,” the apparent increase of the 
tatio, seems to be consistent with the empirical fact,” that |Jer/Jr\?=1.3720%, though 
the extent of the agreement is not clear, because the theory involves some ambiguities. 

However, all the previous investigations” are based on the assumption that the 
nucleon can be regarded as a free particle. Actual nucleon exists in a nucleus, and the 
meson cloud should be a little suppressed by the exclusion principle or changed by ex- 
changing mesons with other nucleons. If such effect is not negligible, the type of the 
beta-decay of a free neutron should be different from that of a neutron in a nucleus. 
The relative magnitude between the Gamow-Teller and Fermi coupling constants, which 
has been obtained by the analyses of ft-values of a neutron and the other light nuclei, 
probably includes both of the two types of the mesic corrections, one of which exists. 
even in the beta-decay of a free nucleon and the other comes from the many-body effect. 

Though it is very difficult to estimate the latter effect quantitatively as well as the 
former, one can obtain an estimate by noticing that the effects have an intimate relation 
to the quenching effect for the nuclear magnetic moment as shown in this paper. The 
corrections to the nuclear magnetic moment due to the interaction between the nucleon 


current and the electromagnetic field have recently been calculated,” using the dispersion 
relation and assuming that only one pion is exchanged. 


Analogously it is possible to. 
calculate the many-body corrections to the beta-interaction, also assuming that the one- 
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‘pion exchange makes a main contribution. This assumption implies that it is sufficient 
to replace the z-component of an isotopic spin operator, 7;, which appears in the calcula- 
tion of the nuclear magnetic moment, by the other component, 7, or 7 in order to 
obtain our results. 

The numerical results in section 3, in which we adopt the Fermi gas model as 


nuclear model, will be discussed in the last section. 


§2. The many-bedy mesic corrections and the exchange 
magnetic moment 


Feature of the mesic correction is that the lepton parts of beta-interaction do not 
suffer any effect, so that the relativistic transformation properties and selection rules of 
the nucleon parts are invariant. (It can never occur, for example, that an additional 
pseudoscalar interaction appears from the scalar interaction.) Furthermore, so far as the 
non-relativistic approximation for the nuclear matrix element can be allowed as the first 
approximation, there exists no difference between the mesic corrections to the vector and 
scalar interactions (or the axialvector and the tensor interaction) in the allowed beta- 
decay. 

In this paper we assume that the main part of the many-body effect comes from 
the one-pion exchange process as shown in Fig. la. From the general theorem” it can 
be proved that the effect of the exclusion principle to the meson cloud in the case of 
no-spin exchange can be legitimately calculated as the exchange terms of the higher-order 
process such as Fig. la. For simplicity we assume the nucleus to be a closed core and 
a freely moving nucleon around it. Although the first assumption is sufficiently good, 
our results may be changed 
considerably if we take into 
account the possible strong two- 
body correlation in the realistic 
wave function of the nucleus. EE Bee 
Therefore the numerical results 
of next section should not be 
taken too seriously, but should 
be understood qualitatively. 


a ie Fig. la Fig. 1b 
The nonrelativistic nuclear 8 & 


matrix elements for the allowed beta-decay are conventionally written as i (Fermi type) 
and |7.0 (Gamow-Teller type). If we replace 7s in the vertex operator of Fig. la by 


T,, we obtain the corresponding process as shown in Fig. 1b, which represents the exchange 
nuclear magnetic moment due to the electromagnetic interaction with the nucleon current. 
Recently Fig. 1b has been studied by Hara and one of the authors,” who related the 


electromagnetic interaction to the pion-nucleon scattering data using the dispersion relation 


and assuming the validity of the static approximation. As will be noticed at once, only 
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the Gamow-Teller nuclear matrix element \re9 is corrected by the mesic effect in this: 


cannot lead the nucleon to the isobaric 


paper, because the Fermi matrix element, | 7z, 


“ 


state of total spin 3/2 and isospin 3/2, which is prevalent in the cross section of the 


(m, 7) scattering or (7, 7) production. 
§ 3. Detail of the calculation 


By the same method as in reference 3 the correction terms corresponding to Fig. la 


to the transition operators of the Gamow-Teller type in the allowed beta-decay are written. 


S 
m= > My 
DP = (+r) {POOF — (a 4-0 | file) Be 
MO = (+2) 4 (+0) fie) My 
MP = (fw) 4 PEE EO) — (ot) | fala) We 
MP= (e— 7) 4 (G0) f(x) Me 
MP= (ex), {PO z x9) _ (atx) f(x) M 
and 
MO= (eX 7) (XO) f, %) Di, (1) 


where the lepton part 


Mi =1d.7s7u Pr (axialvector) or Det tury Pr (tensor) , 


3 re 1 Neo 


< x 


and fi® =fa() =2fh(~) =— (2 197) gA(0)( . ae 


fi) =fu®) =fo) = — (2/157) 94 (0) ~ (2) 
(x=pr: 2: the pion rest mass). 


We put h=c=y=1 throughout this paper. The value of A for zero-point energy, 4(0), 
has relation to the total cross section, 7,, of the pion scattering in the pure (3/2, 3/2) 


state as follows : 


tes] 
y 


A(0) = (5/187) | (dp/a,”) o33=0.33 (3) 


0 


where the small contributions from o,, and o,, are neglected. 
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For simplicity we first assume for the nucleus to consist of a closed core plus an 
odd nucleon, only the latter of which takes the responsibility of the beta-transition. Then 
the transition matrix element due to the many-body effect is given by 

Gy MI) = 33) Gi) Ke) FY) H@) 


X MP (1) 4.) dx,dx,, ie 


where the transition occurs between the states, ¢, and 4/,’, and i), represents a state in: 
the closed core. 

It is quite easy to estimate the direct terms of the right-hand of eq. (4), which 
turn out to be an order of magnitude larger than the exchange terms. The sum of the 
ordinary Gamow-Teller nuclear matrix elements and the effective one-body corrections. 


arising from the direct terms are expressed numerically as : 


e 


Mora +t Man= 9 (1 — 0.084) \ Bale (5) 


Eq. (5) shows that in the Fermi gas model the Gamow-Teller matrix element is quenched. 
by about 8% owing to the one-pion exchange effect. 

For the exchange terms of eg. (4) the numerical calculations can be carried out in 
a tedious but straightforward way. Then we can express the result as 


Pe 


Wexcn= — 0.00849 | t, (3k (-k) /Pr—9), (6) 


where we assumed, for convenience of the calculation, that the magnitude of the 
momentum, k, of the outside nucleon is just equal to the Fermi momentum, py. If 
we take off this simplifying assumption, the numerical factor 0.0084 will be slightly 
modified. Smallness of We.xcn means that the previous treatment for the mesic effect is. 
almost unchanged, although we take into account the Pauli principle. Apparently Wexcn 
depends on the property of the odd nucleon state. To be acquainted with the dependence 
on the angular momentum of the last nucleon, we calculate the Gamow-Teller nuclear 
matrix element explicitly in the case where the last nucleon has the total angular momentum 


j and the orbital angular momentum il. 


qj'll We | D a (jl Mora Mair + Meexen | D) > (7) 
Cj" || Mall 7) = — 0.084 (7"|| Movall {5 (8) 
(2j—1)/2(j+1) for j=l+4 


Gil Meexen ll /) = 0.0084 Cj | Dtorallf) (9). 


(2j +3) /2j for (=. 


and 
(j-+ 11] Mexenll 7) = — Cjl] Meceall j+1) = — 9.0042 F411] Morall j) for F=j+4 
(10) 


where (j/||Mtl|j)’s are the reduced matrix elements. Table 1 shows the numerical results 
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; i i model 
for the total corrections to the ordinary nuclear matrix element in the shell 5 


(jl! Morall 7) d 


Table I Values of 4’s, which is defined by (j/| Mtl) =A—A4 ¢ 7 WWQoralli) 


ae é 51/2 | p3/2 P1/2 dsy2 dsrs 
— | 
si/2 8.4 (%) | =a a wa "a 
p32 | = 8.1 8.8 — = 
pile — | 8.8 5.0 
d5/2 “5 — ) — 7.9 5 
dj2 — — | — 8.8 6.7 
gfe, Grek WEE Ln ee) ee eS eee 


It is also interesting to show to how much extent the beta-tramsition can occur in 
the /-forbidden case by virtue of the exchange terms. Similar calculations to the case 
of eq. (9) or (10) lead to 


(j+1||Mexen|| 7) = —0.0084 3, Gry ?) (1) 
where 
jHlit3 
j+1=L,—-32, 
if we put the initial and final orbital angular momentum to be l, and /,. It is clear 


that the above matrix element is about 10° times small compared with the ordinary 
superallowed transition. There are a number of allowed transitions having anomalously 
long lives, of which the most remarkable example is C¥%—»N“(0* 17) with log ft=9.0, 
in spite of the fact that log ft= 2.9. for He®—>Li®(0* >1*). Though some of these 
transitions may be caused mainly by such mesic effect, it is quite difficult to decide whether 


the transition arises from the mesic effect or the mixing of the impurity states. 


$4. Results and discussions 


In this paper we have studied the many-body mesic corrections to the allowed beta- 
transition. The main results are the following two points : 

i) The numerical results in the one-pion exchange and the Fermi gas model show 
that the effective Gamow-Teller coupling constant in heavy nuclei should be about 8% 
smaller than in the free nucleon, whereas Fermi’s one is not changed. It is already 
known that the mesic corrections in a free nucleon lead to the apparent increase of the 
Gamow-Teller coupling constant relative to Fermi’s one. Ina nucleus the above tendency 
should be somewhat mitigated by the many-body effect. Namely 9,7 should be on the 
average a slightly decreasing function of mass number. The future precise measurement 


with the ft-values of neutron and the other light nuclei will give an insight into this 
problem. 


in eS 
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ii) The exchange terms, which make the allowed-type beta-transition possible also 
in the /-forbidden cases, turn out to be small as shown in (6) and (11). Probably 
the discrimination between the mesic effect and the configuration mixing is impossible. 
However, another aspect of the two-body mesic effect is that the two-particle transition 
can occur in addition to the one-particle transition, if we stand on the view-point of the 
independent particle model. The order of magnitude of the former type transition 
probability is maybe the same as 6.4X10~* times the superallowed transition. It seems 
to be possible for some of the allowed unfavoured and U-forbidden transition to be ex- 
plained by the two-particle transition. 

Of course the numerical values in this paper should not be understood quantitatively, 
since the following additional effects must be also taken into account : 

i) The expectation value of the two-body operator is quite sensitive to the existence 
of the two-body correlation in the nuclear wave function. For example, the deuteron 
model probably enlarges our effect. 

ii) Validity of the static approximation is quite doubtful if two nucleons approach 
too near. In reference 3 the detailed calculations were carried out, the higher Fourier 
components being cut off, but the order of magnitude of the numerical results is invariant. 

iii) As the next approximation the two-pion effect should be taken into account. 
The similar problem with the exchange magnetic moment has been treated” by the 
fourth-order perturbation method, which : ‘ 
suggests that the two-pion exchange effect | 
is not so large compared with the one-pion \ 


exchange so far as we take the cut-off 


momentum not so high. | Tr 
iv) Considerable contribution to the 
nuclear magnetic moment comes from the | 7 
. > . * 
ts’ type process as shown in Fig. 2a. o 
Cee sewd as é Fig. 2a Fig. 2b 


If the following process, 
m+ 357° +etty, 


exists even in a slight amount,” although the 7—z decay is not observed by virtue of 


the smallness of the z*—7z° mass difference, another many-body mesic effect of Villars’ 


type as shown in Fig. 2b will appear additionally. 
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Note added in proof 


After having written this short note, we received Nuclear Physics Vol. 6, in which we found a paper 
written by Bell and Blin-Style with respect to the mesonic effect. Their paper seems to be based on the 
idea similar to ours, but do not involve the theoretical treatment. Though the meson-theoretical work is 


referred to in it, it seems to be somewhat different treatment from ours. 
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A dynamical theory of nuclear distortion leading fission, is presented. This theory is based on 
the unified model (the collective model). But, since the spontaneous fission is, so to speak, an 
irreversible process, the present theory is treated reflecting this speciality. 

Starting from this dynamical theory, a semi-empirical formula of spontaneous fission half-lives of 
even-even nuclei is derived. This formula can well reproduce the general tendency of spontaneous 
fission half-lives of even-even nuclei. The essential features which affect half-lives are discussed. 

Finally, a reasoning for the asymmetry character of fission and the orientation in calculating the 


difference of half-lives between even-even and even-odd nuclei are given. 


§ 1. Introduction 


The theory of nuclear fission has developed in a close relationship with the theory 
of nuclear structure and nuclear reaction. Bohr and Wheeler’ attempted to disclose the 
mechanism of nuclear fission by making use of the liquid drop model. In those days, 
the liquid drop model of nuclear structure was regarded to be significant as a compatible 
model to the compound nucleus theory of nuclear reaction. Then it was considered that 
nuclear fission should be treated with the same model. 

In recent years, the nuclear shell model” succeeded in explaining the ground state 
spin and magnetic moment of nuclei. It was pointed out by Mayer” that the abundance 
of fission yield is closely related with magic numbers. 

More recently, it was shown that the unified model”? taking imto account both 
single particle and collective aspects of nucleus can explain the nuclear quadrupole moments 
and the properties of excited states of nuclei with great success. In this connection some 
authors» discussed the relation between the shell structure under the deformed 
nuclear surface and the fission barrier, and the stability of odd parity deformation. 

The success of the unified model in explaining nuclear structure, together with the 
recognition of well organized single particle modes in nuclear reaction” revealed the 
inadequacy of the liquid drop model. And it is becoming more and more apparent that 
the liquid drop medel is not efficient enough to explain nuclear fission. Inadequacies of 


applying this model to fission phenomena have been recognized especially in the following 


two aspects. This model cannot explain 
(1) asymmetric fission, and 
(2) systematics of spontaneous fission half-lives. 
Frankel and Metropolis’? showed that the liquid drop model with uniform density 
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; : i _ It can also be 
is more favourable to the symmetric fission than the asymmetric one 


shown that the asymmetric fission is probable under certain soar distribution ies in 
the liquid drop model." But, on this standpoint, it can be said that any liqui op 
model’? contradicts with the requirement of the parity conservation. Let us take an 
example. The ground state of even-even nuclei has even parity. If we wi Bi: 
the asymmetric fission of this type of nuclei with the liquid drop model, we are force 

to assume that odd parity states become mixing to even parity states in the course of 


fission. This is unreasonable. 


sh to explain 


Huizenga” pointed out two systematics in spontaneous fission half-lives : 

(2-1) half-lives of even-odd nuclei is longer than that of corresponding even-even 
nuclei, and 

(2+2) spontaneous fission half-lives of the even-even isotopes correlate with mass 


number regularly (see Fig. 2 of the present paper). These regularities are two-fold : 
(2-2a) for an isotope series of a given element, spontaneous fission half-lives of 
even-even nuclei go through a maximum with increasing mass number, and 
(2-26) roughly speaking, spontaneous fission half-lives < of the even-even nuclei 


correlate with mass number A by a simple relation, 
c7=c-exp[—kA4], 


where & is a constant independent of A and Z, and c is related to A and Z simply. 
The liquid drop model fails to explain these systematics. That is, this model predicts 


that the fission half-life becomes longer with increasing 4 for fixed Z. 0° This result 
contradicts with the experiments. 


Thus it is clear that the liquid drop medel is too simple to treat nuclear fission. 
On the other hand, it is apparent that the shell model is also too simple to treat 


nuclear fission, because the collective aspects of nuclear matter are essential in this 
phenomenon. 


On the contrary, one can understand at least qualitatively (1)? and (2-1), on 
the basis of the unified model. And one can expect to explain (2-2) on the same basis. 
Before entering into discussion about (2-2), however, it would be necessary for us to 
notice that the fission phenomenon includes two sorts of problems. One is a sort of 
problems which can be treated well by the kinematical considerations. Angular distribu- 
tions of photo-fission' and induced fission belong to this sort. The other is a sort of 
problems which require the dynamical considerations. Explanation of features in spontaneous 
fission belongs to this sort, because this problem is essentially a dynamical one, being 
concerned in an irreversible phenomenon. In this view-point, several attempts so far 
presented,” ~” which were based on, so to speak, the “ static”? considerations, seem to be 
unsuitable in applying to the spontaneous fission half-lives. This will be criticized in 
more detail in the last section of the present paper. Rather, we think that the problem 
should be treated with the time dependent consideration, as was done by Hill and 
Wheeler? (their Appendix Fig. 34). It was these circumstances which has motivated the 


present work. In the present paper, consequently, our attention is mainly concentrated 
to (2) and especially (2-2). 
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In § 2, a dynamical theory of nuclear distortion leading to fission is presented. In 
this section, first we shall discuss the reason why the consideration presented by Hill 
and Wheeler, is applicable to the half-lives of spontaneous fission. The derivation of the 
equation of Hill-Wheeler is made in Appendix of the present paper. Secondly, the 
consequences of this equation in applying to spontaneous fission half-lives are presented. 

Developing the discussion made in § 2, a semi-empirical formula of spontaneous fission 
half-lives of even-even nucleus is derived in § 3. This formula can reproduce the general 
tendency of spontaneous fission half-lives of even-even nuclei well, as will be seen in 
Fig. 2 and Table 1. The essential features which affect the half-lives are discussed in 
§ 4. Finally a reasoning on the asymmetry character of fission and the orientation in 


calculating the difference of half-lives between even-even and even-odd nuclei are given. 


§2. Physical image of the mechanism for spontaneous fission 


In this section, we will discuss what plays the important role in determining the 
spontaneous fission half-lives. This argument leads us to use the equation of Hill- 
Wheeler. For a while, we will mainly 
argue the problem of even-odd nuclei for E(a) 
simplicity. But the same argument is also 
applicable to even-even nuclei. This will 
be discussed in more detail. Then, we 
will derive the expression giving the half- 
lives of spontaneous fission for even-even 
nuclei. (eq. (11)). 

In the even-odd nuclei, it is confirmed 
that the picture of the single particle B 
moving in a deformed potential well is awa 
nicely fitted to the experimental result of 
nuclear spectroscopy.” Schematically, the les 
energy surface for the single particle in a 
deformed potential is represented in Fig. 1. se a 
Considering time dependently, the Fig. 1. Schematic representation of the effect of 
deformation leading to fission upon single particle 


iti ili at a crossin oint ices 
transition probability perckoa ee energy levels. ay indecates the equilibrium defor- 


between two single particle ani =o surfaces, mation, where the total energy is minimum with 
a and b, is decided by three factors”; the respect to w. 
transition matrix element,, H,,, the time 
rate of deformation, &, and the difference in slope of the two unperturbed potential 
energy surfaces, d(E,—E,) /da. 

The jump probability, J.,, is very small when H,,, is very weak, or when @ is very 
large, or when d(E,—E,) /da is very large, or when two or three of these are appreciable 
simultaneously. At some crossing points, H,, vanishes owing to the selection rule. @ 


is very large near the equilibrium deformation and is very small near the turning point 
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in semi-classical sense. About d(E,—E,) /da we will discuss later in this section in more 
detail. 

The jump probability, J,, becomes very large when the a 
constructively. And, at this point J,, is much larger than J,,. We name the crossing 
point satisfying the above condition the typical point aime 

Near the point P, the jump probability, Ja, will be effectively contributed by the 
coupling of two levels, E, and E,, though J,, may generally be influenced from many 
other levels. 

Following Hill-Wheeler’s expression, the jump probability, J.», is given as follows,” 


bove three factors contribute 


jst ae (1) 


where 


ae Aen” (2) 
had(E,—E,) /2da 

It is expected that the deformation corresponding to P will not be larger than the 
one corresponding to the fission threshold in a semi-classical sense. However, after pass- 
ing the typical crossing point P once, the single particle energy E(@) is fed into the 
surface oscillation energy. Then @& becomes larger and larger, the jump probability becomes 
smaller and smaller at further crossing points and the deformation proceeds on. Conse- 
quently, the half-lives of spontaneous fissions seem to be decided by the jump probability, 
Jay, at the typical point P. Hence the half-lives are given by 

t= fe. (3) 

In the above argument, we assume the single extra-particle for simplicity. This 
situation, that the single particle moves in a deformed potential, is justified by the recent 
experiment." This guarantees the validity of the adiabatic approximation adopted in the 
above arguments. 

In the case of two extra particles, there are no direct experimental proofs of the 
validity of the adiabatic approximation as yet. But under the following situation, we 
‘could expect the applicability of the same approximation. According to the shell model, 
the nuclei in this region are characterized by the large j, the spin of the single particle 
-state. Following Mayer’s argument,” the pairing energy for the extra particles is large 
for large j. This will cause the situation that the direct coupling between extra particles 
-over-whelms the particle-surface coupling and consequently the extra particles are coupled 
to the surface in the same manner as a single particle (BM**) II, c)). 

Hence the argument of this section will be applied without a serious modification 
also to the even-even nuclei. That is, the above discussion does not set the limits to 
the even-oddness of the nuclei. 

In what follows, we limit our attention to the even-even nuclei, since in this sort 
-of nuclei slow and regular variations of the various experimental quantities are expected 


whereas in the even-odd nuclei abrupt changes of these quantities are expected from 
nucleus to nucleus.!” 
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In the region considered, |H,,,|? also may not depend so much on mass number, 
etc. Since in the very neighborhood of equilibrium @ is very large, J.,, becomes so small 
that the slippage does not occur effectively. In other words, & is very small and it is 
taken to be constant with mass number etc., at the crossing point slipping effectively, 
irrespective of the order of magnitude of the zero point energy of nuclear surface oscilla- 
tion. Accordingly, the probability of slippage is affected only by the difference between 
the slopes of energy surfaces E,(a@) and E,(a@). This can be estimated roughly as 
follows. Expanding E(@) to the second order of a, we put 


E(a) =E©+E®.a+tE®.@?, a 


As for the mass number dependence of E™, we adopt the one given by Bohr and 
Wheeler.” 


hen, E® oc A**. (1—x), 
= ZIM / RY Aran 5 
(Z*/A) crit, = 44.65 (using the value given by Fermi’). 


And then, we may assume the following relation. 
E,(@) —E,(a@) =p (A) - (K+ K®-a+K®-a A’¥. (1—x)), (6) 


where K, K“ are constants, i.e., independent of the mass number, and (A) is the 
function of A. Assuming the nucleus as the potantial well with finite depth, ¢(A) can 
be found to be a decreasing function of A, since the difference in the slopes of the 
energy levels of unperturbed single particle states become smaller with increasing mass 
number. And K is not always smaller than K™ since K™ and K® correspond to the 
difference between E“ and E{?, and E® and E§, respectively, though E™ itself is 
generally much smaller than E® in the neighbourhood of the closed shell. The region 
under consideration is far from the closed shell, so E“ does not vanish. 

In order to estimate @ in eq. (6) from experimental quantity, we use the relation 
E,,.0c1/J, where E,,,, is the energy of zero point oscillation. And J is the moment 
of inertia which is related to a as JocM,-a’, where, M, is the effective mass for nuclear 
capillary oscillation. 

Eee (172) “(the experimental value of the energy of the first excited state). 


can be estimated from the equation 


Hence we get 

a= Ky/1/ME,,=K/A- E20, 0) 
where K and K’ ate constants with the change of mass number. Expanding eq. (2) to 
the power of a, we get, 

A,» 1 a) K® 2134 — 8) 

Gis | ob jee (KD EAK Ora AAA, x)\)). ( 
~ ba Vp (A) 

Using eqs. (1) and (3), 
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log,,= = log io (1 me a 


~ ene (9) 


Because J,y> Jaa, G is very large. Therefore the above expression is appropriate. 
Further in expanding the last expression, we take some G,, expand it to the first order 
of G—G, and get 

log, t= A=B-G, (10) 


where A and B are some constants. G, is the average value of G at the typical crossing 
point over the nuclei under consideration. Combining eqs. (8) and (10), we get finally 


log,,7 as follows. 
log, T=U—B" - [9 (A) 7°? + C[p (A) J7°”?- ATM (A 2) EZ, (11) 


where, 2, QB’ and © are positive constants and determined approximately as adjustable 
parameters with experimental values of known three pairs of half-life and the energy of 
the first excited state. Now there are 13 nuclei for which both the half-life of the 
spontaneous fission and the energy of the first excited state are measured and so these 


13 nuclei are studied by means of three adjustable parameters. 


§ 3. Derivation semi-empirical formula 


W. J. Swiatecki™ gave the semi-empirical formula giving the fission half-lives in 
connection with mass number, but his formula does not seem to stand on a sufficient 
theoretical basis. In this section we will derive the semi-empirical formula following the 
arguments developed in Section 2. 

For (A) of eq. (11), it is tempting to assume a decreasing function of A, having 
a form of p(A)=7°/A? (7; constant). This would not be so unreasonable in the 
region considered here. Then, eq. (11) becomes 


log,,T=AU—BA+CA"- (1—x) -EZY, (12) 
where B= -7. 


The above form well indicates the typical features of the experimental results of 
spontaneous fission half-lives. In the case that the second term is large compared with 
the third term, —really it will be shown in later numerical calculation that the ratio of 
the third term to the second is about 15/100, —the main features of the half-lives are 
determined by the first and the second term. And the third term appears as the correc- 
tion to the formers. Actually, it is well known that the relation 


logpt=U—B-A (13) 


represents the gross behavior of the mass number dependence of spontaneous fission half- 
lives. This is the relation which was stated in (2-26) of the introduction of this paper. 
The experimental values of energies of first excited states, being adopted here as the 


measure of the zero point energy, have some regularity, as is well known.” And we 
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can find that this regularity log, st (years) 
has some correlation with ; 
experimental data for half- 
lives of even-even nuclei, that 
is, in U and Pu the energy 
of the first excited state goes 
through a minimum with 
increasing mass number for 
given Z. The nuclei show- 
ing the above minimum 
correspond to those having °* 
the maximum point in the 
curves of spontaneous fission 
for even-even nuclei (see 
Pigs 2s 

Adopting 29{=102, B= 
0.441 and ©=5.76 for A, 
% and ©, respectively (E,,,,. 


is measured in unit of kev.), 


230NNos4 238 242 246 250 254 258 
the result is given in Fig. 2 Mocs 


and Table 1. In this result, 


the calculated values of Th fission. The dotted curve represents experimental results. The 
and U are rather small in solid curve represents the calculated value using eq. (13). 


number A 


Fig. 2. Mass number dependence of half-lives in spontaneous 


comparison with the experi- 
mental results. However, we may conclude that this result represents well the general 


tendency of experimental data. 


§ 4. Discussion and concluding remarks 


In spite of much simplification, our model seems to be able to reproduce the complex 
experimental results consistently. If we want to study this problem in more detail, we 
have to estimate |H,,|*, and the contribution from lower levels. At present, we have 
not enough knowledge to calculate these factors in detail. These knowledges should 
include more reasonable estimate of effective mass of capillary oscillation than that 
predicted by the present liquid drop model. Further, it is necessary to have more infor- 
mation about the structure of matrix elements and so on. The factor (A) in the 
preceding section must be estimated by using the deformed potential with finite depth, 
though we assumed here an inverse square proportionality with mass number for it. 

However, taking these into account makes not only the calculation much complicated 
but also the insight ambiguous. As a result it is considered that the conclusions in § 2 
and § 3 indicate the essential quality of this phenomena. Namely the half-lives of the 


spontaneous fission are decided by the probability of slippage in the crossing point of 
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Table 1. Hialf-lives of spontaneous fission 


| —_ er logjat (years) 


4 | 4 | ae . “Cousin oe Experimental values 
90 | 230 | 52.4 16.3 17.2 
232 49.0 16.0 18.2 
92 | 232 | 47.0 14.1 13.9 
234 43.5 14.4 16.3 
| 236 45.2 13.9 16.3 
| 238 | 50.0 12.9 15.9 
94 | 238 44.1 11.0 10.7 
| 240 42.9 , 1a 11.1 
| 242 | 44.6 10.5 10.9 
96 242 42.3 7.9 6.86 
| 246 / 42.9 7. 748 
| 248 | 43.4 4a 6.62 
98 | 250 | 42.0 4.1 4.18 


E; the energy of the first excited state. 
For Th, I. Perlman and J. O. Rasmussen, Alpha Radioactivity, Encyclopedia of Physics vol. XLII, 
Nuclear reactions III (1957), Springer-Verlag (Berlin). 
For U, P. C. Allen, Phys. Rev. 105 (1957), 1796. 
For other nuclei, E. K. Hyde and G. T. Seaborg, The transuranium Elements, Encyclopedia of 
Physics, Vol. XLII, Nuclear reactions III (1957), Springer-Verlag (Berlin). 

t; half-lives of spontaneous fission. 
For Pu2!2, Mech et. al., Phys. Rev. 108 (1956), 340. 
For Cm28, Butler et. al., Phys. Rev. 108 (1956), 965. 
For others, those values are used which are collected by A. Ghiorso and cited by W. J. Swiatecki 
(Phys. Rev. 100 (1955), 937). These values are in accord with E. K. Hyde and G. T. Seaborg 
(I.c. p. 306) except Th***, 


‘the energy levels, at which the three factors H,,, & and d(E,—E,)/da appearing in 
(2) work constructively, and especially influenced by the difference in the slopes of the 
energy levels with reference to the deformation. ‘5 

Our treatment is different from those attempts”~’ to calculate the potential barrier 
for fission by the static considerations and to estimate the penetration probability for it. 
In these trials, the probability of slippage in crossing points is always assumed to be 
unity so far as the whole levels are closely occupied by particles in order from the bottom 
-of the potential well. In these models, it is assumed that the unoccupied levels never 
come down below the uppermost occupied level in any stage of the deformation. But 
such an assumption cannot be justified and within such a restriction we may not expect 
to obtain good results in calculating spontaneous fission half lives. Naturally, at the 
initial stage of the transition probability J, from a level 4 to another level 6 varies with 
energies of zero point oscillation, which show some regular variation upon mass number 
‘as remarked in the previous section. Consequently, J,,, should not always be equal to 
unity. In our treatment, on the contrary, this situation is taken into account properly. 
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Finally, we will discuss the asymmetry character in nuclear fission and the difference 
of half-lives between even-even nuclei and even-odd nuclei. First we shall discuss the 
‘asymmetry character. The level distance between the closed shell and the upper levels 
of the unperturbed (i.e. not so deformed) nucleus is very large. Consequently, the 
crossing points of closed shell and the upper levels appear in far distorted form compared 
to the initial distortion. This circumstance gives rise to the following two situations. 
One is that the number of crossing points near the initial stage of the distortion will 
be much reduced. Hence, at the decisive stage of spontaneous fission which is expected 
to lie not so far from the initial distortion, the closed shell nature is not so disturbed. 
The other feature is that near the closed shell the difference of slopes of concerning 
levels is large in the neighborhocd of the crossing points. This makes the transition 
probability J, small at these points in this region compared with that at crossing points 
in other region. As a result, the levels below closed shell tend to conserve the single 
particle nature in deformed potential. That is, even in the middle stage of the process 
these levels may conserve the closed shell nature fairly well. If the collective oscillation 
of odd parity is raised in this stage, as a result of excitation of several particles coupled 
together with odd parity, the asymmetric character in fission will appear. Hereby asym- 
metric fission and closed shell nature will be strongly correlated. In such a stage, how- 
ever, we must take into account many concerned levels and the adiabatic approximation 
may not be justified. In this case, the quantitative treatment of fission asymmetry have 
much complexity. 

As Newton’ pointed, it may be considered that the difference in half-lives between 
the even-even nuclei and even-cdd nuclei reflects the following circumstances. In heavy 
nuclei, generally the total angular momentum of single particle is so large that the pairing 
energy is also large. And two particles are easy to couple together. Then, in even-even 
nuclei each particle is able to jump over from the initial level to the next level with 
either the same parity or the different parity from the initial one. On the other hand, 
in even-odd nuclei, when one particle outside the core carries the parity of the nucleus, 
the requirement of parity conservation forbids this particular particle to move to another 
level with different parity. Of course, this argument is more complicated when several 
particles are coupled together and the states resulting the coupling of these particles 
include both even and odd parity. However, as is mentioned in § 2, it is confirmed 
experimentally’ that the single particle moving in a deformed potential is the suitable 
model in even-odd nuclei, this argument is to be taken as reasonable. According to the 


above considerations we can understand that the lives of even-odd nuclei are longer than 


those of even-even nuclei. A corresponding character of even-odd nuclei is found in 


alpha decay.”? Quantitative discussion of the problem requires, however, the more 


detailed knowledge of matrix element H,,. 
In conclusion, authors express their sincere gratitude to Professor M. Kobayasi for 


his kind advices throughout this work. 
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Appendix 


In this appendix, we will derive Hill-Wheeler’s equation (Fig. 34 in their Appendix) ,” 
starting from the total Hamiltonian of the system. 

Now we consider the point, Q where two single particle energy surfaces (a and 6 in 
Fig. 3) are crossing each other. Here the contributions from other levels (c in Fig. 3) 
are assumed to be neglected. 

Now, the Hamiltonian of the total system is given as follows, 


H=H,(*) +Him a) +H(a) +H'(s, a). (A-1) 


Here, 
H, (x) + Hine (x, &) +H (@) (A- 2) 


is the Hamiltonian proposed by A. Bohr,” a 
where H,, is the Hamiltonian of an extra Ee) 
particle, H(a@) is that of the surface oscilla- 
tion, including up to the square of a, and Q 
Hi,(x, @) is that of the interaction between 
the single particle and the surface oscillation. 


H!(x, a) represents the residual interactions 


including all the remaining part, which was i Se 


neglected by A. Bohr. H’ (x, a) includes 3 ce : 


the third and higher powers of a, besides 
the interactions between nucleons probably, 
which could be expressed as the body oscil- ; 
P ; Y Fig. 3. Schematic representation of single particle 
lation, and so forth. Then, this part may energy surfaces. 
include other characteristic coordinates, for 
instance, the one describing the body oscillation mode, etc., but generally it can be 
written by using x and a’ because of the conservation of degree of freedom. 
Inclusion of H’(x, @) in (A-1) is necessary for a system which undergoes. 
spontaneous fission. The Schrodinger equation of the total system is 
oy 
ib —— Hae de A:3 
ry (A-3) 
We will solve the above equation with time dependent perturbation method, initial 
condition being 


Wf. @) -fa(a). (A-4) 


This initial condition is suitable under nearly adiabatic conditions. This is the so-called 
cos stationary state wave function,” which is already perturbed by H(a@), but @& 
is treated as though a were fixed. /, is the eigen-solution of the equation 


(Hy + Hin &)) Has 4) =Wa(@) Pal &). (A-5) 
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The functions ¢/,,(x, a) form an ortho-normal set with respect to x for all values of the 


parameter a. f,(@) is the eigen-solution of 
F(a) -fa(@) =&,: fa (Ee)- (A-6) 
In order: to solve eq. (A-3), the following interaction representation is adopted. 
Geen Mp tied art (yn a 2) 
Fie Ae Bint i2 | Pyle US pt Aino te) (A-7) 
Then eq. (A-3) is transformed as follows, 
ib (0/0t) P7=H" P!, (A-8) 
and the initial condition becomes as 
C't=—o)=f. (A-9) 
Near the crossing point Q, we can put 1 as 
Pi= $a Gala, t) +%r9o(@, 4), (A-10) 
where 
In (@, d= | on" (x, a) P(x, a, ddx. (A-11) 
Introducing eq. (9) into eq. (7), we get the simultaneous equations, 
16 9.=HaJat Har 
1D y= Fy Jat Air, (A-12) 
where 
Hel | dad (H(a) +H" (x, &)) Pp dx- ef Mm(@) —Wala))th, (A. 13) 
As for the diagonal part of eq. (A-12), we put 


Ayin= FAs (H(a) +H" (x, &)) pn dx=E,, (@) (A-14) 


and as for the non-diagonal part, we put 


yin = | db, * H! (x, a) p,dxei Wm(@) —Wn(@) tlh = H,,,. (A-15) 


mir 


Then eq. (A:12) with eqs. (A-14) and (A-15) is the equation of Hill and Wheeler. 
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——- Thermodynamical Properties* 
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An attempt is given for calculating the thermodynamical quantities with respect to the vibrational 
frequencies of the imperfect lattices. Using the method of contour integral the expressions for the 
difference of vibrational free energies between imperfect and perfect lattices are derived. These are 
given in two forms available at high and at low temperatures. 

As the examples of the practical calculations using these formulae, the n-dimensional simple lattices 
with defects are treated. Especially the simple cubic lattice as well as the linear chain are investigat- 
ed in detail. Although in this example a simplified model of the actual defect is considered, it may 
be enough to discuss the general features of the thermodynamical properties of the imperfect lattices. 


§ 1. Introduction 


The calculation of the free energy of the crystal with defects has been done for the 
Schottky or Frenkel defects, although it was only a phenomenological one. Recently the 
investigation of the influence of lattice imperfections such as impurities, holes, and dislo- 
cations on the physical properties of crystals seems to have been one of the subjects of 
great current interest in solid state physics. The thermodynamical properties of the 
crystal with defects also have been investigated by several authors.'~” 

First, Kirkwood and Stripp’’ calculated the influence of the Schottky defects on the 
vibrational part of the crystal partition function by means of perturbation theory. Their 
evaluation of the partition function seems, however, to be fairly approximate because only 
the first term of the expansion in terms of T~', T being the temperature, is retained. 
In the problem of the evaluation of the effect of defects on the vibrational partition 
function there may appear two difficulties: the first, the troublesome problem of the 
evaluation of the influence of defects on the frequency spectra of the lattice vibration 
and the second, the difficulty to carry out the summation over all the normal modes. 

The problem of the first kind has been treated by several authors”? '~' . Most of 


* This work was started while one of the authors (T. Tanaka) was at the University of Maryland 
and was completed after he returned to the Kyusyu University. Quite recently we learned that the same 
problem had been taken up by the Maryland group (J. Mahanty, A. Maradudin and G. Weiss) after T. 
Tanaka came back to Japan and that a similar result was obtained (Prog. Theor. Phys., this issue). Since 
the two groups had developed the formulation rather independently until they obtained the final results, 
it was hoped that the two results be published in the same place. 
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them have limited their considerations only to the one dimensional lattice, while the 


general theory developed by Montroll and Potts” is applicable to the case of three 
dimensional lattice and, actually, a general expression for the additive functions of the 
frequencies of normal mode vibrations of the crystal lattices has been given without 
perturbation. The second difficulty, however, seems to remain still unsolved. Montroll, 
Potts, and others calculated the zero-point energy by means of replacing the sum by 


the contour integral on the frequency-plane ; nevertheless, the difficulty of integration 


2,6,7) 


seems to have limited the range of possible calculation remarkably. 

On the contrary, in the case of the vibrational part of the crystalline free energy 
the method of contour integral seems to be more convenient than the usual perturbation 
because, fortunately, the special form of the function of the free energy (logarithmic 
form) makes possible the integration of the contour integral on the frequency-plane. 
(See § 2). Although our method is, of course, not applicable to the investigation as 
to the effect of defects om the other parts of the crystalline free energy, it is rather 
useful so far as the vibrational free energy is concerned. . 

In §2 the expression for the imperfect lattice is derived by the method of contour 
integral. Since this is of the integrated form, we can find the final result at once if 
the characteristic function—that is, the secular determinant for determining the normal 
frequencies of the imperfect lattice divided by that of the perfect one—is given in any 
way as the function of frequencies. This expression is, however, of the form that is 
‘convenient at moderately high temperatures. Therefore the expression available at low 
temperatures is derived next. This is of the form such that it consists of the change 
of the zero-point energy from the perfect lattice and the temperature dependent cor- 
rection terms against this. Then, the methods for calculating the characteristic functions 
are reviewed in §3. The method of Montroll and Potts” is employed in this paper for 
evaluating the characteristic functions. 

As the examples of the practical calculations using these formulae, the n-dimensional 
simple lattices with nearest neighbor interactions are investigated. The discussion is con- 
fined mainly to the case of impurities including isotopes. In addition only the harmonic 
terms are, for simplicity, retained in the potential, that is, the effects of thermal expansion 
and any relaxation of the structure about the defects are ignored. In §4 the change 
of the vibrational free energy due to one defect will be considered. The calculations 
are done for the cases of the simple cubic lattice, monatomic and diatomic linear chains. 
The interactions between defects and that between a defect and the boundaries are investi- 
gated in §5. Then, the influence of many disordered defects on the vibrational free 
energy is treated in § 6 and is shown that the method of Montroll and Potts” for de- 
termining the characteristic functions is also useful in some degree to the case of 
many disordered defects. 


§2. General formulae for the free energy of the imperfect lattice 


First, we consider the quantity : 
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G=Si9(o,) : (2-1) 


> . : 
where Y(w)’s are arbitrary functions of frequency w and the summation is extended over 
all the normal frequencies ws of the lattice vibrations. In the present paper we shall 
use, for convenience, as the frequency and temperature the reduced (dimensionless) 
quantities 
w=w/o, and T==27 (xT /bw;), (2-2) 


where w, is the maximum frequency of the perfect lattice and « is the Boltzmann 


constant. Since the normal frequencies w,’s are the solutions of the secular equation, 
D(o) = (const) -/] (w*—w,7) =0, (2-3) 
? 


we can write (2-1) in the form 
t ot : 
G—— 4 J (w)d [log D(w) | (2-4) 
271 cr 
as we get from (2.3) 


Alege Do se ee Bee (2-5) 


j or o+a,! 


In (2.4), the contour C’ is such that it includes all the positive zero points of D(w) 
in it and not any singular point of D(w) and g(w). Therefore, when the secular 
equation of normal frequencies of the perfect lattice is given in the form 


D,(w) = (const) - /] (w?—w,”) =0, (2-6) 
d 
we obtain as the difference quantity JG of G between the imperfect and perfect lattices : 


Jeane Ge. “4 7(o)4 | log Big) | (2-7) 
Ti Sg D,(@) 
In (2.7) the contour C is similar to the C’s in (2.4). Hence, if the function 9(w) 
is regular on the positive half w-plane we can take as C the following; that is, it con- 
sists of the positive semi-circle of radius R and the vertical diameter along the imagi- 
‘nary axis on the w-plane. Since d log(D(w)/D,(w)) is given as the difference of two 
quantities ; one is given by (2-5) and the other is the corresponding quantity of the 
perfect lattice, we can prove that the contribution of the semi-circle to the integral (2.6) 
becomes zero as R tends to infinity if the function g7(w) behaves~|w|*(2> X20) as 


|w|—> 00 ; namely, 


1 1 a 2W 
Cease -=———,, 
wo—W; Ow CO sin gees Dy 
9 092 

v 20 oO; —W,;- iL 

Pi EE EE ‘ ete 9 Sys 0: pe OAc 3 
p 2g : 2 p 2 2 

wo — 0; ww; (w —w) (w — (7 ) ow 
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} g(w) OC = Vdo~!2l 50 (if 2>XL0) as jw|>~. 

R w wo 

Then the contribution of the integral in (2-7) comes only from the integral along the 
imaginary axis if the function 7(w) has such a property. Furthermore, noticing that 
D(w) and D,(w) are the function of w®, d log(D(w) /D,(w)) is an odd function on 
the imaginary axis. We can, then, conclude that only the odd part of 7(w) contributes 
to the integral in (2.7) on the imaginary axis. That is, 


aaa eile ey d D(iw’) | F 
AG= | D ) load ! _— dw’. 2-8 
Ti Jo Lg Go’) ] ail ee, G's) 


If we put 7(w) =4hw, we obtain as the difference of the zero-point energy JE, 


4E,=— fen | aa [108 eo d ldo. (2-9) 
27 0 du! - D, (io) 


Eq. (2.9) is the same as that obtained by Montroll-Potts.”? (Hereafter we refer to it 


as (I).) 
Next, we consider the Helmholtz free energy of the system of harmonic oscillators. 


Since it is given by 


F=+3"ogl2 sinh( = aw; Is pei sete Ag (2-10) 


eed: 
Bs bo, T 


we may put 


J (w) = Flog 2 sinh a »)| : 
Then, since 


9 iw’) = : 1 log | 2 sin 70 )Aiz(k+ - ) \ w'=0 


for kn |w!| < h+1)z, 
ap ( )z, 


we obtain 


Lg (iw) baa = + iz (k+- 4) / B ; wo [0 


for kn< ol < (kR+1)2. 
lad (k+1) 


Inserting this into (2.8), we get 
4F=F—F, 


23 c d [ pin) 12, (A+R d D(iw! 
See lo do! —.S% | a lal? 
28 Jodo: © D, tio’) 1 Bass ie dae” Dey lt 


(2:13) 
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Eq. (2.11) is integrated at once and we obtain the final result : 


AF = A(0) —A(o0) +5 SVM AUT) —ACRE DT) (2-12) 
where 
Biss. D(io’) : 
lay) = eB io’) Gel)” (2-12a) 


Eq. (2.12) is of the form available at high temperature, because the second term becomes 
zero in the limit of high temperature and so can be regarded as the correction to the 


first term at high temperature. This can be easily seen by expanding D(w) in terms 
of D,(w). From (2.3) and (2.6) 


i (w’—w,”) 
I (w* —«;”) 


D(w) _ ns a wo; 
shee eae aes 22(of 04") — 


—+-+-]. (2-12b) 


= 


Since we may put w=ikT in order to insert D(w)/D,(w) into the second term of 
(2.12), |w| becomes large at high temperature and the second term of the above ex- 
pansion is of the order of 1/|w|*, the third term of the order of 1/|w|*, etc., which are 
of the order of 1/T*, 1/T’, etc., respectively. Then the second term becomes zero as 
T tends to infinity. 

Although (2.12) can be used at all temperatures, it is very inconvenient below 
ordinary temperatures ; because, at low temperatures, the zero point energy is dominant. 
Then it is desirable to rewrite (2.12) in the form in which the zero-point energy appears 
as the first term plus correction terms. Fortunately, this is done very easily using the 


Euler-Maclaurin summation formula : 


a 


Sig) =| 9 eee +9(c0)] 


+3i—1y fg) 9" (C0)}, 13) 


where B, is the Bernoulli number of the /-th order. If we put 


g (k) =k A(kT) —A((k+1)T) | 


in (2.13) and perform the partial integral,” we obtain 


AF = AE, + 2825) T® a,(T) (2-14) 


20 i=! 


* In the derivation of the expression (2-14) we regard the quantity A(AT) as that which is normalized 
so that the constant factor of (2:12b) is 1. The expression (2:14), however, can be used for the quantity 
A(kT) which is unnormalized. The authors are indebted to Mr. I. Mannari of Kyoto University who 


pointed out an error in a previous derivation of the second term of (2:14). 
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with 
= Z (i-1) B 1 Aw» 0 — Aa» (T)] 
a, (1) =) ison =: (0) 
oye 1 jet) oe (2-14a) 
(2l+1) ! 00) (2/+2)! 


where JE, is given by (2.9). Eq. (2.14) is the low temperature expansion for JF and 
is the second result we have wished to derive. 

If the series S’A(kT) converges, the second term of (2.14) becomes the expansion 
in terms of T. It is easily seen by the rearrangement of the series "97 (k); namely, 


Si ® =SHAET) —A(k+)T)) 


k 


3} A(T). 


k=1 


‘Then, inserting this into (2.13) we obtain 


Yj aay | HOA Ay T* a! (2-15) 
27 1=1 
-with 
a jo Z B, (27—1) 
af == (— {ye AGO). (2- 15a) 


(21)! 


It is worth noting that (2.12) is of the integrated form. In other words, we have 
only to calculate the characteristic function D(w) /D,(w) at any rate, and the difference 
of free energy JF is obtained at once by (2.12). The second term of (2.12), however, 
converges very slowly in general. Hence, it is more convenient to expand it in terms 
of the power series of 1/T*. In most cases it is performed easily. 

Finally, we consider the quantity JC), that is, the difference of the heat capacities 


between the perfect and imperfect lattices. It is obtained at once from JF by the follow- 
ing expression. 


dCp=—27 eT (dP \, (2-16) 


bw, 


From (2.12) and (2.14) JC) is given in the form of series. When JF is given in 


the form of series in terms of 1/T? or T, JCy is, of course, obtained in the same 
form. 


S$ 3. The characteristic functions 


In this section we shall review briefly the method of calculating the characteristic 
functions developed in (I) for the convenience of the later calculation. 
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In n-dimensional lattice, the time-independent part of the equations of motion of 
the atom at the lattice point s= (5, 5,°:-, 5,) ate given by 


ood Mauet SV PE =0 OP yo yey) G-1) 


ayy ee, 


Pes 1 9 Q : G : 
where u,= (u,’, u,7,---,u,”) is the time-independent part of the displacement vector of 
the atom s from its equilibrium position and /'(s, s’; k, l,) is the component. of the 


“force constant tensor ”’ 


-If we use as the characteristic functions, D(w), D,(w), the determinants elimi- 
nating all the u, from the equations of motion of the imperfect and perfect lattices, re- 


spectively, these are given in the following forms : 
D, (w) =det|a(s, k; s/ 1) |, 
D(w) =det[a(s, k; 8’, 1) +4 (s, k; 8’, 1) ], (372) 
bk l==1j 25-0 7S 
( $s’ 1,27" N 
The detefminants D(w) and D,(w) are of the nN”-th order when the total number of 


the lattice points is N”. If the number of the defects is small compared with that of 


the normal atoms, we can expand D(w) in terms of D,(w) and obtain 


D(w) 
D, (w) 


=det[ Oy i:,2 + >14 (s, k > s’’,q) a’ (s’, l > 8”, q) ie (3 : 3) 
(s!1,9) 


The right-hand side of (3-3) is the determinant of p-th order if the maximum number 
of non zero 4 in the row of D(w) is p and the summation of (s’’, q) is taken over 
all the values of s’” and gq for which 4(s,k; s’’,qg) 40. The a-'(s,k;8”,q) is the 
(s, k; s’’,q) element of the inverse matrix of D,(w) and is given by the following ex- 


pression : 
N 
ce (s, k 3 oe l) =) [4 (dy) ig Uy y (4;) us (9;) ? (3 ss 4) 
oj=1 


where 4’(¢,) and ui,(¢,) are the eigenvalue and corresponding eigen function of (3-1). 
We now consider, for simplicity, only the nearest neighbour interaction. Then 


eq. (3-1) becomes 


007 Mosy59-++5, Usyso°-5n 
n 
+2) Fp | Us 897-5 — 1-5 PMs, 50-5 PL — YE 5, 59°05 p°*- 5, |= O, (3*5) 


where 74 is the force constant associated with the displacement of the atom parallel to the 
k-th coordinate axis. Fortunately, each equation of (3-5) has the same form, for different 
values of / and in the equation of motion of the /-th component of u, there appear only 
the /-th’s and not other’s. Then, if we rearrange the rows and columns of D(w) and 
D,(w) appropriately, they become the product of n determinants of N“th order; that is 
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D,(w) = Dy (w) =I deta, (6, 8") 

D(w) = es (w) =II det[a,(s, s’) +4,(s, s’)], (3-6) 
where 


9 5 » 
a,(s, 8’) =a'a 7,0, er 


m ee a 
+e Saya! (Os, g= DL FOsf} sp $4 207, si) IT Osx, Ske 
i=1 k= 
(RAL) 


Therefore, the difference of free energies of imperfect and perfect lattices is given in the 


form 
AF=S\4PI 
ja 
and the result of § 2 holds for each JSF’. Since we have only to discuss any one of 
AF”s we hereafter omit the afhx j for brevity. a 


§4. Free energy of the lattice with one defect 


4/1) n-dimensional monatomic lattice. 


We consider as the defect the following three models: (see eq. (3.1)) 
A) Impurity 
Mo I beled bir 
B) Isotopic impurity (Isotope) 
M,—M,'. 
C). Hole 
yr’ (only one place) when n=1. 
‘Clearly, B) is a special case of A). Since it may be physically improper to regard the 
hole as an impurity which is characterized as the limiting case of M—>0 and 7—y", we 
temporarily take the above medel only in the linear chain. 
Now we consider the case in which only one impurity exists in the perfect monatomic 


lattice of n-dimension. Then (2n+1) equations become different from (3.5). These 
are 


[wo 77M 5y.0-5p 15-5, — Gir oe" +e tte + ats +27,) |- U sysptle--s, 


ea Ti" [Mesyerrsycte Lowney oh Way ogy ste Tein bade, | 
([Ak) 
+(e: Usy--spt2s, + TK * Us, sys, |=; (k= it 2,---n) 
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[0700 77M" 54--- 540-5, — 27x | a Usy- SK Sy, (4 : 7) 
kal 


n 
/ 
+ Sor [U.S t= Hac ae Sede | 103 


The characteristic function D(w)/D,(m) is, therefore, the determinant of the 
‘(2n+1)-th order of which the form is (3.3). This is rather complicated to treat and 
‘so we consider the following two cases: (i) the simple cubic lattice and (ii) the isotopic 
impurity in n-dimensional lattice. 

(i) Simple cubic lattice 

In this case the characteristic function is the determinant of the seventh order given 
in Appendix. Since the total number of lattice points, N*, is very large, the sum in 
(3.4) can be approximately replaced by the integral. (That is, neglecting the order 
‘of O(1/N).) Furthermore, if we use the Born-von Karman condition, eq. (3.4) becomes 


« pl cos [ (5 —5;') 0,]d0, 
ee. (s, s’ 3 iv’) = aN i = = : (4 2) 
) w?a7~pM+ S)27,(1—cos4,) 
t=1 


-Eq. (4.2) is rewritten by the simple transformation’? into the following expression. 
do (a, 8. fw —\ax exp[ —x(>}27,) (1 +20”) \/TTs,—s/(2mx). (43) 
0 Z=1 Z=1 


Hereafter, we consider the isotropic case for simplicity; that is, 7,=/.=-"-7,=7-. Of 
‘course, our result can be extended to the non-isotropic case by means of a slight modi- 
fication. Then (4-3) becomes 


*_P,(s—s'; w’), (4-4) 
27 


a~*(s,8' ; iw’) = —- 


where 
Pi(s—s'; w')= "ax exp [—nx(1+ 2”) Js, — sf Co (4-5) 
v0 l= 


The characteristic function given in Appendix is calculated easily using these 


functions. The result is 


i / it Ss 9 
ri 1 = (A; + fs) ge =O! (e,+ 22— fs) i 
9 (tO 


pAlb —=€ (2nw’*) P, — 06! (7, —4p.— \s ) 
ain ed" { (2n"P,) (o;— —40,— i) +6 (2n0”P, ) Pes i (4-6) 


+60” (0,2) (1-42 rs) | 


where 
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€=1—M’/M, 0’=1-7'/7, Gm 
and 
Pi=3{P, (000 5 w!) —P, (0015 o')}, 
f2=% {P, (001 ; w’) —P;(011; w’)}, c-5) 
Ps=4{P, (001 ; wo”) —P, (002; w')}, 


P,=P,(000; wo’), P,==P;(001; ow’). 
Furthermore, remembering the following properties of the function P(s—s’; w’) ; 


2nw’ P,,(s—s'; ow’) 0 (w’—>0) 
(4-9) 
/ 
2nw"P,, (s—s! 3 w!) 305.41 (v' >) 


P,(000; 0) =0.4990 


P,(0013; 0) =0.1658 


(4-10)* 
P,(002; 0) =0.0778 
P,(0 113; 0) =0.1045, 
we obtain for the free energy difference, JF, owing to one impurity 
1 “ : yee G(kT ; € 0’) 
AF =—=— log G, (€,0”) +—> ik log—_——_*__*_ (4-11) 
eon” pee Sage Ted) 
where 
G,(€, 0) = (1—€) 1 (1— 0") (1— 0.21060") #(1 — 0.18390")*® 
Chrys (4-12) 


D,(ikT) 


At higher temperatures the first term of (4.11) is dominant, so the vibrational free 


energy decreases owing to one impurity if it satisfies the following condition ; 
G,(é, 0’) AL 
or, from (4-7), 


M ad rid ri 
a * (02804 O22 CBee) (OAL kites 1. (4-13) 
/ 


Now, consider the expansion of the second term of (4.11) in terms of 1/7 a 


sufficiently high temperatures. Since from (2.12) and (2.12a) w'=kT in the second 


term of (4.11), «w! becomes large at higher temperatures. Therefore, the integrand of 


The values in (4-10) are in accordance with reference 11). Although these values seem to have 
the errors of some 0.5%, it does not exert any serious influence to our results. 
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(4.2) can be expanded in terms of (27/Mw;° w”) = (2nw")~' and easily integrated 
using the orthogonality conditions of the cosine functions ; namely, 


lc 1 
P,(s=s'; w!) =( me ie »)4s(3—s') OD) 


2nw 


where q; is defined by 


alas) =(4 ) |: Sips $1 (1 cos # Npiicane Nest “(dena 


m=1 


The special values of q, for the small values of ; aad (s—s’) are shown in Table 1. 


The expansion of the second term of (4.11) using (4.14a) is, however, very tedious: 


Table 1. 


qo(0---0)=1, 9, (0---0) =n, go(0---0) =n? + 3n, 


q(0---1) =0, q(01)=—4, g.(0--1)=—n, 


though straightforward, so we consider the case in which € and 0’ are small. Picking 
up only the first order of € and 0’, we get the following expression as the high 
temperature expansion of JF. 


if y 
Ah -—_loo GA eyo’ 
28 og o( ) 


+3 (e-209|¢@) ato) +0 +32) (,) =~] 


(eB) 
(4-15) 

Next, we consider the low temperature expansion of JF. It is obtained at once 

inserting the characteristic function, (4.6), into eqs. (2.14), (2.14a) and (2.9). The 

integral of JE, however, seems to us too hard to give the exact integration. Therefore, 

we consider only a special case in which € and 0’ are small. Then the integrand of 

ME,, i.e. log (D(w) /D,(w)), can be expanded and integrated, because the integration. 
of P(s—s’; w’)” with respect to w’ is of the Gaussian type ; namely, 


\ (2nw"”) {P,,(s—s’ ; w’) "do! 
0 


1/2 
= oY “(= a pi? (s—s’), 
/ n 


(4-16) 


where p,2(s—s’) is a constant number given by the following integral : 
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os Abode 
cE ITs, —s/ (Xp) d%. (4-17) 


ro 
5 —t/26 
pip (s—s’) =) Jat ac aad bat 


Therefore, JE, becomes, for small € and 0 


boy 


dE, = 


i [€-2nw’ P,(000 ; wv’) +60" iP, (000 ; w ') —P,(001 ; w’)} |do’ 
0 


= 204 [ $3 (000) +36 : a te 4) (990) — p73 (001) } }. 
e-={— os pokes slew, ois A (4-18) 


vl 


Then, the low temperature expansion of JF for small € and 0’ becomes as follows: 


=f =(4 ie p% (000) +30 (= = 3 (000) —pi72 (001) } J 


dF 
ie SS T® a/ (4-19) 
27 Z=1 
where the second correction terms are calculated easily from (2.15a), namely, 
22-1 
A Sesame: | 2 ) {¢ (2m) P, (000 ; «) 
(21)! L\a (4-20) 


+60’P,(000 ; w) —60’P,(001 ; o)} | 
(w=0) 


(ii) Lsotopic impurity. 

The difference of free energy JF due to one isotope in the simple cubic lattice is 
obtained at once by putting 0’=0 in the expansion of (i). However, we can discuss 
not only the simple cubic lattice but generally the ndimension case, because the charac- 
teristic function for one isotope is the determinant of the first order 

The characteristic function for one isotope is given by 


D ae % , 
iw") =1—€2nw"P,(0---0; w'). (4-21) 


D, (iw) 
Inserting (4.21) into (2.12) and with (4.9), we obtain the difference of free energy 


‘due to one isotope ; namely, 
dF = + log 1 +2 Sik log prea e, \. 

2p 1—€ Pk 1—€2n(k+1)*T*P,(0:--0; (k+1)T) 
(4-22) 


The high temperature expansion of (4.22) can be performed in the same way as 


was done in (i), while (4.14) is simplified in the present case. It is necessary to 


calculate only ¢;(0:--0), which is given by 
3/2) | 
Se a LE | ee (4-23) 


y 


q)(0-0) = 
a=0 (j—24) 1G) ™ 2" (Hy) ! 
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where the summation with respect to /4,, is taken over all the sets that satisfy the follow- 
ing condition : 


n 
ap! 


vi 


A=: 


_ 


m= 


‘Then, the high temperature expansion of JF is given by 
1 LF; 2 : —1*\/}) 
4F=—| | 1—e ey eS) ,(0-:-0) (2; —-- | 
2 og ( ) 234i )e¢ i( aan 


= Glee 9 Dontlae &C (2) sa (tte (4) oa), 


+e(4 +24) =f} (4-24) 


In the simple cubic lattice, Eq. (4.24) becomes, of course, the same expression as 


(4-15) in which we put 0/=0 (7/=7). In one dimensional case eq. (4-24) becomes 
i B) 4 —1\" 
7) wien | | 1—€ ae —e€>) ss) 2; (< ) ||. 4-25 
g 8 (135) at € (27) ae (4-25) 


It can be easily proved that (4.25) is the same as (4.34) with (4.35b). 


The low temperature expansion also is obtained in the same manner as in (i). 


Inserting (4.21) into (2.14), 


Ne fe" log[1 — €2nw"P,(0---0; w "Ide! ee ens: Sul ay (4-26) 
T l= iat 


The integral in (4.26) is performed at once from (4.16) if € is small. Then, the 


result is 


dp ten. 4 ar (00-. 0) +225) SIT ay! (4-27) 
27 


where @,'’s are calculated from (2.15a) and (4.21) straightforwardly, i.e., 


a= (| 7) ke Be ei de oe ae 


Eq. (4.27) is to be compared with eq. (4.19). 
From (4.22) and (4.28) the vibrational free energy decreases on account of one 


isotope if the following condition is satisfied. 


7 


€<0 or Ts, (4-29) 
M 


Eq. (4.29) correspords to (4.12). When the normal atom is replaced by the heavier 


isotope, the normal frequencies of the lattice vibrations may be decreased. So eq. (4.29) 
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is expected naturally. 


4/2) Monatomic linear chain. 
Monatomic linear chain is the only case in which the integral of the Green function 
(4.2) can be performed exactly. Therefore, we shall investigate this particularly. 


The Green function of Monatomic linear chain is given by 


Is—s! 


a‘(s,8'; iv’) = eel (s—s’; w’) wee oe ; wo >0 (4-30) 
Phe fA Sorat 4 
where 
z=[V 0” +1—o’f. (4-31) 


The characteristic functions of one impurity, one isotope and one hole are calculated 


easily and the results are (See Appendix II.) 


Bi ! 
Ae Rae Nea ; isotope, (4-32a) 
Da se) V1+0” 
D OE as 0 ne 
EAeSLIESS (5B) = ; hole, (4-32b) 
D,(io' ; °) ¥ t-+a- 
se hie Y , 
aol ae! 2 i Ba (24 20) 
D, (io! ; € 0) V 1+o” 
/2 eee 
4+46e—¢__[Y1+o0” —o’] (4-32c) 
V 1+o” 
pee te Aa oe = [1+ —o’f ; impurity, 
Z wo- 
where 
/ , — 0 
€=1—-M'/M, 3=1-7/7r'=——- (4-33) 


1—0O 


In eqs. (4.32b) and (4.32c) we have used the fact that JF is invariant even if the 
constant factor is multiplied to the characteristic function. It seems to be that from 
(4.32) € in the case of isotope, 0 in the case of hole and (€+2n0) in that of im- 
purity are corresponding to each other. This may be due to the properties of our models. 
as defects and, therefore, holds in other cases. 

Inserting (4.32) into (2.12) we get the following expression of JF. 


fm 4og G.(e,.8) to Se ee 
28 : = °G(e. 8: (k+1)T) aM 
where 
(1—e)7" ; isotope, 
G, (4.0) 5 (i=) ; hole, (4: 35a) 


(1g ere) ; impurity, 
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and 
kT f 
i ¢ = / isotope, 
VY 1tRT? 
[eae ee ; hole, 
V 1+kT? 
1 oat 20 ee (4-35b) 
G(¢,0; kT) = Volt 
R eT? a : 
V IteT? ete 
: ad es re eae 
+40? (1 — €) —~—_7p1+ FT? —kT 
(1-€) er + ] 


; impurity. 


From (4.35b) the high temperature expansion of JF is obtained at once. 


il y ee 1 
JaF=——los G,(€, 0) + —— Tt Wes for. | = | 1, 4-36 
‘The first several values of f, are shown in Table 2. The expansion for 1 > (2G,—1) 7 
>0 is also obtained in the same manner. For the impurity, eq. (4.36) is correct so 
far as the first order of € and 0 are concerned. 
Above the ordinary temperature the first term of (4.34) becomes dominant. Then 


the condition that the vibrational free energy decreases owing to one defect is following. 


G;(€,40)=1 
or 
(M’/M) >1 ; isotope, 
Co. ae ee! ; hole, (4:37) 
CMM NTL). 2 ; impurity. 


This corresponding to (4.12) and (4.29). 
The low temperature expansion is also obtained easily by inserting (4.32) and (4.34) 


into (2.15); namely, 


Vi 4E,+>) apt (4-38) 
‘with 
[ bor et fine 17 2G aye bw;,; 2G,—-10 
FA ge chy . (4-392) 
ple Gon lo 1—y/1—2G6, are ee 1—2G,>1, 


ari qoL—2G, 1+//1—2G, 
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Fig. 1. AF due to isotope or hole in the monatomic linear 
chain. 
(1) :--M’= (2/11) M, or 7’ =(11/2)r. (2) ---M’ = (2/3) M, or 
7 = (3/2) r. (3) ---M’= (4/3) M, or 7’ = (3/4) 7. (4) ---M’= 
(20/11) .M, or 7’/= (11/20) 7. (5) ---M’= (20/9) M, or 7’= (9/20) 7. 
(6) ---MW’=4M, or 7’= (1/4) r. 


and the first several values of @,’ are 


PS at oh Oe [3 Cet —2( So? ) ne (4-39b) 


1 > 


nC * 4) Vath 
The expansion (4.38) is that for small € and 0 and, for the impurity, the values of 
eqs. (4.392) and (4.39b) are correct so far as the first order of € and 0 is concerned. 


In Fig. 1 we have plotted IF due to one isotope or hole. The deviations from the 
classical treatment of Kirkwood and Stripp'’ are seen clearly. 


Table 2. 
1 
mol (Gua\F 
f= G-1EQ) 
fa—- (G)4-2) (Ge) 


f= - (G,—1).{15+97G,—1) +21G Soe 6) 


where G, is given by (4.35a). 
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AC; /K 


Vie Tees Hp) 1/3 1 3 5 7 


18 
Fig. 2. Change of heat capacity due to one isotope or hole in 
the monatomic linear chain. 


(1) ---M/= (4/3) M or 7 = (3/4) 7 
(2) ++: MW’ = (2/3) M or 7/= (3/2) 7. 


Finally, the change of heat capacity due to one defect is obtained at once from 
(2.16). We have plotted the schematic graph of JC, in Fig. 2. 


4/3) Diatomic simple lattice. 


It has already been shown in (I) that the equation of motion of the atom in the 
diatomic simple lattice reduces to that in the monatomic simple lattice by the follow- 


ing transformation : 


nh 
PR goes SV y+, P/2 
ed Fs (id 5 pt—2> , aa ! us, 
k=1 


M, when 5} 5; is even 
= a (4741) 
M, when 3S} s, is odd, 


61 


x 


and 
wo wo; M*= a te t+[ (w?02M,—2>57,) (v0? My—2>371) )"" (4-42) 
Furthermore, the optical and acoustical bands are obtained by 
wo o{M* =2.S) 74(1—c08 Q,) - (4-43) 


In the case of linear chain, for instance, eq. (3.5) becomes 
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[wo M,— 27 |ugs +7 tos +1 +7 ¥2s-1= 9 


(4-44) 
[ wv 2M, — 27 \los 41 +7 vo6+2 + Tes = 0 
which, by the transformation 
Ug =[e? wz? M,— 27} I thas 
Uses = [2 2M, — 27 Fue 15 
reduces to the equation 
[ ww 2M* — 27 lu, +7241 +7 U1 = 9 (4-45) 
where 
w? 02M* = 27 +[ (wo 2M, — 27) (wo 2M,—27) }”. (4-46) 


Eq. (4.45) reduces to the equation of motion of monatomic linear chain if M* is replaced 
by M. 

Then the Green function a~'(s, s’) can be calculated in the same way as that of 
monatomic lattice. The expression of z, however, becomes different due to the fact 
that M* is the function of w. Then it can be easily proved that the Green function 
is given by (4.31) if z is replaced by 


= ral [ (o?0 2M, + 27) (w°o7M,+ 27) Li 
~[(o'w2M,-+27) (o%o2M,+ 27) — 47°}. (4-47) 


The characteristic functions for one defect are calculated in the usual way and the 
results are 


(o%e "My +27) (w% aM, +27) 


er 1/2 
a ae es] oe coum] —— | 
2 (w?0,°M,+ 27) (@°«2M,4+ 27) —4° 


9 


(4- 48a) 
; isotope, 
SAR ee ae) fee Me a (4-486) 
Dd, (iw 4 0) | (w?w°M, + 27) (w?w? + 27) _ 47°] a 


x [ (wo 27M, + 27) (w?w 27M, + 27) Te 


; hole, 
where 
M, ; 
€,=1- 2 ; é6=1— 7 7 
2 M af 


Therefore, JF is given by 
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G(€,, 0; kT) 
2 


4F—— log G G, (€, 3) + ee J Parris (4-49) 
with 
Cha ; isotope, 
Cole, 0) =i (1—d)7 ; hole, (4-50a) 
eo) eas) se simpunty, 
and 
€, 
= ie Myo ET 02M, PT? £27] 
G)(e, 0; T= ; isotope, (4-50b) 
La as —_—_| (07° M PT? +7) (o2@M PT? +7) —7"] 
\ ; hole, 
where 
A(kT) =| (oP MRT? + 27) (oP MRT? + 27) —47?]}!” 
[oe MET? +27) (of MRT? + 27) 1, 
pe cali ee dbadlced : (4-50c) 


Ze M, M, 


The generalization to the n-dimensional case is straightforward but tedious. 
Comparing (4-50a) with (4-35a), it may be clear that the sign of JF is controlled 
by the sign of €, which in turn is controlled by the relative magnitude of M, and M,. 


2 
Therefore, the defect whose mass is between M, and M, is apt to be replaced by the 
lighter atom. At low temperature, however, such a simple discussion may not be possible. 


Actually, the calculated expression of JE, has not such a simple form”. 


§ 5. Interaction energy of defect 


We define the characteristic function owing to the correlation of m defects, K,,K,,---, 
K,,, 28 follows : 


Cam COM K,, Kets K,,) 


Ms wild : ml 
EB = 3) GT (o! ) Kup Kian eae] 
eee K;) aE (5-1) 


where 
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Dio! 3 Ky+Ke ++: +Kn)_ (5-2) 
D, (iw! 5 K,+Ke+-+ +Kn) 


Gols KAKA + Kn) = 


is the characteristic function when m defects are present. Summation with respect to 


(js jostttsfe) is taken over all the i sets. Especially, when m=2, 


G(w'; K,+K) 
G(o' ; K,)G(a’ ; K) 


Gi(w'; K, Ky) = (5- 1a) 


From this definition the characteristic function of m defects, (5.2), is expanded as follows. 


Gil Kb 4K) HG (a's K)IL43S2 SG; a eae 
¥) I Jd 
ee t=2(jy-Ju) 


Now, since the interaction energy is defined as 
AF?(K,,» “sho = 4F (K,+K,+--- +K,,) 2) dF (K;), (5 -4) 
ja 


we get from (2.12), (5.1) and (5.3), 


ART (KKK) => AF (Kj, Kj--sKjd (5-5) 

with 
Ar} = 7, [Al(0) AY} $1 Ai ET) —AN(k+1)T)] (5-62) 

and 


é Sas Gi (ow! 3, Ki: X_) 
Al (w’) =log}| 1+— ie |. (5 - 6b) 
lc Selb eaae(@ pRGy oR) 
p=2(jfi--Jp) 

If we consider the isotopic defect, it is easily proved that the first term of (5.6a) 
vanishes. In the elements of the determinant G(K,+K,+-:--+K,,; w’) there exist those 
not included in the determinants G(K;; w’)’s, #=1,2,::-,m. Because those elements 
have the form that is proportional to w a~'(s,s’; iw’) where ss’, they vanish at 
w'=0 and w’=co from (4.9) and so the first term of (5.6a) is zero. This means 
that the interaction energy of defects is as small as can be regarded as the correlation 
to AF of one defect and that vanishes at infinitely high temperatures, which may be ex- 
pected physically. In fact this can also be proved in the case of the impurities of the 


linear chain. Then, it may be concluded that the first term of (5.6a) is zero in general 
and (5.6a) becomes 


1 = fi I 
AF f= i ak [Ai (kT) — A? ((k+1)T)]. (5-7) 
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5/1) Interaction energy between two defects. 


First, we investigate the correlation between two defects. In order to know the 
essential properties of interaction it is enough to consider only the isotopic case, although 
the impurity can be discussed in the same way. 

In the monatomic linear chain the isotope and the hole are equivalent and we can 
treat them simultaneously. If we put 

€ for isotope, 
jee (5-8) 
| 0 for hole, 
the characteristic function due to the correlation of two defects calculated from (5.1a) 
is given by 
12 
wo” ar Alp (Set 1 
ov Ly Te —w' | |sy—so|-FI+ 


1+” 


/ / > 
{i-K 5-2} 1-K5 2 | 
Vy, 1-0"? VY 1--o? 
where the minus sign is taken when the defects K,, K,, are of the same sort and the 
plus of the different sort. Inserting (5.9) into (5.7) we get the interaction energy of 


Gi! (o's KK) = # KK, (5-9) 


two defects, 


AF{=S} logh1+G} (kT ; K, K)]} (5-10) 
k=!1 


where we have used the fact that the absolute magnitude of G,’(K,, K,; w’) is less than 
unity for all the values of K,, K, and w’. Also using this fact it is clear that the 
sign of JF,’ is decided by the sign of G,’(K,, K,; w’), that is, the sign of (#-K,K,). 
In other words it is controlled by the sign of (+K, K,) according to whether the inter- 
action between two defects are attractive or repulsive. This is shown in Table 3. 


The result involves that derived in (I) as a special case. 


Table 3. 
AFL K, K, 
f (=) M/<M(y>r) M, <M (7! a iy) | 
same sort 

| (a4) Mo MGy <7) M,! > M (7! < ) 

attractive ! i 
iC) M/<M ie <7 
| | different. sort 
\ =) M/>M ee Saif 
(> M/<M(y/>7r) My! > M (j-s' <7) | 
| same sort 
iG) M/>M (7! <7) M,' <M (j.' >7) 

repulsive | Sie 4 
i < ie > ij 
Bow ; : different sort 
\ Gs) M/>M fear 
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Although the above discussion holds at higher temperatures, it is also valid at low 
‘temperatures. It is seen at once because the interaction part of JE, is from (2.14) and 


(5.9), 


dei= oa logl1-+G,"(w! s Ky Ke) Wo!’ (5-11) 
DT ado 

‘Clearly the sign of JE, is determined by that of the integrand of (5.11). This integral 

was calculated in (I) in the special case of small K, and K,. That 1s, 


\ 


(2+ 4s+1F0 ) 
4 


boy, | 2(4s+1¥1) 


JE}= (+ K, K,) : 
Ar \16(4s+1-F1iyf—1 


+o 4+ (45+1¥1))}; 5=|5,— 50], 


while, at high temperature (5.10) becomes approximately 


JFI~ pie FK RS (4541 )(, fe? cee 


These results are at once extended to the case of isotopes of n-dimensional mona- 
tomic simple lattice. The correlational characteristic function of two isotopes €,, €s is 


given by 


Gio! ; €,, €,/) = (ee, «) [ 2nw"P qe . w’) F ew 
{1 —2nw" *P. (Oe<-0,§ w’) €,} {1 —2nw’ °P.(0---0; Wa) €} 


-where P,(s—s’; w’) is given in (4.5). Also in this case it can be proved 
Oe, 3) [a1 and, then, the discussion about the sign of JF 7, is retained. 

The high and low temperature expansions can be done in the same way as in 4/1) 
(ii) altogether. The results are : 


dFJc~ + log [1—€,€,'0 (4 , 
v= og | x (48) (5 1) i 
s=Sl5—s/], (5-14) 


3hw 7 9 
dE,'’= —€, €4' 3 y" Q(s—s 
nr TV aGt ( 2n Pore at 3 


° (‘co et +t) n 
pip=const = |" \4x — 
J0 


(xx!) 9? 5 Le Ts;—sy (x) Tsy—sa/ (*') (5-15) 


Eq. (5.15) is the result of small €, and €,’. 


For a pair of isolated isotopes at a separation R= (>‘(s;—s,')’ *)'251, we find the 
following distance dependency of JF, at high temperatures, 


dF *=[P.(s—s' ; T) PooR-“- exp[—Moz*T*R), 
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while at low temperatures according to (I) 
dE? co Ra enc 1) 


Obviously the correlation of a pair of isotopes at high temperatures is more short-ranged 
than that at low temperatures. This result is to be compared with the result of Stripp 
and Kirkwood" that the correlation of a pair of Schottky defects is proportional to R-° 
at high temperatures. 


Finally, we calculated the characteristic function due to the correlation of three or 
four isotopes. The results are 


12 Re ty 
Gia! ; K,, Ke R21 Ini es ot) (5- 16a) 
asiet ian €,2nw”P, (0 a0) 3 w’) 


Gi (o's K; K,, ks K,) 


€,2nw”P,, (s; LTS wv’) 


poke a a ee 5-16b 
(rings) (éjkl) 1— €,2nw"P,(0-+-0;0) 


The summation of the right-hand side of (5.16b) is extended over all the permutations 
of rings. The minus sign is taken when the ring is united, while the plus is taken 
when the ring of four defects are devided into two rings, one of which is composed of 
two defects. From this the interaction energy due to three or four defects can be 
calculated at once. The magnitudes of interactions are O(G,’(w'; K,, K,)*”) and O(G,/ 
(w’ ; K,, K,)*""), respectively. Remembering G,’(w’; K,, K,)'”<1, it is clear that the 


strength of higher order interaction decreases as the order increases. 


5/2) Interaction of a defect with the boundaries. 


We have used the Born-von Karman condition till now. However, if we use other 
boundary conditions, for instance, rigid or free boundary, there appears additional term 
in JF, which we may interprete as the interaction term between the defects and the 
boundaries. 

Understanding the essential features, it is enough to consider only the monatomic 
linear chain. Since the wave functions and the eigenvalues are determined so as to 
satisfy the boundary conditions, Green function a'(s, s’) also depends upon them. Neg- 


lecting O(1/N), we obtain them by simple calculation. (See (I).) 


js—s?] 
rome a EKG 
i See 

; 1 |js—s/| ____|s+s| 
at(s #3 iv’) =) — + fe > RB. (3°17) 
iid Ve WE 
|s—st| \s-+s/—1| 

Vs +z ae ie 


v z—z" 
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Now, we define the interaction energy of a defect with the boundaries as 


AFL »p =4Fpy.—4F 2.x- (5-18) 


Bs De se 


If the defect is an isotope or a hole, from the simple calculation like Appendix II eq. 
(5.18) becomes 

p= Di logli+G4.5 kT; K)] (5-19) 
aoe k=1 F.B. 


with 


1 P+41/2 


1 pe euliai tad ae ee 
Grz=t K 77 cia 2 
FB. V O) 1—K’ 


(5-19a) 


/ 
wo 


V1+0” 
mere s is the distance between the defect and the boundary. Again, since G!(R’; w’)| 
xe 1s 
(4F%,,) 20 when (+K’)=0, 
BR. 


or 


M<M (7'>7)---RB. 
attractive (JF’< 0) “| 
M>M (7 </y)---F.B. 
M'>M (7'<7):--RB. 
repulsive (JF’>0)--: 
M <M (7'>7)-:-F.B. 


The high and low temperature expansions are given by usual calculations ; namely, 


AF, 5 ~i log 1+ K’C (2s— 1 Ps E )( 1 sara (5-20) 
ra. p "earth gitar 2 nek ; 
I es , bw, 1 y 

JER eek ; (K’~0). (5-21) 
FB. Ce et 


Eq. (5.21) has already been obtained in (I) for the case of isotope. 


$6. Free energy of the lattice with irregular 


distributed defects 


We shall interprete the lattice with irregularly distributed defects as such that each 
atom of the perfect lattice is accidentally replaced by a defect with a certain probability. 
Then the characteristic functions corresponding to such lattices are obtained by following 
two processes ; (a) fixing the total number of defects and taking an average with respect 
to all the possible relative configurations of defects and then, (b) taking an average with 
respect to all the possible total numbers of defects. 
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When the total number of defects m is of the order N, the characteristic determinant 
becomes that of the order N, which is, of course, impossible to calculate directly. As 
was noticed at the end of 5/1), however, the higher order correlations may be neglected 
if the correlation of defects is not so strong. Therefore, we have only to calculate the 


first several terms of (5-3) and may put approximately 


(G(o'; K,+K.+--+K,,) = <11G (a! 7K, exp> | [ >)G,@s Kj, Kp) >|. (e-0) 
t (jade) 
where the notation ( ) means the average of type (a). Since the function G(w’ ; K;) 
is independent of the relative configurations of defects it is invariant by the average of 
type (a). 
From (2.12), (5.2) and (6.1) we obtain 


where 


EK, Kee Kn) = SGT Kj,---Kj,) )- (6:3) 
io (Je) 
In (6.3) we have rearranged the series with respect to k. 

Henceforth, we shall treat the monatomic linear chain with the isotopic defects (or 
holes). Our results obtained here can be easily extended to the simple cubic lattice and 
the impurity, although the calculation becomes tedious. Also we shall consider, for sim- 
plicity, all the K,’s to be equal to K. 

Now, we shall take the average of type (a). First, we suppose that the probability, 
with which each normal atom is replaced by a defect, is independent on the configura- 
tion of other defects; namely, it is determined only by the concentration of defects m/N. 
Then, remembering G,’(w'; K;, K,) is given by (5.9) and 


4 —1 N- 
(ot etemaly= iss (2S 4 sree 


(95%) 2 sj=1 As 


tiem lag eee 


we get the average value of (5.9), 


€>).G2(a) Foes Ry) 


(j, &) 
=m(™)(—K’) ow? / awe), z Pel 
Piy-sseri 


aca 
V ct +o” 


2 


Leg (6-4) 


Substituting (6.4) in (6.2) and (6.3), we obtain 
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(AF (mK) ) =m{ JF (K) + © (AEM) \] (6-5) 


where (1/m)(4F’(mK) }) is the mean interaction energy per pair of defects and is given 
by 


(1 AF! (mK) ) 


#T"_/{74EE 


a yo ae ea 
=k \-Oy ai FS 1 ETP? =4r}'| mains: 


Eq. (6.6) is the first term of the expansion of interaction energy with canine to con- 
centration. The next terms are obtained by calculating the terms t=2, 3,---in the right- 
hand side of (6.3). It may be clear that the interaction energy per triangle of defects 
is O(m?/N?) which is small compared with that of pair, when the concentration m/N 
is much smaller than unity ; because the pair’s is of the order O(m N). 

Next, we suppose that the probability, with which a normal atom at a distance s 
from a defect is replaced by a defect, is proportional to exp(—4F(s)éw,/2=T). This 
assumption seems to take account of the correlations between defects approximately. 


Then 


—1N 


(S} discal) =—N(2) = ran POs + Sexe’ fi?) 


mm > ere PARI SAR yy Se DOE. 28 


=> q! (1—z?**) g=0 gq’! 1-7" 


where f is determined by the condition : 


a3 S (SAD aoc fe es 


The summation (6.7) may not be expressed in the closed form, but it is sure that the 
value of (6.7) is between m(m/N)z and m(m/N)z/(1—z) and may nearly equal to 
m(m/N)z. Hence the average value of (5.9) becomes in the present case 


RPT? 

1 m ae /1+kET? 

pivaled 74 | q ~ ‘dl tat 2 “\ l 7 Ohno — 

a be (mk) ~(™ )« KR) > eer [\/1 +RT?—kT}* (6-8) 
1— I =—S>S= 
| V FR 


which has to be replaced for (6.6). Clearly, if the other factors are not considered, 
(6.7) may be more probable than (6.6). 


We shall, then consider to take the average of type (b). If the mean concentration 
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Bal N is much less than unity, we may approximate the binary distribution to the 
Poisson’s. Then, : 


JINGlOn en ON ens) m, 
RGCi in. KS pris 0)" [6 (0's K) Peap [3X3 Gis 6))) 
‘ TOE fr 


~~ / 1 el 
= exp|m,G (w ; K) exp {2 <2) Gs (ow) ee ea’ my | (6 : 9) 
0 & 


where we have replaced (1)//m)(>1G,'(w')) by (1/m) (S1G.! (w’) 5 


contributes to the avarage where mm, and also because exp (m,~'(>1G’,)) is very nearly 


because it mainly 


m=mo 


equal to unity. The notation “bar” means the average of type (b). Use of. (6.8) 
and (2.12) leads at once to 


(AF (m, K) ) =m) dF (K) + 4F’(m,K) >| (6-10) 
with 
ee) = kT k+1)F 
OF CK) = 1.9 aS eM M1) 
22) uP a )ky Rao ae Saal Cee" 
and 


0 


CAF? (m,K)) = pat (kT ; K) fexp (CSG ET +-m,K))) —1} (6-12) 


—G(kT+T ; K) {exp ((3)* GIRT+T ; m, K))) —1}], 
mM) 
where G(w'; K) is given by (4.32) and (1/m)(>'G’,(w'; m, K)) is obtained by put- 


ting m=m, in (6.4) or (6.8). 
The high and low temperature expansions are also easily calculated. The results are 


CIF (m, RY) =F mK S77 @D(T) |+CP Gu) (6-13) 


and 


(HE, (m, K)) = ner mK CAE aK) (6-14) 


where we have omitted the detailed expressions of (IF !(m)K) ) and ( dE,’ (m,K) > because 


they are of rather complicated forms. 
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Appendix I 


The characteristic function of the simple cubic lattice due to one impurity 1s the 


seventh order determinants given by 


Diiw’) _ 
Dy (io!) 
1+, 4, =p, 4, — pa’ 4, — py A, — pz 4, — ie ds —ps' 4, 
= 4s 1 + ,'4, — fre! dy — ie! Ay — pA, — fis! dy — fie’ de 
— 02 4, Pe 4; 1 +p 4, — pz 4, — fs 4d, —p,' 4, — py! A, 
A, A, wh 1+/1, A, A, A, (A-1) 
— pe 4; — 2! 4, — pa! 4, —p'4, 1+ 4, — 2d, — p74, 
— pe’ dy — 3! , — fs’ A, — fro! dy —fie'd, 1+’, — po’ de 
— ps! 4, — 0,1 4; — fie’ 4, — (ry A, — fi’ A, —p'4, 1+ 9,'4; 
with 
=a '(0 0 0)—a*(0 1 1) 
(2 =a'(0 0 1)—a "(0 1 1) 
p3 =a '(0 0 1)—a "(0 0 2) 
where 
da" {0 0 1) =a~"{s, &,) “with s;—s,=s(G. 0 2) 
etc., 
and 
4.=7l =7(1 7) (i=1, 2, 3) 


: 
A,= Mo” or eae (Ss, $s) +2(4,4-4,) p.' +24, (03 — 9") 
(j, &, =1, 2, 33  cpedic) 
A’ =Mo” wy € a= ss <<) +2(d, + 4,+ d,) (o,'— 4’). 


If we consider the isotropic case 7,=72=7,;=7, ~;’s (i=1, 2, 3) reduce to —7f4’s given 
by (4.8) and all the 4,’s and A,’s become the same respectively. Hence (A.1) can 
be readily calculated and the result is (4.6). 


Appendix II 


In the monatomic linear chain the equations of motion (4.1) become 


Effects of Defects on Lattice Vibrations 355 


Oop MATa 1) lustnshy tea) G0 
[w?o?>M — 27" lu, +7! Ug .4 +7/us_1=0 (A: 2) 
[ww2M— (7 +7’) Jussi tru t7/u,=0. 
Hence 4(s, s’)’s of (3.6) are given by 

4(s—1, s—1) =4(s+1, s+1) =70=4, 

A(s—1, 5s) =4(s, s—1) =d(s, s+1) =4(s4+1, s) =—4, 

4(s, s) =Mw"o? € a*(s, 5) +270=4,4+24, ue 

4(s, s’)=0 otherwise, 


and the characteristic determinants are 


D(iw’) = 
Dy (io") 


ta Ae OE Mea WSS Gao Dak 4 Ge Derbi d) 
bexp—=2 k=s—2 k=s—2 | 
| 


e+] s+] s+1 
Ds A(s ka" (s—1, &), 1 +246, k)a~" (5, k), Si4(s, kya" (s-+1, k) 


k=s— k | 
| 
s+2 s+2 $+2 | 
pas ik 1, k)a~* (s—1, hk), 14 (s+1, ka (5,k), 14+>)4(s+1, ka '(s+1, &) | 
i k=s k=s 
(A: 4) 
Inserting (A.3) and (4.31) into (A.4) and after some calculations we get the following 
result : 
Dio’) w! 9 i ileal ily eesieaaes 
= =] — (€ 4+ 26) - 4+ 406 Yh? =! 
D, (iw) \ 52) +0” V 1+” M4 Smet 
30 vw” (ae : 
+ 40 (1 ea 2 a a ly 1t+w?— w' | 3 (A: 5) 


This is the same as (4.32c). The characteristic functions due to one isotope is obtained 


at once by putting 0=0, (7’/=7), in (A.5), and the one due to one hole is also 
Yin. g 


obtained in a similar way. 
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Regarding the pion as quantum of the lowest exciting mode of nucleon-antinucleon pair from the 
vacuum, we get a normalized wave function for the pion and use this wave function as an auxiliary 
variable for the treatment of the many body problem. Then we perform the canonical transformation 
and obtain an effective pion-nucleon interaction. The result is equivalent to the customary meson 


theory. 


$1. Introduction 


Various phenomena related with the low-energy p-wave pion-nucleon interaction were 
explained successfully by Chew and Low” by using the cut-off Yukawa theory without 
nucleon ‘recoil. But, problems which seem to have a close connection with the s-wave 
pion-nucleon interaction, e. g., s-wave pion-nucleon scattering, anomalous magnetic moment 
of the nucleon, the radius of the nucleon core, the electron-neutron interaction, etc., are 
still unsolved on the basis of a satisfactory theory. Probably, the virtual nucleon-anti- 
nucleon pair plays an important part in these problems. If an adequate method to treat 
the nucleon pair is known, these problems will be solved successfully. 

On the other hand, in the composite model of the pion proposed by Fermi and 
Yang,” the pion is regarded as a bound state of a nucleon and an antinucleon and, 
therefore, nucleon pairs play the essential role.* Already Machida” attempted to treat 
the pion-nucleon scattering on this model. We also here regard the pion as quantum of 
the lowest exciting mode of nucleons and antinucleons from the vacuum. In section 2 
we consider the eigenvalue problem for the bound state of nucleon pairs using the new 
Tamm-Dancoff method, and identify this state with the pion. Then we obtain a normalized 
wave function for this state and use this wave function as an auxiliary variable for the 
treatment of the many body problem as in the case of plasma oscillation.” In section 
3 we perform the canonical transformation and obtain an effective pion-nucleon interaction. 

We proceed along the treatment of the plasma oscillation formally, but, of course, 
there is a sharp difference physically. In the plasma oscillation repulsive, long range 
force plays the essential role, while strong, attractive, short range forces are concerned 
with our problem. Nevertheless the point that makes our treatment possible is that 
interactions between particle-antiparticle pairs play the main role. 


* A composite model of this sort was considered for strange particles by many authors.®) 
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§2. Eigenvalue problem for bound state of nucleon-antinucleon pair 


We assume that the interaction Hamiltonian between nucleons and antinucleons has 


the following Fermi type, 


il i ae Bs 
H! = 2S) | GO rahe) I Wrst! ©) 
4 M? =1 ¥ 
XP O)r5 7h (Fr) —P W) rt! (r)) dr, (1) 
where 
¢ = br 
p=( a and (i =17, C= ; ) (2) 
Yn — 5 
are the nucleon and antinucleon field opetators respectively. <,’s are the usual isotopic 


spin operators. We assume the above interaction Hamiltonian because it is the simplest 
one which produces a bound state of the nucleon-antinucleon pair corresponding to the 
pseudoscalar pion. Of course this point interaction makes our treatment impossible as 
it stands, therefore, following the usual treatment of the pair theory,” we shall resort 
to a cut-off in momentum space chosen of the order of magnitude of the nucleon mass. 
We shall not write this cut-off function explicitly, but use this when necessary. The 
binding energy of the pion is about 2M (M is the nucleon mass), so that the relative 
momenta of the nucleon and the antinucleon in the bound state are expected to extend 
to 2M. However, higher cut-off momentum makes our mathematical treatment more 
difficult. Hence we consider here a fictitious model and cut off the relative momentum 
at the nucleon mass. 

Now we represent the interaction Hamiltonian in terms of the creation and annihi- 
lation operators and pick up the terms which contribute to the nucleon-antinucleon 
interaction predominantly, neglecting higher powers of p/M. 


Then we have 


1 £ a) r * . ak 5 “ 
ee eS (2) t 02 ~) : ~2 \ 
ny aoe \ {= 4p’ +k//2 Bp, oo 8— (p/—k//2) + (p+! /2) Co 02 Mp/—k’/2} 


3 i ps" 0: i p. Yi 
A apie k’/2 Bos ps bo (p"+k//2) + b- (p!’—k’/2) Cos 0, ap! + k//2} dp’ dp" dk’, 
(3) 


where 


Of 27 t4) Oy (4) 


and a,* and 6,* are creation operators and a, and 6, annihilation operators for the 


nucleon and antinucleon respectively. 
Next we consider the state which has one more nucleon-antinucleon pair than the 


true vacuum and treat this state by the use of the new Tamm-Dancoff method. For 
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Hamiltoni d — B! 

this purpose we calculate the commutator between miltonian an ae + i ns ink 
48 . ‘ 

iM Ss aes t the excited states which have two 

b= (p—k/2) 9f 5 (p+k/2) Sp, 2 4p—k/2- If we neglec 

or more pairs compared with the vacuum, and neglect nucleon or antinucleon exchange 

terms between two pairs because it has smaller spin and isospin statistical weight than 


the non-exchange term, then we have the following commutators, 


p2* J 
a, uo Bip Oe Gmee 


ee: a ce ap’ es eB eee —kj2) + b= (p’+k/2) or Os ap’ —k/2}> 
(5) 
ee (p+k/2) CG 0s ay —k/2> FY") 
a | dp! {a BB got in tara ee 
(aia MM. ip’ +k/2” ss (p’ —k/2) (p’+k/2) “ps 0; 4p’ —k/2 
(6) 


In order to obtain the eigenvalue equation for the one pair excited state, we set up 

; : Oi Ble ight d pf ~J 
equations of motion for the operators —dp4k/2B,,5,6—(p—k2) 22d 6— (pik 2) Cp, ps 
dp —b/2 in the Schroedinger representation. If we denote an eigenstate of the total 
Hamiltonian H(=H,+H’) corresponding to eigenvalue E by ¥ and the vacuum state 


which has the lowest eigenvalue E, by %, we have the following equations with the 


aid of the commutation relations (5) and (6), 
(F, (H—E) (—451 4/2 Bp, ibm (p—bjay) Fo) 


= (Ep+kj2t Ep—kj2—E+ E,) OF, — 45+ k/2 By, ssbee (p—k/2) Y) 


= = x 


( {* * . * 9 

: hop Ps" } Pa” Ps +] p. > 
Gin) ppg ts, dp'(L, (—4p’+k/2 Bos p,o— (p’—k/2) += (p’+k/2) Cosa; 4p’ —k/2) Po) 
AF (7) 


(¥, (i— Bye a kin) Grape 4p — kj2 Po) 


= (—Ep + kj2—Ep-—kj2—E+E,) (g, (p+k/2) we ah k/2 T.) 
sey \a ROD (mga nee ee 6° Cc) at y 
(27)3 Pore a L 4p’ + k/2 0s 04 (p’ —k/2) 0 (p’ +k/2) “3 0, 4p’ — k, 2) Poy 


=0. 


(8) 


From these formulae we get the eigenvalue equation 
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{1- oN af E ees —k/2— Ok Ep+k2tEp—kj2t+oK 


x i px” J oo 1 / 2 
= \ dp (P, (4p + k/2 Bp, p,— (p/—h/2) +O" (p/ + /2) Coy os Ap’ — e/a) Yo) =0, 
(9) 


that is 
= eho mee), a 
Be M E aaron —kj2—%; = Ep+kjt+Ep—Kj2+or 
where Ehren 
This eigenvalue equation has two sorts of solutions. One is 
Op Ep +kj2t+Ep—kja> (11) 


which corresponds to a free nucleon-antinucleon pair and the other is 
OO, 
OR~V E+E, (12) 


which corresponds to a bound state. The latter is obtained by integrating (10) and 
solving the resulting equation with respect to w,. In doing so we neglect terms higher 
than quadratic with respect to k/M in square root of (12). If we cut the integral in 
(10) at the nucleon mass and take 7=6.6, we can adjust #2 to the pion mass. 

We regard this nucleon-antinucleon pair in the bound state as 7-meson, and replace 
the fermion pair by a boson field effectively. We seek the explicit form of the wave 
function representing this bound state. This may be put in the following form, 


i(k) =A, ©, (13) 
A =P ao pt bi Bo, ppb” (pki) +B B” (p+k/2) Sp, 4p —k/2)4P, (14) 


where ap and f, are chosen so that v! (k) becomes a properly normalized eigenfunction 


of the oa Hamiltonian. That is, using the equations 
(Hy +H’) 2(k) = (E,+o,) 8 (k), (15) 
(H, +H") Yo=E, % (16) 


and the commutation relations (5) and (6), we can see that a, and 3, are written as. 


gas INE 
PSE RES 3 oC) 

+k/2 —k/2 k 

pth] p—k| (17) 
icon ok Nk 


Po. Ep+kj2+Ep—k2t ’ 


where Nz is determined from the normalization condition. Namely we can proceed as: 


follows. (15) can also be written as an operator equation 
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(+H), 4.]=onAb- (18) 
Taking the Hermite conjugate of the above equation, 
(H+ H"), Al=—on Ae (19) 
From this we have 
(H, +H’) 4, Vy = (E,—0K) Ae Po- (20) 
Since Y, is the state of the lowest energy, this equation can be satisfied only if 
A, €,=0. (21) 
Using (13), (14) and this condition, we can fix the normalization as 
Crick), Vik) =(%, dp Ay ) 
=(%, [Ais Abe] fe) 
=40,,0(k—K’) | (\ap|?—|8pl*) 4p 
=0,,0(k—K’), 


where we make the same approximation that leads to commutation relations (5) and 


(6). Thus the normalization condition is 


4 | dp (ap |— |p|?) =1. (22) 


If we determine Nk from (22) and (17) using the cut-off procedure at the nucleon 


mass, (14) and its Hermite conjugate can be written approximately as 


F 1" i os" Oy +t 2 
A eee fu SS d 7 Fo +k/2 Be gt Cpt “p Oe (sand 2) Cos “p—k/2 23) 
Jina of aan E +E a E E pos 
V 204 0; 02. p+k/2t&p—k/2—k p+kj2t&p—kj2+k 
po” t 09" oy +i pe 
1 <P fea ee 2 oF Ae ers a 
Ans > {4 p—k/2 "0, 02° —(p+k/2) | (p—k/2) ~p, p2 Bie), (24) 


V 20K pipe En+ k/2t+Ep—k/2+oK En+k/2t+Ep—ki2—Ok 


§ 3. Effective Hamiltonian 


Making use of the pion wave function obtained in the preceding section, we seek 
for an effective Hamiltonian for the pion-nucleon system. To this end, on the total 
Hamiltonian 


‘ * * 
H=\Ep(qp 4p +bp bp) dp 
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1 CTAN 
@n)5 = ove | ¢r dp, dp.dp.dp, 


Oni 02 —i(Pi—po)r___p,* o2* ei (Pit P2)r 
x [4 01 02 Fp © as os Bh py Ops & ; r 


os o eH ay, e\PitP2)r SD 1 fe eq Per) 


pr p24 Do, 5p é 
x [ap pairs Ne ae ey Bi anit —i(ps+pa)r 
+ 5p, Ch actor bp, D ts fi eid (25) 
where 
Beane #2) (9. — ean = = 7.)D,D, 
SN Mt rege Bap) Der 


Dy=V (Ep+M) /2E,, 


and the subsidiary condition 
Ay U (5=0, (21’) 


we perform the following canonical transformation 


OES SO HlP9700 [y- : [iS, [15,0] +-~, (26) 
where 
a? y i p02" 
Ee es) ee e Sif dkdp— =| + k/2 Fem (p—hi2) OO (p+k/2) Ci ts a Ey 
(27) Ms 20% Ep+kj2t+Ep—kj2—% Ep +k[z+Ep—k/2t+oR 

pi" f = pen” 

Sp —k/2 Bo 9,°— (p+ k/2) b” (p—k/2) Soros vethrni2) 4f') (27) 

En+k/2+Ep—kjat+K Ep+k/2t+E, —k/2—k 


We first consider the transformed form of (21’). Making the same approximation that 


leads to commutation relations (5) and (6), we obtain 


01 +1 p2 
bm (p—k/2) ON spn) 


| e+ = ask ja i ay nO” penn 
(27) 3? / 20, E ptk/2teEp—Kj2teK  Ept+k/2+2p—k/2—-on 


1 : 
s\(2 ee 5 Ai gi = 0, 
a aes VK, \a Pe Srreey —k/2 +R)” Tree Nees eon s v 
(28) 
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In (28), we adjust the value of f so that the term containing Aj is eliminated, i. e. 


1 
eee fe bin sialon or 
ey OK (E erates 5 k2toR)? (Ep+kj2+Ep—kj2— bk) 


From this equation f is fixed as f~14.6 with the use of cut-off as before. Then (28) 
becomes 
1 a * eee psn , pK Cor Sp+ke) gr 
0. (29) 
Lie ork +kj2tEp—k2tok Ep+k/ ot Ep—k/2—k / 


Next we transform the total Hamiltonian by (26). Before we do this, we rewrite H’ 


as follows, 


1 gf moons 
H~ (anye err | dp,dp,dp, 4p.) (1]+ I] + (TI) 
7 ip v 


+2((IV]+[V]+[V0+[ViI)}, 


whete 
[ 1 Jaap apt dp dp, Dp, pzDpsp,9 (Pi— Ps + Ps— Po 
i a a ay Kaien in 
[uur ]=oe ies dP Ch os Coup, 9 (Pi + Ps + Ps + Ps) 


fe ee oY Bi 8 . 
Liv i= 4p, 4p; Op b a Dp, ps B 030, ? (Pi P2 + ps+ Pr) 
[Vv j=a a) "pe a “ a’, Ty ci 5 
Py [Ps 0 0: ~ps.04? (Pi — P2— Ps— Ps) 


r — 0 ps td a”: 

[VI ]=—2p, bp, 69, 4p-Dp, o:Dosex® (Pi —Ps+ Ps— Po) 

: = a o* ps ae 2 

[Viljs— bp. » Ops 4p, Pee Cos p, 9 (Pi + P2— Ps— Ps) - (30) 
Here we write only terms relevant to the following calculations. We transform (30) 
by (26) and pick up terms which correspond to the free meson Hamiltonian (we denote 


this by Hy), 
I ps : : 
=a li S, [i S, Ay) 


1 


1 
4(Epy katt -kbia Se 
ke P ; | (Ep+ki2t+Ep—kj2—R)* 


1 or 
+a aw Ai A 
(Ep+kiz+Ep—kj2tor)® k ks 
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+ bs, [iS, 


of ee on "| ¢p.dp. dp, dp, ((11]-+[11]) ]] 


oe 
(27 


< ny reas 
S| dkdpap! 
i P Ey 20; M? 


1 1 
x ; 
Eptkj2tEp_katon Ep +hj2t+Ep—Kj2—oK 


1 1 rye 
a | 4i Aj 
Ep+kj2+Ep—kj2—OK Ep/+kj2t+Ep/—kj2ton)) * 


aa BS eS: oe = \_s | ep.dp.dp, dp, 2[VIt]]] 


>— pee J S| dkdpap! “ oS 


(22)* 20}, Me 
x | : Wie ey seo 
En+k2+Ep—k/2—k Ep’ +k/2t+ Ep —k/2— Ok 
J 
5 


Al. Aj. 
pt+kj2+ Ep—kj2+ or area eee pals 


Making the approximation E,,~M and using the values, f=14.6 and gy=6.6 and cutting 


off integrals with respect to p and p’ at the nucleon mass, Hy, becomes 


Hy~>}| dhe, Ay, Aj. (31) 
j 


Then we pick up meson-nucleon interaction terms corresponding to the p-wave 


interaction and s-wave interaction separately. First the p-wave interaction terms are 


[iS, Aaa =a e | dp. dp.dp.dp,2 1V]+[V)] 

= a a (20 ltd Lee 

A Geecre saaraer {4p 5) 2. ( 2 pee re En): dp, dp; 0 (p;— p,—k) 
44a, f, Ein Ee) _dp,dp.0 (p,— p.—k)| Af 


1 DP Oy ier (p+k/2) FP» Ee “yy 
En+kj2teEp—kj2— Ok {ap, 7 wa En+kj2t+M ES Pa)? Pe P34P2 (P3— P» ) 


: oP — oP _dp,d; i) } Ak 
10 Ce Gaston, E,, ae ip, dp, 0 (Pp, — p.+ k) 
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* 2 (_o:p__ 9: (P+k/2)) 4 ap,.dp,3( k) 
da wT. P Ps? P aed if ae 
aL En+k/2t+Ep—k/2— 1 1G es a Ep+kj2+M a 8 oo 


oP, OP» ,—k) } Al 
+44, ate mine, Pay) oP: dp28 (Ps— Pa) 


1 (* _( op: 9: (p—k/2)), gy dp,o(p,—p,+k 
a ee ed Ee — / P; Ps? (Py PsA ) 
En+kjzt+Ep—kj2t+ Ok | Pi 7 E,,+M Ep—Kj2+M ) Ps 


+4a,, 1G: ca aman Pe) 4p, 4P: P29 (Pi— P+) | aj, |. 


Making the approximation E,~M, the above expression becomes 


Fin~ 


a = S| épap. 


L “Gn Vv 20, 


Xa, =,(3-k) {— Aj 0 (p,—p.—k) + 4), 3 (p,—p. +k) } 4, (32) 
where 
2 (2/3m*-p/M Foy" 
47 4n 


~ 0.08. 


Here we denote the p-wave interaction Hamiltonian by 


H,,. and use the values, f= 14.6 
and g=6.6, and the cut-off procedure as before. 


Secondly, we pick up s-wave interaction terms, 


[iS, Hy] 


=f - S| dkdp (Ep +kj2t+Ep-— 2) 


0" J p2* 
x F _ Sp +k/2%o, o2?— (p—k/2) is “e (p+k/2) C} 


Ai o- k2\ 43 
Ep+kj2t+Ep—kj2—-%k 


Epik ee kj2 TORK 


a Bi ppt bf 
+{- p— 2 Bp, 9,9— (p+k/2)_, & (p—/2) Syn netp+k2| 4p), (33a) 
Eps kj2tEp—Kj2ter Ep+katEp—Kj2—ek 


[is 5 eae 3 ih 


fg os 
- dp, dp,dkd, 
Kame) eh 7 in | es PY: = 
x| tp ao bp, id aN Be 
ed 2 12 3 Al r) et 1 Pies Ps 
es ae ki2 +R (p, + p—- y+ 


En+kj2+Ep—kj2—k 


x Ay, O(p,+ pith) |, 


(33b) 
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pomeafice 
is, ae | ¢p.dp. ne 
ip 


3 SITS, | dpdpdhdp— 
eg 


(27) 9? M? V iy 
62 CI | g@ bP1 Ci gi 
| Pi “01 02° P2 ip Py 01 02 P2 
x = —=Aj O( pet pst ky + - ; 
VED + Kj2+Ep—kj2—K vee En+k/2+Ep—kj2+%K 
x Af, (p+ p.— Ik) |, (33c) 
E esta ay \ dp, dp,dp,dp,: 2 (vit | 
Tl ip 
ad sx? 2! jo uip chap 
(27) 9? *- M? | Pere Py ode 
r : pi 2" p J 2 
x | ans Bi 02 ps 0(p,+p.—k) — bn oF 02 “p2 
En+kj2+Ep—kj2—©k - Ep+kj2t+Ep—k2+%R 
x8 (p,+p.+h) | Af 
A” Bs pee Or Cl gle 
Py- 0102 Po rn) +p,+k) Pi__ i102 P2 
e Ep+k2t+Ep—kj2+k Bins Ep+kj2t+Ep—kj2—%R 
9 (p+p.—h)| Ay |, (33d) 


La[iS, (8, HJ} 
Z 


ae {S| dedp dk’ ap! ot em (Eh pe Be as) 
ij "UK CRW 


1 1 
eggs eo E p+kj2tEp— k/2—%k 


k k/ \ 
Xap 4424 pw 24y AL (p— Pp eae) 


1 . 1 ~ = 
ze Ep 4h j2t ep —h//2t On Ep+kj2t+Ep—k2+0k 


* * fie , k k’ : | 34 
X apm iFtp— hin Ae AL9( PP +5) Gt 


+i ‘Ys [ss, ie = | dp.dp.dp, dp. | | 


0 
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i 2 
> (ony? WD 


(27)? / 20% V 20x 


i! 1 
a ee i 
epee —K2tOn Ep+wjat Ep’ _h/2— OW 
9 ; k’ 
x {ap £55; 4p + W/29(Pi—P -k+—) 
+a, GS z ap— n/29( Pi— Poe * +i) 


1 1 
+ = ‘= BS, Ghee = —— = 
En+k2+Ep—Kj2—@k Ep te j2t Ep —Kj2t OK 


k’ 
x {a 5); 4p'4 w/29( Pi— p’—k+ -) 


+4; Po fee 1/29(Pi— ee E +w)h|, 


Ufc g 

Atay . \4p.dp.dp, dp, (10 | 

; is, [, oe Gay | 4p. 4p.dp,dp,( 011] | 
ie si | dled dk’ dp’ d a 5A, Aj 
eae M° Be P Ta eee 20, 1 x 


1 1 
En+k2t+Ep—k/2+ Ok Ep +k j2t+Epy— Kw 2—-OKn 


* f j I’ \ 

x {ay y —W/27i7)4p.9( Pa P’ +k) 
pea 

+a ay + k/2 7) T4p,9( Pa Pk + + )t 


I 1 
En+kj2tep—Kj2—OK Ep ae j2tEp—ejatow 


: | re: 


* r k ie 
a Gy thin 7 3 Pe (pa p—k +~)} 1, 


+|is,[is, Cae | apap. dp.dp,2{V11}| | 


il 
m)* V ‘20% Vv = 


ie Z| ; g 
>—4+__3}—*_ | dkdp' dk’ 
(any? M p' dk Bers 


| dkeap, di! dp! dq— 2, 5 Ai, Aj 


Ay, Aj 


(34b) 


(34c) 
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1 1 
Epikj2t+eEp _K2+oR Ep wptEp Kt oy 
Pf \ 
x {Bay —h//2° 55; an, 3(p.— Pp’ — ah) 


ae 
ee rae T,) dp, 0) (p:—P;—k—F’) f 


— 1 : i < 
En+k/2+Ep—kj2—K Ep! +h /2t+ Ep —K/)2—OW 
: cabo 
x 8a, 5)Fi4p/ + /29(Pi— P gam) 
* s ; 
—4,, (4—7,7;) dp, 0 (p,— ps tk —k) |]. (34d) 


Making the approximation, E,~M and using the value 7=6.6 and the cut-off procedure, 
(33a), (33b), (33c) and (33d) become 


Sei fone J mpi J as0 jay 02 pe 
(any? ae D| thee LB 2 Boo = (p—k/2) 4k 2" (p—K/2) So, 00% +h/2 4k) 
a jp on™ jhe 8 j 2 
— 7 (a5 1/2 B92 (p+ki2) Ak 9 (w+ki2) Com p—ki2 4b 1+ 

x~y~ —0.05 


and this is negligible compared with the others. This fact shows that the above trans- 
formation separate the 7-meson mode in good approximation. (34a), (34b), (34c) and 


(34d) become in the same approximation 


te 1 py aS EN Wen hee 
Gye 523 dtd eas SET ay tp Ayy Ajo (p’ +k’ —p—k). 
(35) 


Here we denote the effective s-wave interaction by H,. 

In this way, if we choose suitable parameters, we can obtain the free meson 
Hamiltonian, effective p-wave interaction and s-wave core interaction. (31), (32) and 
(35) show that the effective pion-nucleon interaction obtained from our interaction 
Hamiltonian (1) is equivalent to the customary pseudoscalar meson theory with pseudo- 
scalar coupling. Although we obtain our results on the basis of one pait approximation, 
the above conclusion will hold even if we take account of higher configurations, because 
in the new Tamm-Dancoff method the lowest approximation includes at least partially 
the effects of infinite number of virtual nucleon-antinucleon pairs and the cut-off procedure 


makes contributions from terms neglected here rather small. 
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In this paper, several methods of analysing vibrational properties of crystal lattices with defects: 
are developed. Integral expressions for additive functions of normal mode frequencies are derived 
following the work of Montroll and his collaborators. It is shown that the Helmholtz free energy 
can be evaluated at high and low temperatures without performing the integrations. The methods 
presented are valid for lattices of all odd dimensions, although specific results are presented here for one- 
dimensional monatomic and diatomic lattices. Using a method similar to that developed by Lifshitz, 
it is shown that the properties of a lattice with defects can be expanded in a series of powers of the 
concentration of defects. The coefficient of the nt" power depends on the properties of a lattice with 
n defects. Examples of such expansions are given. An exact expression for the frequency distribution 


function of a monatomic linear chain with an isotope defect is given. 


§ Ll. Introduction 


In the study of the effect of defects on the thermodynamic properties of crystals and. 
order-disorder phenomena it has been customary to ignore the vibrational contribution to 
the free energy and to study static effects only. This reduces the problem to a purely 
configurational one. Many methods of solution for these configurational problems have 
been devised. However, it is known both experimentally and theoretically that the 
vibrational contribution of defects and disorder to thermodynamic functions of crystals is. 
not negligible. We therefore present in this paper a number of methods which have been 
found useful in studying defect problems. We illustrate the use of these techniques by 
calculating several quantities of interest for the one-dimesional lattice. Although there is 
very little physical significance to such quantities as the specific heat or free energy of a 
one-dimensional lattice, the mathematical tools that we use in the following calculations 


can be applied to three-dimensional problems as well. The one-dimensional calculations. 


* This research was supported by the United States Air Force through the Air Force Office of Scienti- 
fic Research, Air Research and Development Command under Contracts No. AF18(600) 1015 and AF 18- 


(600) 1315. 
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have the virtue of mathematical simplicity yet serve to illustrate these methcds whose ex- 


position is the primary purpose of the present paper. Calculations for the more realistic 
three-dimensional cases will be described in future papers. 

The earliest calculation of the effect of large numbers of defects on the vibrational 
properties of crystals is due to Stripp and Kirkwood”. They determined the effect of 
lattice vacancies on the vibrational partition function to second order in perturbation 
theory. They showed that the vibrational contribution effectively lowers the energy for 
vacancy formation so that the equilibrium vacancy concentration should exceed that of a 
static model. They also showed that this effect is more pronounced at higher tempera- 
tures. Furthermore the inclusion of the vibrational contribution approximately halves the 
“melting temperature ” calculated according to the hole theory. It is difficult to estimate 
the quatitative accuracy of their results since the use of perturbation theory neglects the 
contributions from any impurity frequencies associated with the defects. These may be 
of the same order of magnitude as arise from the in-band frequency shifts. The present 
work is an extension of methods first introduced in defect calculations by Montroll and 
his collaborators?”°*”™, and takes into account both in-band frequency shifts snd impurity 


frequencies. 


§ 2. Mathematical preliminaries 


In this paper we will be concerned with the evaluation of various additive functions 
cof the normal mode frequencies of a one-dimensional lattice. These functions can be 


expressed as 


S=S! f (5) (2-1) 


‘where the normal mode frequencies @,, are the roots of a determinantal equation which 


we write as 


\M(w) |=0. (2-2) 


‘The zero-point energy of the lattice is given by S when f(z) =$b6z; the Fourier transform 
-of the frequency spectrum corresponds to f(z) =N7' exp (i @ 2), where N is the total 


number of degrees of freedom in the lattice ; and the Helmholtz free energy corresponds 
to f(z) =hT In {2 sinh (6z/2kT)}. 

Using a well-known result from the theory of contour imtegration we can express S 
in an alternative way as 


$= 1 \ f@din \M(z) | (2-4) 


‘where C is any closed, counterclockwise contour which encloses all the zeroes of M(z) but 


none of the poles of f(z). 


When the lattice contains defects we can wite M(w), the matrix of the secular 
«determinant, as the sum of two square matrices 


M(w) =M,(w) +0M(o) (2-5) 
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where M,(w) is the matrix of the secular determinant for the lattice containing no 
defects, and 0M(w) represents the perturbation due to the addition of the defects. We 
can now express M(w) in a form more convenient for our purposes by writing 
M (a) =M,(@)|I+M,"'(@) 0M (a) | (2:6) 
=M,()|I[+ D(a) |=M,(e) g(v). 
These equations define the matrices D(w) and 4(w). Since the determinant of a pro- 


duct of two matrices is equal to the product of the determinants of the individual ma- 


trices, we have 


[M(w) | =|M,(@) || 4(@) | : (Ze) 


Thus the impurity frequencies are roots of the equation |4(w)|=0. Substituting the 
factorization of M(w) into Eq. (2-4) we obtain finally the following expression for JS, 
the change in S due to the presence of defects in the lattice 


AS = 


ar | f@din 401 (2-8) 


The change in an additive function S which results from a single defect at some 


point @ is 


d5,= + | f(2)d In |4,(2)| (2:9) 
271 c 

and might be referred to as the “ self-S”’ of the defect. The “‘ interaction -S” of a 
pair of defects located at a and /9 respectively is defined as the difference between the 
S of a system of two interacting defects and that of a pair of isolated defects and is 


given by 


a: |Ja.5(Z) | j 
AS. oe | f@ dn old’ (2-10) 

It can be seen from eq. (2-8) that the matrix 4(w) is of central importance in 
the calculation of any of the quantities JS. We shall also show that J(w) can be 
used as a moment generating function and also to find the low temperature behavior of 
the free energy. Thus we can often circumvent the integration implied by eq. (2.8) in 
finding high and low temperature expansions of thermodynamic properties which depend 
on the normal mode frequencies. 

The formulae eqs. (2:8) —(2:10) are independent of the type of lattice being 
analyzed and its dimension. However, the application of these techniques to defect 
calculations requires a sufficiently simple model of a crystal lattice so that the elements 
of the inverse matrix [cf. eq. (2-6) | M,"'(w) are readily obtained. For this model 
we will take a simple cubic lattice with nearest neighbor interactions, both central and 


non-central. The difference equations for the amplitudes un, ny ng are given by 


Mow Un, no ng tr (un, +1, no ng — 2un, no ug + Uny—1, ng ng) 


372 J. Mahanty, A. A. Maradudin and G. H. Weiss 


+72 (un, no+1, ng— 2Uny no ng t+ Un, 3—1) tg) (2s 11) 
SPI (Unynons+ Lb 2Un non, + Unynon3 — 1) =0, 
where M is the particle mass and the 7's are the force constants. Under the assumption 
of cyclic boundary conditions the elements of the inverse matrix M,~'(w) are found to 
be® 
cos 278- (i—j) /N 


-_ sits (2-12) 
~ ae Mo? —2>"7,(1 —cos7s, /N) 


The form in one dimension is exactly the same except that there is only a single s and 
a single force constant. In the limit as N—>co, the sums appearing in eq. (2-12) can 


be written as a triple integral : 


cy 1 \\j cos (i; — j;) Q:Cos (ty —fe) gscos (, — j;) g.do,dy.de, (2-13) 


a3; mes 
cs He Mw? — 27, (1 —cos¢,) — 272 (1 —cos¢,) — 273 (1 —cos¢s3) 


where the integral converges for w >w,> and w°<0, where wy is the largest normal 
mode frequency. 

In one dimension it is possible to evaluate the integral which appears here with the 
result that + 


ROPE OOM fax era 2? ac, 


Met, WV 14F 
where w?,=47M—' and f=w/w,. The matrix elements a{>" for an alternating diatomic 
chain in which the even lattice sites are occupied by particles of mass M, while the odd 
sites are occupied by particles of mass M, are” 


w<0 (2-14) 


aD = — Ws oe [ (oP +*)'? (a2 + 0°)? — (of+ ee?) 4 
——— ee oe an — 
Aap loly w+oy, (w,@.) “~? 


wo <0 (2-15) 
where w7=27M,"' and w,2=0/+ 0°". 

The contour C which appears in eq. (2-8) is usually taken to be the D shaped 
contour shown in Fig. 1. It can be shown is most cases that the ¢ontribution from 
the circular part of the contour vanishes in the limit R->0o (however, see the appendix). 
The resulting integrals are usually quite difficult to evaluate ; hence we will outline some 
alternative methods for finding high and low temperature expansions for thermodynamic 


quantities. 


We begin by rewriting the expression for the free energy in the form 


: bw, b et 
aes ae na sinh ne = a So;+«T 2 In(l1—e «T) 
=E,—«T >} I,/ . 
K — jn (2-16) 


+* 


In the rest of the papep, a;j‘“) will be denoted by g(m) where m=|i—j|. The g(m) are the 
Green’s functions discussed by Montroll and Potts”. 
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where E, is the zero-point energy of the lattice and [, is given by 
T= Sle— ho j/eT , (2517) 
j 
‘We may now write for the change in the quantities 4E, and 4I, 


ABS Fel" fdin |AGF) | 
— Heel" FO(F Daf, 
i a= tf si o( Bee. f oped 


= 1 ["sin (anf) 2D af (2-18) 


where we have used the notation (f) = 
MAS} and @,=nbw,/kT and 
the fact that 2(f) is an odd function of 
f. At low temperatures a, will be large. 
Hence, if we integrate the expression for 
JT, by parts successively we find the follow- 


ing asymptotic expansion 


sea (2020), 2°00) 


Qa a> Gas 


nv n n 


(2-19) 


where we have assumed that (f) and all 
of its derivatives vanish in the limit f—> 
co. When we substitute this expansion 
into the expression for JA we find that 
the change in the free energy is given by 


co i+ Tap? +2 
Ee es (on me ae 
7 n=0 (2n+2)! 
(0, —iR) 
gen ea " B. 2-20 
Figure 1. Contour of integration in the w-plane. x (0)( chin n+l ( ) 


where B, is the n’'th Bernouilli number. This expression is not limited to one dimension 
and shows that the low temperature change in free energy is related only to the behavior 


of Q(f) at the origin. 
The high temperature expansion of the free energy is conveniently obtained by the 


oe bw 
method of moments. When the condition ——<27 we can expand the free energy as” 
K 
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B,( =). h 


bw, Ke “toe ° 
= = In Sy 1)* Slo (2-21) 
= Oly mal 2 ( ) an (2n)! j 


where, again, the B, are Bernouilli numbers. The first term can be written in the form 


a DO & SP a ae 
In = -— Inf] «o; (2-22) 
2 nay ty ee 
pres 
=| +— In |D(0) | 
n() Z 


where |D(0)| is the determinant of the matrix D(0) whose eigenvalues are the w”'s.* 


We also have the relation 


* 
So =N/e 


Zn 
Fel 


=I: | D0} (2-23) 
where 4, is the 2n’/th moment. We therefore obtain 


J4="" in | 4D (0) | KTS) ee a (- : ) NAie (2-24) 


1 (2n) (2n)! 
where Ndt,=Tr|D(0) \*—Tr|D,(0) |". Hence in order to find a high temperature 
expansion of the free energy we need only find an expression for the moments. We shall 
show that this can be furnished if we know the function Y(f) introduced in eq. (2-18). 

By definition 2(f) can be written in the form 


Glee OtyF 5 y=(0) —y; 

2) =— > nl — iS = 2wals > st. 2-25) 
f dful ert e a (w*+¥5) (w*? +4; (0) ) ( 

where the »,(0) are the frequencies of the unperturbed lattice. If w is large enough 

we can expand this in a Taylor series in powers of f~', 


223) 5 (o)" pa m(oy} =2N So "bn 
f r= a ow; fo 3 n=l wo; | Bi ez 


Hence we see that Y(f) is really a moment difference generating function. If the 


coefhcient of 2(—1)"f~©"'” in the expansion of 2(f) is given by a@y,,, then we can 
write for JA, 


(2-26) 


sA=* In| AD(O) |e TS 
Sant 


b 
(ey oe (2-27) 


In the remainder of this paper we illustrate the methods presented in this section 


* The matrix D(0) is obtained from the matrix M(w) by dividing every row and column of M(a) 
by the square root of the mass appearing in the diagonal element corresponding to that row and column, 
and changing the sign of each term. The resulting matrix can be written as D(0) —w? I which defines D(0).. 


' Note that D(0) is independent of w:. 
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by calculating the changes in thermodynamic functions of one-dimensional lattices due to 
the presence of several kinds of defects. 


§ 3. Effect of defects on the thermodynamic properties 


of a monatomic linear chain 


In this section we consider the effect of the following kinds of defects on the vibra- 
tional free energy of a monatomic linear chain: (i) a defect which changes both the 
mass of a given particle as well as the force constants associated with the interactions 
joining it to its nearest neighbors; (ii) a single isotope defect; (iii) a pair of isotope 
defects; (iv) anomalous force constants; and (v) a lattice vacancy. All other ther- 
modynamic functions can be obtained from the free energy by well known relations and 


hence will not be treated here. 


(i) A General Impurity Atom 
If a normal atom whose mass is M is replaced by an impurity atom whose mass is 
M’, and the force constants 7 of its interaction with its nearest neighbors are altered to 


7’, the defect determinant |J(w)| is found to be 


Pe (O19 0) 9-2) i eg) eg Cay eg (2) 


aries) mente (0) =e C1) ly Set) eg (aad) 1 29 (0) — 6.9 (0) 89 1) 
—&9(1)—&9(2) &lg9(2)+9(0)]+269(1)1 —€,9(1) =€,9 (0) 
| (3-1) 


where €,=7’—7 and €,=y7—j'(M/M’). This determinant factors into 


|4(e) |= {1- G7) [9 (0) —9 (2) Ii fig aps wg (1) 
/ 


+( met —M jo" g(0)}. (Ge) 
\ 7 


The change in the zero-point energy of the lattice is given by 


dEy=—"2\ "fa nlf) | 


= Pel" fa In| {1+ (P—J) (v ‘1+f?—f)*} 
an 


bw, (° Q vipa) (3-3) 
=| flint + (PTY 


‘c -1) a on 


/y and Q=M’/M and have taken the D-shaped contour describ~ 


as Pa 
where we have put P=7’, 
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ed in Section II as our path of integration. These integrals have been evaluated for 


the conditions P> be a a - P) >0, Q>0, with the result that 


saat re ePY | 
pye \ Bay = 


A ade A-BS 
(B—a)(a@—-1) (a—f) (8-1) 


27 
a (A— Ba) un (ease) 
CRON Ger ae ae ena he \, a 
3(A— BP) 1 1+)/1—f? 
of I —— =k ee Fr ‘ : 
'(a—8)—f) V1—-# a 3 ) Fy 
where 
P+/ P Sire SSB) 
pe UES eet a lh 
2Q—P—PQ 2(2-* —P) 
Q 
pally AL ae ee rolph 
Cais) ag eee Bi —Q 
20 T= PQ 2(2-— —P) 
Q / 


d 


The function 2(f) = if In |4(if)| is readily evaluated from (3-2) in the two 


limiting cases of small f and large f with the results that 


il = 
2(f) > (PQ—3P+2) +O(f*) (3-5a) 
20 CREE Noh 0-8 


These expressions in conjunction with eqs. (2.20) and (2.26) lead immediately to the 
following low and high temperature expansions for the change in free energy : 


AAT) no = 4B a = ae 
‘J pices ot ee | (3: 6a) 
4SA(T) n2,=— KT In a +«T | z |=(e+ )=1 |e) 

2 Q P42 Od eT - 


= aot ea wk hs dee (3- 6b) 
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(ii) An Isotopic Impurity 

An isotope impurity is one which alters only the mass of the atom at a particular 
lattice point and leaves the force constants unchanged. It is the special case of the 
general defect considered in (i) which results from putting P=j7’/y=1. The expressions 
for the changes in zero-point energy and free energy follow from eqs. (3.4) and (3.6) 
on making this substitution with the results that 


dey ="22.4 (1 8)” cos“! —8) -=} (377) 
2 2 


ae 


€ 7? 
AND pei 4b eo —( aes 


7 bw, 


2 (3-28) (ee eo} (3-82) 


e) 7, 


SAT) 7, = = In(1—€) +«T [20 1 + ol 


“assol 2G 99) 3 ey} 6-80 


where for simplicity we have put €=1—Q. 


Hower, the principal reason for discussing this case separately is that an alternative 
expansion of the temperature dependent part of the free energy can be obtained in the 
following way. 

The expression for JI, eq. (2.18), becomes in this case 


4, == {, sin a, - In [1—&f(1 Hf) df (3-9) 


which for small €(\E|<|) can be expanded as 


Hit é eee df ak fsin On faf+0(é) 
0 (1+/7)*” z Jo (1+f?)? 
4 ee —ancosteos dl — ae LE ELON CE 
az Jo 4 


If we multiply this result by —kT/n and sum from n=1 to 00 according to eq. (2.26), 


we obtain 
&bw,("? cos 4 di 
4A(T) = dE,+ 4 a ih ell Cos olny 1 
> bw, 1 3 (3-10) 
ae pe be . hw (a? 1 Bee OK )s 
ehorlar __ 


an expansion which is valid over the entire temperature range. The integral appearing in 


. . . ion. . 
eq (3-10) can be evaluated as a series expansion in powers of T of I 4) otnustacas 


easily, by numerical integration. 


(iii) A Pair of Isotope Impurities 


es ee 
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If two like isotope defects with masses M’=(1—€)M are separated by m lattice 
spacings in a linear chain the defect determinant |4(w)| is found to be 


Ao) \=1 —e& Mo Fm) — 
| (a) | 11 +EMw*g (0) F 


~~ 2 fe / = aa O(e& (3 . 1 ) 


in the weak defect limit, |€]<1. The change in zero point energy to this approxima- 


tion is given by 


sepa els (PEP HD 


ana Jo 1+f? 
asa Sim _— (3+ |ml) +0 (4+ |m|)| (3-12) 
47 \ 16m*—1 


where 


TQ=4 In (2. 


The temperature dependent part of the interaction free energy can readily be obtained 
using the methods of Section II. In the weak defect limit 2( f) becomes 


; a) Gt pas 
Q =—2M' 0, [ f2g°(m ; if) ] 
(f) a [f?9(m ; if) | 


Bes OBS iloe Phe rhs 
Car 2 & Sen 


__4\m| fi 
a—f): 


Evaluating the small f expansion of 2(f) we obtain for the low temperature expansion 
of the interaction free energy (to O(&) ) 


A ae cia Oy 


a a \3 
AA Giode, = e{24\m| 2 (2 
a 90 \ bo, 


167° Sj aly 

+4. 16 |m|° + 24m? + |m| of : 3-14 
189 kiN : are ( 
This type of expansion is mainly useful for small separations between defects. For large 


|m| we can use the asymptotic form for the Green’s function, eq. 


g (ms if)~— > etm, 
4yf 


so that here 


c2 (a 

Cc . 5 . 
Ap = | Sih Meal Gi 

ae ) 


Jl 


Vibrational Thermodynamic Properties of Lattices with Defects, T: The Linear Lattice 379 


2 3 
= 1284, | ° (| uc Pas ey }. (3 - 15) 
(ages .6m) e816 (@,,? + 16m”)? 
Multiplying this expression by -— and summing over n from 1 to co we obtain 
finally 
AA(T) = dE, — 18 hey, (=) {4\m °F (et 
7 bw, bw, 
b 
—F bnl( Ow, a) F(a) (3-16) 
bw, 
where 
1 {37 7m ~~ cosh mx | 37? 1 
Fi (x) = +43 + 27. 
; 2 ha 4x2 sinh? xx 8x* sinh? zx ( 2) 
pe eee are ae 
ee sinh 7x x° j 
BiG) sade oped alee gps SAR sinb stad |p) 
DMG 2x8 esimla ex 2 sinh 7x xe 


At high temperatures we need not use the moment method in this case as another 
method becomes available to us. In the weak defect limit 


A " ces ae 
—- \, sin & df ig « if) af 


Be “A “fg? (m; if )cos a, f df. (3-18) 
If we multiply this by ee dd interchange the order of summation and integration 
n 


JA(T) = 48,0 ee |" {>} cos a, f Lf? ms if af 
TU 0 


n=! nN 


= dE, — Peek |, frosa("eet— 2nj)—* Earn 14 Pfs 


(3-19) 
where we have used the relation 
21 €0s nx=1>)0(x—27 j) +d. (3-20) 
The 0-functon converts’the integral into a sum with the result that 
4A(T) = —€xT >} = ree sl Lee ig ake (A531) 


where 
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fe ATRL 
; bo, 
For T > ee. the sum can be expanded in inverse powers of T: 
27 
5 bo, B, (m+1) Boy st (bo, 
JA(T) = —€«T (2) i mm sm 
( ) K ‘¢ gimti (4m)! 2*"*!(4m+2)! (> ) 
4 4m’ + 11m+8 Paes bw, ee A (3-22) 


8 21 (4m+4) \x«T? 


Thus in the weak defect limit the first non-vanishing term in the expansion of the in- 
teraction free energy in inverse powers of temperature depends on the separation between 
the defects. This result can also be established in the general case of |€!< 1 either by 
the method used to obtain eq. (3.21) or through the use of the moment-difference 
generating property of (f). This qualitative feature of the high temperature interac- 
tion free energy had been noted previously by Montroll, Maradudin, and Weiss. 


(iv) Anomalous Force Constants 


The next defect we consider is a single anomalous force constant which is the result 
of replacing the force constant 7 for the interation between two neighboring atoms by a 
different force constant 7’’. 


Anomalous force constants have some features in common with vacancies in the lattice. 
We shall first consider the case in which there is a single such defect in a monatomic 
linear chain. 


The defect determinant in this case is easily shown to be 


|4(w) | =1—270,9(0) —9 (1) ] (3-23) 
att 
where dai The change in zero point energy is given by 


bw, {° (f-V1+f?) 
dE, = ee Bes 2 d| F ha SW et As ees 
; 27 a a : a ae Ses 


=f] jie een es oe) ) =| 
= —— — COS s —_—— 
27 Ly 1+20 1+o 2 


NX . te . 
where 0 is taken positive, meaning thereby that the anomalous force constant voy 


is stiffer. 
From (2-28) 


Atzeld Ate, eeepatlg, F a PGi feed) la 
3 \ df Vf+1 f 


al ) ik sin @, df 
otf V1itfto itp) Ga 
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This integral can be evaluated in a power series of 0 when O is small. For this case 


sin a, f a (ery § d 


41,=— ~(' sin Oo df= =i) \; | 
a+fy? CLE (Se een a -bi)? 
+0 (9) 
= b ary] e601 c05 6 d+ e ||; en aecosteos : dF ayer 
T 0 oe , ‘ 
+0(0*). 


Multiplying this by cee and summing from n=1 to ©, we get 
n 


4JA(T) = =e St eee, ae -=| 
27 Ly 1 +20 1 ae 2. 


me /2 
Q bord _ cos Og 


0 hy, [KD cos® ___ 1 


y ae [2 
7 0 


ee cose i IC 


St : | -0(@") : (3-26) 


4 ee AL ae 1 
The high and low temperature expansions are obtained without much difficulty, and 
they are given by 


Loe wes =e he cml =o == | 
20 1+ 20 1+0 2. 


r p ; cae (3-27a) 
oe (a) =(=5)|2 -(<) | aa +3} 


bw,\? (34+0)0 /bo,\4 
AAD) 7p = 3 el Inia + T {2 Ye) eee ae S28 eee’ +o}, 3-27b 
Beans aaeed 1A i-iee KT 11520 a ( ) 


A more interesting case is that of two interacting defects of this type. 
The correction determinant for this case is 


4(w)|\=14+ 0? Mo? 7 g(m)|g(m—1) —9(m+1)] 
Ba ree (1-27 9[9(0) 9)" 
~14+0Mwy9 (m)(g (m—1) —9 (m+1) | (3-28) 


when 0 is small. Here m is the number of lattice spacings between the defects. The 


> : P . i 
expression for the zero point interaction energy 1s (for 0 small), 


me | 


AE Perl" dln |1—0o 
=e ptm (1 


_ bo, lo Eis f2 A\m| (3-29) 
pe ie ams ht ee e 
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citi _ bo, & 1 
2m = 16m?—1 


Proceeding in the same way as we did to obtain equation (3.21) we get the follow- 


ing expression for the interaction free energy : 


ACL tT —V1+077?—Cj™"}. (3-30) 
IV AO 


The low and high temperature expansions of the interaction free energy are (in the 


weak defect limit) 


/ = KL ogg LY Se 4.49) 
LC eaamere tere ert (3-31) 


and 


x9 bo; \'"" (Bom _ (2m+1) Bay +1 bw, \ 
WCE) apd «T(! = oe m1) ) 


2xT 2 (4m+2)!\ 2«T 
4 (6m+1) Bam -2 : ( bw, ) ee | (3-32) 
4 (4m+4)! \ 2«T | 


where the B, are Bernouilli numbers. 
Again it should be emphasized that it is not necessary to restrict the calculation of in- 
teraction energies to the weak defect limit. The methods described here are completely 


general. 


(v) The Lattice Vacancy 


We can apply the results obtained for the anomalous force constants to a calculation 
to the self free energy of a lattice vacancy. We proceed in the following way. A vacancy 
can be characterized by setting the defect mass M’ in eq. (3.2) equal to zero while at 
the same time stiffening the force constant 7’ associated with its interactions with its 
nearest neighbors. This choice for M’ and 7’ resembles a vacancy in a crystal lattice if 
one examines what happens when an atom is removed from a lattice site in a monatomic 
lattice. The repulsive force which kept its neighbors in their places is removed so that 
a tendency exists for them to fall into the hole. We characterize this tendency by cou- 
pling the neighbors of the hole to the hole’s lattice point by a stronger than normal 
spring (7/>7), or by an alternative model in which the neighbors on opposite sides of 
the hole are coupled to gether with a spting constant 7/’=})7/>7;. 


If we set M’ equal to zero in eq. (3.2) and denote the defect force constant by 7” 
the equation for the impurity frequency, |J(w)|=0, is 


(7’—r)[9 (0) —9(2) ]=1 (3-33) 


which is readily solved with the aid of eq. (2.14b) to give the result that 


0 


Ova . cu a 3 . 34 
SGA Sd 
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where. Pa (77/7). > 2. 


The contribution to the change in free energy from this impurity frequency is 


bor, P 
Kl 2 P= 
ie P a 
Se Per eee es rm 2-35) 
el 
Ik 3G 


The contribution to the change in free energy JA, arising from the small shifts of 
‘the in-band frequencies can be obtained as follows : 

(a) We insert an anomalous bond with force constant 7’ into the lattice. The 
change in free energy 4A, arising from the in-band frequency shifts is the difference 
between the total self free energy JA,;,. given by eq. (3:27) and the contribution 


from the localized mode whose frequency is given by the solution of 
2G” —7)[9(0) —9 (1) J=1. (3-36) 

‘This impurity frequency is found to be 

(PE 


Doe ew 2 Beis) 
Pilibgesco2Ri an \ 
where P’= (7/7). Thus we find that 
AA,OUT) Z4A(T) 24, — 2 $0 | ===" 
b a.f.c. 2 L ae 
tial, WP" 
oer VOPR ei 
ee 
—«T In : : (3-38) 
) W 7, 
a eats 
LG 
(b) We now remove one particle from the chain with 
/2 g Pie sin 0 
ap vee ee Mec | ln (lgeeat add, (3-39) 
8 Tw VO 


‘since this is the total free energy per particle in a long chain. 

(c) Without changing the topology of the lattice the anomalous bond is chosen 
to span the vacancy so that it is extended over two lattice spacings rather than the usual 
one per bond. The choice P’=% P is made to correspond to a consistent model of a 


vacancy. The total self free energy of a vacancy in a linear chain is then given by the 


sum of (3.35), (3.38), and (3.39): 


| ba 
DAC) cas SAAD) ate. ba, 


Plat [2eP 7 
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bor, te 


—*T ("in (1—e eT dd. (3-40) 
0 


— 
aod 


Low and high temperature expansions for this energy are given by 


_ bw, £. aay 2=P rag "—2| 
AACE) pe =< | ( P te 
a; en _ =e | “eT ) 
a nat {-—(74 NG os, a ae oa \ bw, 
ay Add he Ft (tan, “iar 
4A(T) wp. = ¥ KT In io? +x] =: ia ) ( a) 
salt te ool ebay aM. (3-41b) 


46080 «KT / 


§ 4. Defects in diatomic lattices 


With the aid of the “‘ M*” transformation and the expansions developed in Section 
II, the calculation of the effects on the thermodynamic properties of diatomic lattices 
becomes almost as easy as the corresponding calculations for monatomic lattices. 

We consider first the case of a single defect which alters both the mass of a par- 
ticle and the force constants associated with its interactions with its nearest neighbors. 
We assume that the particle at the origin whose mass was originally M, is replaced by 
a particle of mass M,(1—€). Furthermore, its force constant 7 is changed to r({1—y). 
We consider only the case M,>M, since the discussion in the case M,< M, follows along 
identical lines. 

The equations for the time-indepent parts of the displacement components are given 
below. 


M,(1—€)w*u(0) +7 (1—7)[u(1) —2u(0) +u(—1)]=0 
M,w*u(1) +7[u(2) —u(1) ]—7 (1—) [u (1) —u(0) J=0 (4-1) 
M,w*u(—1) +7 (1—7)[u(0) —u(—1) ]—7 [u(—1) —u(—2) |= 


“ce +> 


We now go over into the “v” notation described by Montroll and Potts. The 
transformed equations are 


M*w*v(0) +7[v(1) —2r(0) +v(—1) ]= (M, cat — 279) Se ate v(0) 
oO 27 . 


M*o*v(1) +7[v(2) —2v(1) +v(0) J=n7 {v(0) — at ae 11) (4-2) 


My" — ne 


M¥otv(—1) +7[0(0) —20(—1) +0(—2) =r {v(0) — ae v(—1) 


so that the determinant |J(w) |=|I+M,~'(w)0M(w) | becomes 
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t +900) Fie Ae —9J (0) er 9 (1) €.—9'(2) € =9 (1yerrg (2) 

| 9(—1)&:—9(0)€  1—9(—1) €,+9(0) €,—g(1) €, —9(0)€,+9(1) €, 

| Oi 2); —9 (1) €, —9(—2) &.+9(—1)¢,—g(0) €,51 —9(—1) €,--9(0) €, 
(4-3) 


where 


., Mo S27 ie 
7M, wo 7)" 


el A/) (4 3 4) 


— 1/2 
Se Si See 
(mM, wo Por) ir 


This determinant factors as follows : 
|4(@) |= {1+9(0) €:—9 (2) &} {C1 +9 (0) 5-29 (1) €&)) (1+9(0) €, 
291) 4-9 (2) &) —29(1) 9 ©) €5)- (9. (2) 9 (1) eg +9 (0) €) }, 


(425) 
1 
In the special case of an isotope defect (€,=€,=0, €¢,=—M,ew° EEG pee b ie, this: 
(Mw? —27)'" 
reduces to 
M,w?— 27)" 
|\d(w) |\=1—M, ew Dae 9) : (4-6) 


The use of (4.5) for finding changes in thermodynamic quantities leads to cumbersome 
integrals. The essential features are there in the isotope case; so we will consider the 
isotope case alone and use (4.6). 

For the evaluation of the low and high temperature expansions for the changes in 
thermodynamic properties due to the presence of the defect we need to evaluate the 


small f and large f expansions of the function 
oot ae nls) 


respectively with | a Using the expressions for the /(k) given in eq. (2.15) we 
Oo; 


obtain for the leading terms 


2(f) —{— i oe 
1 


a 
(O} g 2( 0 205 — Of 0, + oso, 4 
ee = 2 <A 1 1 i erage (4-7) 
iS A (OF Of os 


Tea aiuapeer pails: c)+| (4-8) 
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- where 


SMES Ps 2 
at (1ta)?+ 22> 20°=2) 
4 4 


1 
Canal | (#+32") (7? —a?*) ge i lias, : 
8 


We can immediately write down for the leading terms in the low and high temperature 


“expansions for the change in free energy : 


(AA) py = SE, +E Sof SE )+ ies Bei ) = pat (4-9) 


60, bw, 450, \bw, ¢ 


(4A) TA ee < In(1—<€) 


+ T|—— (722) 1 eB ER fos |. 
ae 2880\1—e =1—€/\«T 


(4-10) 


‘It remains only to evaluate the change in zero-point energy, JE,, which is given 
‘formally by equations (2.18) and (2.15) as 


f £2 2 \1/2 
4E,= — ae dal oaaect = (f+? ie (4-11) 
27 yo V14+f?\ff+eé 


‘This integral has been evaluated by Mazur, Montroll and Potts and by Montroll and 
Potts” in two limiting cases, M,=¢~* M,, © <1, and M,=(1+7)M,, 7<1. To terms 


4 


linear in € and 7, the results are 


-log 
JE\= } bo, {t[—}+ =e) ohne +sin7'€) ]+0(e*) +--+} (4-12) 
1 1 3 6 . 
dE,= A bo,|— as = (1—€*) “7 (742 sin“! €) 
) 
{24/2 In(1+)/2) —\/2e"+ t-+2 sin te foto a 
gz(1+e 2) V 1 (= ee sin €) |+ f 
(4-13) 


The more interesting case of a defect pair can be treated in the same manner as 
‘the monatomic case. We shall consider here the case of isotope defects only. There 
care three distinct possibilities : 


(a) Both the defect atoms are on heavy sites, M,’=M,(1—€«,) 
(b) Both are on light sites, M,’/=M,(1—«,) 


(c) One is on a heavy site, and one on a light site. 
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As in the monatomic case, only terms up to € can be evaluated without difficulty 


so that the treatment we give here is valid for € small 


(a) Defect pair on heavy sites : 
The correction determinant in this case is easily shown to be 


|d (iw) |=1—e comm” To JM (m, iw)? (4-14) 
Oe OF 


where ¥ is given by (2.15) and m in the number of lattice spacings between the defects 


The expression for the interaction energy is 


a ' 
AE, 2" dinfi—ete MyM (oe )9 cm m)?| 
~~ €, ede F T+o;y = gm)" dz (4-15) 
2m JO EN 


where we have integrated by parts and retained terms up to €,. The integral can be 


evaluated with some approximations” and the value is 
al fz. 2 89 
AE, 1 poy € 76MM, ‘ (4-16) 
3 4x (M. i+, ) Ray 2 2 


m 


The low and high temperature expansions of the interaction free energy are obtained 


as before, and are given by 


ree KT Y ie 
\bw;, 31 of 


4A = JE eee 
( dz Mig 0 = AS we 
i ys ie (4-17) 
bw, 
ey? (w Wo) “| 4 Don (yy 2 Bory 2 Oy ae * ; 
=e Se s (4-18 


(b) Defect Pair on light sites. 


The correction determinant for this case is 


|4 (io) [1 — efM Myo" ) gy (m, iw)? (4-19) 
wo +s 
‘The expression for the interaction energy 
dE,=—** oa | Ato yen m) dz 
47 -0 \ 2+," 
LY; l (i FOP es, "pM, M, (4-20) ° 


m 4n (M,+M.,)*°? of 
The low and high temperature expansions of the interaction free energy are 
78 2 
(4A) pe = 4e,— 27 |* "Oy sf “TD \- E [ (2m? —1) w,3+02%0)7 +.0'] 
he 450. \ bw, 315a,! 


8 
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x(EY—) (4-21) 


} 
bow, 4 


9 \ im=2 
Se ie (Os) EF im 2Bomsi0z_(_% _ | « » (422) 


(c) Defects at heavy and light sites. 
The correction determinant for this case is 


| A (iw) | 1— €,€.M,M,o%g (mio) > . (4-23) 
The expression for the interaction energy is 


Hf Rapa a see ee |e 9 (m) dz 


ee 1 7" €,€bM,M, (4-24) 


— mm? 4nz(M,+M,)>? - 


The low and high temperature expansions of the interaction free energy are 


{ —6 vs 4 5 
(4A); > =4E,— oF (rts 23) 7 =| 2m wy, (7 Froid (4-25) 
wx 


T 3t5 ofos \ bo; 
(4A) ii _ KT €,€, (WW) a Bony (22) Bom + 1 wo, rpm BN it Ps \ (4 . 26) 
play (4m) !\ «eT (4m+ 2)! \«T/ } 


§5. The effects of a large number of impurities 


Let us suppose that there are n defects in a lattice and that we are interested in 
<S(n) >, the value of a thermodynamic function S averaged over all configurations of 
n defects in the lattice. By a well known formula from the calculus of finite differences 


<S(n)) can be expressed as 


(S() =3i(7)K4'5(0) (5-1) 
where 

(dS (0) > =¢S(1) —S(0)> (5-2) 
and 

(A*+S(0)) = CA (*S(0))) (5-3) 


If we denote by W(r,,1r.,---r,,) the joint probability that there are defects at rj, To,°:- 
r, then 


(45(0) )= SW (r,)[S(r,) —S(0) 


a) ae 2 Wr, r2)[S(ry 72) —25(r,) +5(0) | (5-4) 
rere 
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where S(r,, r,,---r,) is the value of the thermodynamic function for defects located at 


Ty) Fo.°° Pp. 


We define c, the concentration of defects, by c=n/N. If, as N->co we have 
y iS 


where the S, are constants, and n—co in such a way that c remains constant, we may 
write, using eq. (5-1),” 


’ Se 6 
(S(n) )=S,+S,e+ plies . (5-6) 
If S is an additive function of the normal mode frequencies, as in eq. (2.1) then” 
{ AF ee ey | = x9 ® 
(4800) = 7 _< {FAD (—1) (in |4,(0 |) (5-7) 
LIN IC ja! i 


where |J;(z)| is the correction determinant associated with j defects. 
The problem of averaging is a matter of some difficulty because of the interactions 


between the defects previously described. At sufficiently high temperatures one can ignore 
the correlations so that 


W (ry yarn) =1/(p ) (5-8) 


For example in this approximation the ‘‘interaction-S”’ of two defects is given by 


42 = 2 | 4p (z) | d 
a (| foam VRE (5-9) 
‘One can thus evaluate the first few terms of eq. (5.6) without difficulty. 

For example if we allow the position of defect atoms to be uncorrelated, i.e., if we 
use eq. (5.8) to express W(r,) and W(r,, r.), then we find for the free energy per 
A(T) ) 

N 


particle, ¢ , at high temperatures for a concentration c of isotope defects : 


AT)? yt T In. bw, Es T| 1 eas 1 (Fe) | 
N 2«T 48 \ «T 7680 \ «T 


tsladnent kava 7 Aer? 
ae - aka) on 


23,040 \ 


where we have used eqs. (3.8b) and (3.22). 
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Appendix. Change in frequency spectrum due to a single isotope defect 


The mathematical methods outlined in Section II are general enough to be applica- 
ble to lattices of more than one dimension. However, in order that the D-shaped con- 
tour described in Section II can bz used as the path of integration in evaluating JS the 
function f(w) (cf. eq. (2.1)) must satisfy several conditions. Inasmuch as J(w)) is 
a function of w” only, the logarithmic derivative = In| J(iw)! of |d(iw)| evaluated 

w 

along the imaginary axis is an odd function of w. In order that the integration down 
the imaginary axis should not vanish, f(i) must either be an odd function of w also, 
ot it must have an odd part. Secondly, and as important, in order that the contribution 
to the integral from the semi-circular portion of the contour vanish in the limit as R-> 
co, it can be shown that f(w) must behave like at most |w* as w/— 0c, where @< 2. 
When these two conditions are not satisfied we must either re-express f(w) in an alter- 
native form or modify the contour. (It is these considerations which led to the decom- 
position of the free energy given by eq. (2.16).) In the case of the characteristic fune- 
tion of the squares of the normal mode frequencies which is defined by f(w) =N™~™ exp 
(iaw*), these two conditions are not satisoed. Choosing instead Laplace transforms for 
which f(w)=N~ exp (—a@w) leads to very difficult integrals, but appears to be the 
only alternative in the two and three-dimensional cas2s. However, for one-dimensional 
lattices, a special transformation of the integral repres:nting the inverse matrix elements 
of Section II leads to an interesting method for evaluating changes in additive functions 
of normal mode frequencies in the presence of defects. We give here an outline of the 
method and apply it to the evaluation of the change in the frequency spectrum of a 
lattice containing a single isotope defect. 

If we introduce the transforma- 


tions* 


@— plane Wo 


(A.1a,b) 


the one-dimensional form of the ex- 
pression (2.12) for the Green’s func- 
tion g(k) becomes 


‘ : : f 1 i ys ty * 
Figure Ala. An alternative contour of integration in g(k) =— —$ - L => dy 
the w-plane AT 1 J), =1 (y—€) (y—¢- ) 
(A-2) 


* T. Tanaka, unpublished work. 
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I 5 -E4 s 

= ae for+|€| <1 (A: 3a) 
Looe ST = 

ae imran 3A for \¢ | et (A; 3b) 
/ $ Sal 


The transformation (A. 1a) maps the contour C in the w-plane (see Fig. A. 1a) 
into the contour /’ in the ¢-plane (Fig. A. 1b) where the dashed curve is the unit 
semicircle. The expression for JS given by eq. (2.8) becomes 


| 1 fo) (SSE ———— ee 
dS= | — 2—F—F" | A(z Ce 
27i site BA Sy 2 IE (A-4) 


We now apply this formalism to the evaluation of the change in the frequency 
spectrum of a monatomic linear chain due to the presence of a single isotope defect. 


In this case f(w) =N~ exp (iaw*) and we define JF (a) by 


: if sore 
AF (ol) = | ee 6d in |d4(@)-. A:5 
@= 4 (0) (A-5) 
4F(a@) is related to the change in the distribution function of the squares of the normal 
mode frequencies, IG(w*), by 
AG (u?) = 1 |” dF (a) e-*** dev (A-6) 
PAIS Jv oO 

G(w") is related to the usual distribution function for the normal mode frequencies,. 
g(w), by J(w) =20G (w*). 

The determinant |4(w)| corresponding to a monatomic linear chain with one isotope: 


defect is obtained from eq. (3.2) as 


Figure A. 1b. Transformed contour of integration suitable for one-dimensional 


problems. 
| (w) |=1—€Mw*¢ (0) (A:-7) 


and in the present case becomes 
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a =¢—E") 

JV 2-8 -F") — ) Ti (A-8a) 
Aes : 

SF a iva lane BS ag (A- 8b) 


To carry out the integration of (A-5) we divide up the contour I” into suitable por- 


tions as shown in Fig. A. 1b and use the appropriate value of cede yale on 
each portion. 
a) Integration on /’, and I’,. 
On I, €=(1+0)e* and |€|>1. In the limit as d—>0 
sy OR alae 
(GF) n= I ey ee la +es | (A- 9a) 
271i v0 isin 
On I, €=(1—0)e* and |¢|< In the limit as 0—>0 
(4F) r= wf PO sini2 (0/2) J Inj| _¢ =e = |. (A- 9b) 
27 Jo isin 
Adding these two results we obtain 
(4F) nen=—| eit? sin? (8 PT jd In| a et os] | : 
JO : 1 sin @ 
age i eit)? sin? (9/2) cos 6=4 
a Jo sin’ #+&(1—cos #)* 
= 1 i eit 7? sin2(@/2) da ’ (A: 10) 
20 (1— (1—&) sin*(@/2) 


b) Integration on /'),: 
We divide this into two parts, /’,"? where OX 4 (7/2) and /’,” where 7/24 


mt. In either case, ¢=—(1+0e'). 
For 1, , |€|>1, and in the limit 0-0, 


“ae /2 
1 | exp (iaw,*) (i) dd= —1 exp (iaw). 
Ti. 4 


(AP ies 
2710 
For 1’, ||{1, and in the limit 0>0 
S 5 thi a ; 
CAP) og = exp (iaw,”) (—i)d#= —- ; exp (iaw;°). 


Hence 
ape Ree 
Pn Sy op Gaz), (A-11) 


c) Integration on /’,: 
In this case, by the substitution ¢=1—Jde"*, O<A<z it is easily shown that for 


0-50, the differential coefficient of the logarithm vanishes. The contribution from this 
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source therefore vanishes. 
Eee 


d) Integration along the rectangle /’,: 


The contribution from this part of the contour is simply the sum of the residues 


at the poles of the integrand inside the rectangle. The poles are at the roots of the 
equation 


(ee 


—=0. (A-12) 


There is only one root to this equation in the interval (—1, 0) and this is given by 


i=" with O<€<1. If €<0 there is no root inside the rectangle. We thus 
+1 
obtain 
CAE.) = eitw pl €2) QetS ik (A- 13a) 
=200) c= 0: (A-13b) 
Adding all of these results we obtain 
JF (a) oe |E| \ exp (ao; sin” (G/2)) ap 1 saw 724 tow p2/—e2 ) 0o< E€ ih 
2mN Jo 1—(1—€)sin® (0/2) 2 
(A- 14a) 
Seclae [iseGtes sin’ (9/2)) yg 1 iaw,? 80, (A-14b) 
2nN J01—(1=&)sin' (0/2). 2 


If we substitute these expressions into eq. (A. 6) and interchange the order of 
integration we obtain 
dG (vw?) = —- lé| a(S: sin ey ae) di) 
2mN Jo 1— (1—€*) sin? (6/2) 


OE ee ~O( pi —w) wa (A-15a) 
2N ; N \1-& 


Mah {- O (a7? sin’ (0/2) — 0) gy yoo") e= 0; 
N 1— (1—&) sin? (7/2) 2 
EOS alee at (A-15b) 


where 0(x) is the Dirac delta function. 


Using the following property of the delta function 
si f%) 
| f(@® o[ [9 (x) Jdx= Dap J (x) 
where the summation extends over all the zeroes of (x) in the integration interval J, 
we obtain finally 
|é| oO; 


dG (ow) =— aN Fenaate 24 (@—1) wo 02-0 


w * ; 
a yon at) + ao On a sw") O<F< 1 (A: 16a) 
2N 
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Seael) cess Op PB Pe aty OER ee 
2anN [w2+ (2-1) 0 of o7—-o&% 2N 


The above expressions show that for positive €, 1.€., 2 lighter isotopic defect, the 


inband frequencies are raised by an amount of the order of 1/N and a discrete frequency 


appears at w=". For a heavier defect only the inband frequencies are depressed; 


no discrete frequency comes out of the band. 

Recently Maradudin and Weiss* have presented a method for evaluating successive 
terms in the power series expansion of /(w) in the neighborhood of w~0. The method 
is general and is not restricted to the one dimensional case. The result of their treat- 


ment is that if we can expand 4(w) about w=0 as 
g(w) =a, +a,0?+4,0'+ >> (A-17) 
then the coefficients a. are given by 


(=1)" 20), 


=_\ (A-18) 
Nz(2n)! wo 2" 


where 2° (0) is defined by eq. (2.29). 
In the particular case of a single isotope defect in a linear chain for which 


& 


2(f)= want ; (A-19) 
3 (24+-6) * GOP) 
eqs. (A. 18) and (A-19) give us that 
cK 3) 5\ um = & E> w* 
dg(w) =— Ish ae aC) eat oea ee 2 fenssh, A-20 
( ) Naw, | & ) oOo; \ e 5 ‘ ) 


If we obtain dy(w) from eq. (A-16) by the relation Jy(w) =2wJG(w*) and expand 
this expression in powers of (w/w,) we get the expression (A-20) exactly. 

It may be mentioned here that contour integration in the ¢-plane instead of the 
w-plane can conveniently be used for evaluation of other additive functions of normal 


mode frequencies for a linear chain with defects. 
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-On the Consistency between 
the P-wave Dispersion Relation 


and the Experiment 


T. Tsuchida and A. Kanazawa 


Department of Physics, Hokkaido University, 
Sapporo 


July 14, 1958 


Recently Puppi et al.’? have pointed out 
a discrepancy between the dispersion relation 
for the <~—p scattering amplitude and the 
experiments. In this note we reinvestigate 
the dispersion relation for the (3,3) P- 
wave amplitude with new data available 
now and show that there is a certain dis- 
crepancy similar to Puppi’s case. 

For this purpose we follow the analysis 
by Cini et al.’ in their determination of 
the coupling constant using the Chew-Low 
equation, since the P-wave dispersion re- 
lations?) are reduced to the Chew-Low 
equation with the point source in the static 
limit. Cini’s analysis involves only those 
data at energies higher than 100 Mey (meson 
kinetic energy in lab. sys.) and we can 
now avail ourselves of some data at energies 
lower than that energy. According to Cini 
the dispersion relation for the (3, 3) 


amplitude 1S approximated to the relation 
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| SiN Oy, COS Ong 


— 2X (C7) | (Op = : a On, 


(1) 


with 


m 
3 


ome Ley 


Opie OF 250 --@, 


X(w,) =1P | a LG (eo 
ve v1 
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; : 
+< g(op) +4 gi(e,) | (3) 


where g and w, are the momentum and 
the energy of the meson in the center of 
mass system. f° is the rationalized re- 
normalized coupling constant. If we restrict 
ourselves below the resonance energy, eq. 
We can 
check the consistency if we substitute the 
empirical data for the l|.h.s. of eq. (1). 
In calculating X(w,) we use Anderson’s 


(1) holds with great accuracy. 


empirical formula 
tatiidg==0.2375p°/(1+420.75p) (1.95220), 
(4) 


for the energy less than 300 Mev. The 
numerical calculations are done with great 
care, because there is a principal integral 
in X(w,). The results are shown in Figs. 
1 and 2. In Figs. 3 and 4 we plot the 
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1.0 1.5 20 : 
Fig. 1. a) w7X (we) as obtained from nace: Fig. 3 wf Seta as “yy ths one be ice 
‘cal i i culation. 
cal calculation. b) The experimenta cal ca : weeny 
values*) of the quantity wy (sin 633 cos 633) . values?) of the quantity w,* (sin ds; cos S33) - 


A“ a) f**=0.003+0.009 
. Y*=0.124+0.012 
Pa 


- 0.2 


a) f?=0.033+0.007 
Y=0.074+0.009 


b) #*=0.108+0.008 


Y* =0.024+0.011 
b) f2=0.092+0.005 


Y=0.016++0.007 


Fig. 2. The difference of curve a) and points b) in Fig. 4. The difference of curve a) and points b) in 

Fig. 1 plotted versus wy. Fig. 3 plotted versus w,*. 

a) f*=0.033+-0.007, Y=0.074-40.009, evaluated a) f° =0.108-+0,008, Y*=0.024+0.011, evaluated 
with a weighted least square fit to the points withh weighted least square fit to the points 
below 100Mev (4~G). below 100 Mev (4~G). 

b f2=0.092-0.005, Y=0.016-+0.007, evaluated b) f=0.003-40,009, Y*=0.124+0.012, evalua- 
with a weithted least square fit to the points 


ted with a weighted least square fit to the points 
below resonance (A~M). below resonance. (4~ M). 
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results obtained by replacing w by w*=a 
+k/2M, as is usually done in the analysis 
of the Chew-Low equation. 

As is seen in the figures, if we use all 
the data below the resonance (A~M), 


we obtain 
f’ =0.092 + 0.005 
| with Y=0.016+ 0.007, 
and f*?=0.108+ 0.008 
with Y*=0.024+0.011, 


where the symbol * means the replacement 
w—>wo*. The values are consistent with 
the one authorized today. If, however, we 
restrict us in the energy range less than 


100 Mey, we obtain 
f?=0.033 + 0.007 
with Y=0.074 + 0.009, 
and f*?=0.003 + 0.009 
with Y*=—0.124+ 0.012 
These values are far less than 0.08. The 


above constants are determined by the 
method of weighted least squares. The 
errors are the standard mean square devi- 
ations, which do not involve explicitly the 
Thus, we 


can conclude that for the P-wave amplitudes 


experimental errors in itself. 


there also seems to be a discrepancy similar 
to the one pointed out by Puppi, so long 
as we use Anderson’s empirical formula in 
On the other hand, 
there are many facts indicated that the 


calculating X(w,). 


theory would reproduce the empirical facts 
quantitatively so far as only P-waves are 
Then we 
believe that this discrepancy would not 
lead to the immediate break down of the 


concerned with low energy. 


micro causality, since the integral X(«,) 
seems to be sensitive to the derivative of 


the curve of 0,, and the discrepancy may 


vanish if we use the modified formula, 
increasing more rapidly than Anderson’s 
at 100 Mey., which is the same situation 
as is indicated by Lommon et al.") on 
We wish to thank Professor 


H. Tanaka for valuable discussions. 


Puppi’s case. 
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5 (1957), 1305. 
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102 (1956), 1174. 
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Phys. Rev. 109 (1957), 1750. 
C; Barnes, 7th Rochester Conference. 
D; J. Orear and J. Lord, Phys. Rev. 93 (1954), 
SW 
E; D. Bodansky, A. M. Sacks and J. Steinberger, 
Phys. Rev. 90 (1953), 997. 
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H;; J. Orear, Phys. Rev. 96 (1954), 1417. 
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On the Radiative Corrections to the 
Scattering of Mesons in an External 


Electromagnetic Field* 


Anand K. Bhatia and Joseph Sucher 


University of Maryland, College Park, 
Maryland 
June 3, 1958 


* This research was supported in part by the United 
States Air Force. 
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The radiative corrections to the scattering 
of a spin-zero particle in an external electro- 
magnetic field were considered some time 
ago by Kinoshita and Nambu’. These 
authors used the Duffin-Kemmer formal- 
ism,” thereby also obtaining the correspond- 
ing result for spin-one particles. In view 
of the complexity of their calculations it 
was felt desirable to recalculate the effects 
in question for the more important case of 
spin-zero particles, by a different method. 

Such a method is available in the Klein- 
Gordon formalism.” On computation one 
finds that the second order radiative cor- 
rections lead, in first Born approximation 
in the external field, to a modification 
jul (x) =f (*) +67, (x) of the free meson 


current j,(x). Here 
\: at. a By of / ~ae cy 
aj.) = | dx’ j,(@’)|(—24+2) 


XE (x2) +(log ”—) O° (F,(x—2') 


hae 


+E, (e=2))} = 7 F(a’) | (1) 


where the notation is as in reference 1. 

Equation (1) differs from the correspond- 
ing result of Kinoshita and Nambu (see 
ref. 1, eq. (140) and eq. (65)) by a factor 
of 2 in the term proportional to E(x—x’). 
This discrepancy may be traced (with some 
MEG LISS) of oref. 
1. When this is rectified the calculations 
are in complete agreement, as is to be ex- 
pected on theoretical grounds”. 

One of us (A.K.B.) would like to thank 
Professor R. C. Majumdar for suggesting 
this problem and for helpful discussions ; 
he is grateful to the University of Delhi, 
where this work was begun, for a research 
assistantship.. We appreciate the interest 


labor) to an error” 


that Professor T. Kinoshita has shown in 


this work. 


1) T. Kinoshita, Prog. Theor. Phys. 5 (1950), 473. 
T. Kinoshita and Y. Nambu, Prog. Theor. Phys. 
5 (1950), 749. 

2) N. Kemmer, Proc. Roy. Soc. A173 (1939), 91. 

3) See, for example, F. Rohrlich, Phys. Rev. 80 
(1950), 666. 

4) On computation of eq. (133) from eq. (132) 
of Ref. 1, one finds that the term written as 


“ae 3 (-e) ) 


Al 


(4+ 


should be divided by 2. 
Ref. 1, eq. (136), for the vector case should 


The results given in 


also be changed accordingly. 
5) K. V. Roberts, Phys. Rev. 83 (1951), 188 


On Parity Nonconservation 


Takao Tati 


Department of Physics, Kanazawa University, 


Kanazawa, 
June 6, 1958 


Parity nonconservation is established in 
various decay processes. Questions arise ; 
why do there exist parity conserving inter- 
actions and nonconserving ones, what does 
parity nonconservation mean and is it a 
characteristic to the weak interaction as was 
expected ? I imagine that parity noncon- 
servation has its origin in the stereographical 
structure of the form factor in the primary 
interaction in which four particles participate, 
in a circumstance similar to that in the 
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structure of d and / isomers in stereo- 
chemistry. I have not yet reached to a 
concrete conclusion, but the following fact 
seems to me to suggest the idea and to 
have a phenomenological meaning by itself. 

The experimental results, (a) the asymme- 
try parameter (a) is large in 9 decay, 
decay, 7 decay and A->p+7~ decay and 
(b) parity is not conserved in 3 * decay 
while the up-down asymmetry is very small 
in + decay" (although the observed values 
of asymmetry coefficient (Pa) are not yet 
accurate), can be explained if we assume 
that (1) the interaction of four Fermi 
particles (including Hyperons (Y)) is the 
only primary weak interaction, (2) parity 
is not in general conserved in the transition 
by this interaction but parity is conserved 
in the exceptional case in which there are 
indistinguishable particles—particles to be 
excluded by the Pauli principle—in the 
particles outgoing from or ingoing to a 
(nonlocal) vertex of this interaction. | (It 
is the case when a and ¢ are identical 
particles in b—->a+c+d.) and (3) the 
interaction (N+N->7) is considerably 
stronger than the interaction (Y + Y=37). 
(The assumption (3) makes the discussion 
very simple and definite.) 

These assumptions give the following 


decay processes. 
Lae pet 
Ue al 
7c —> p= i—> py 
i > p+n>p+y 
A> ptn+p->p+a 
S*+—>n+ptnn+7* 


In these cases, no indistinguishable particles 


are present in every stage connected by the 


weak interaction, hence the asymmetry 


parameter is large for these cases. 


Il. a n+p+p 


A iia ie 
Crone 4 


; pee en 
V+ via Seta 


™ ptp+? 
In these case, both the parity conserving 
process (lower process) and nonconserving 
process (upper process) are possible. For 
these cases, the asymmetry parameter is 
expected to be relatively small compared 
with those of I. 


HI. 27> -n+n+pon4+27- 


In this case, there is only a parity con- 
serving process, hence the asymmetry para- 
meter is expected to be very small. 

In view of field theory, it may be 
difficult to construct the desired form factor 
which gives the interaction considered, since 
the interaction Lagrangian must receive 
certain severe restriction originated. from the 
assumed frame of space-time. In this con- 
nection, a new description of the dynamics 
of elementary particles might become neces- 
sary. In field theory, the system of 
elementary particles are represented by a 
system of some operator functions of the 
parameter x and the dynamics is the pre- 
scription of the variation with respect to the 
parameter x, of the state vector or the 
operator functions. (In this sense, the 
parameter x has the meaning of space-time.) 
I think it is reasonable to formulate the 
dynamics of elementary particles without 
refering to the parameter x.” In this case, 
x is a kinematical parameter (relating to 
the momentum of particles) and has not 
directly the meaning of space-time (Contrary 


to field theory, space-time is derived from 
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the dynamics as a statistical quantity), and 
the restriction for the interaction operator 


becomes looser than that in field theory. 


1) Crawford, Cresti, Good, Gottstein, Lyman, 
Solmitz, Stevenson, Ticho, Phys. Rev. 108 
(1957), 1102; Eisler, Plano, Prodell. Samios, 
Schwartz, Steinberger, Bassi, Borelli, Puppi, 
Tanaka, Woloshek, Zoboli, Conversi, Franzini, 
Mannelli, Santangelo, Silvestrini, Glaser, Graves, 
Perl, Phys. Rev. 108, (1957), 1353. 

2) T. Tati, Nuovo Cimento 4 (1956), 75, Prog. 
Theor. Phys. 18 (1957), 235. 


On the Posijronium in Metals 


Hideo Kanazawa 


College of General Education, 
University of Tokyo, Tokyo 


June 27, 1958 


The experimental results indicate that in 
metals positrons reach the thermal velocity 
in a short time and subsequently remain 
free, annihilating with conduction electrons.” 
However, the possibility of the existence of 
positronium in metals, undergoing continu- 
ous exchange of the electron partner may 
not be excluded. Though this problem is 
a many-body problem similar to that of 
the exciton, we consider a system of an 
electron and a positron and investigate the 
Properties of the positronium in metals. 
The properties of the positronium in a 
metal are somewhat different from those of 
the positronium in free space. The dif- 
ference is due to that the interaction between 


an electron and a positron in a metal is 


not the Coulumb interaction, but a screened 
short ranged interaction.”’ 

Choosing the Sommerfeld free electron 
theory as a starting basis, we have, for the 
non-relativistic Schrédinger equation for the 


positronium in a metal 
a an 
(E-e)yazy, 


and m, 


where L=2m,m,/(m,+m,), m,, 


being the masses of the positron and the 


electron respectively. For the ground state 


(S state) (1) reduces to 


where R(r) =¢5(r) /r, S=qr, b=e w/b q 
and 7=Eyb°q°. The screened potential 
e~*/¢ can conveniently be approximated by 
Hulthén potential e~*/(1—e). Then the 


normalized eigenfunction is given by*® 
(7) = {6(F—1) 9/87} "2 
Ae ad oe ee ee 
and the eigenvalue is 
E=—#'q°/4u- (1—e 2/bq)*. (4) 


The lifetime for the two-photon decay is 
given by* 


T=1/Rer" |Y'(0) |? 
=1.25 X107"/8z4,°|(0) |* sec. (5) 


where rp and a, are the classical electron 


radius and the Bohr radius respectively. 


Table 
en ene. Se 
Metal mt/m —E try) rote A/2a9 
Li 145 | 034 267 | O72 
Na 098 014 075 105 
K 093 015 066 1210 
Rb 089.) O14 | os7 | 226 
Cc 0.83 | o16 | 944 | a8 


1 EE ee 
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In the Table we have computed the binding 
energy and the lifetime, assuming m,—=m, 
=m*. The binding energies of positroniums 
in alkali metals are smaller than that of 
the positronium in free space (0.5 ry), and 
the lifetimes are longer except in Li. 

4 is defined by (4) /b(0) =e"". q is 
determined by the relation’ ¢/k,=0.470 
r,/*, where k, is the radius of the Fermi 
sphere and r,= (3/4mna,*)", n and a 
being the electron density and the Bohr 


radius respectively. 


1) S. Berko and F. I. Hereford, Rev. Mod. Phys. 
28 (1956), 299. 

R. A. Ferrell, ibid. 28 (1956), 308. 

2) D. Bohm and D. Pines, Phys. rev. 92 (1953), 
609. 

3) L. Rosenfeld, Nuclear Forces (North-Holland 
Publishing Company, 1948), 77. 

4) J. M. Jauch and F. Rohrlich, The Theory of 
Photons and Electrons (Addison-Wesley Publish- 
ing Company, 1955), 286. 

5) K. Sawada et al, Phys. Rev. 108 (1957), 507. 
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